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CHAPTER 1 



Introduction 

1.1 Motivation 

Four types of force govern the universe, namely, Strong, Electromagnetic, Weak and Grav- 
itation. One of the most challenging and least understood problems in modern particle physics 
is the dynamics of the fundamental particles responsible for the strong interactions. Quantum 
Chromodynamics (QCD), the SU(3) color gauge theory provides the theory of strong interac- 
tions. But even after more than thirty years of the formulation of QCD, it is not well understood 
and not yet solved. The fundamental particles of this theory are called quarks. They can appear 
in six flavors up, down, strange, charm, bottom and top and in three colors, say, red, green and 
blue. The quanta of the non-Abelian gauge field in QCD, called gluon, mediates color interac- 
tions between quarks. Gluons are flavor blind but carry color charges and hence interact among 
themselves. This can be contrasted with the Abelian gauge theory, Quantum Electrodynam- 
ics (QED), which describes the theory of electromagnetic interactions. Photons which mediate 
electromagnetic interactions between charged particles in QED do not have any electric charge 
and hence cannot interact among themselves. This non-Abelian property of the gauge fields in 
QCD makes the game completely different from QED. Again as the interactions are "strong" in 
QCD (the relative strengths of the four types of interaction at hadronic scale are 1, 10~ 2 , 10~ 7 
and 10~ 39 ), the perturbative treatment which is very successful to solve QED, is not applicable 
to solve the theory for strong interaction. The most important and challenging part of QCD is 
to understand and solve the low energy or large distance physics where it shows "confinement" 
and is completely nonperturbative in nature as discussed below. 

In nature, quarks can only be found in colorless bound states, mesons and baryons which 
are combinedly called hadrons. The color symmetry imposed on QCD says that only colorless 
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states can be physically realized and hence a free quark which carries color charge cannot be 
observed. According to this postulate all hadrons are required to be in the singlet of color 
SU(3). It is just a kinematical constraint to eliminate non-singlet (colored) states. However, it 
is believed that the quark confinement is a dynamical consequence of QCD and thus it requires 
more attention and investigations. Owing to the non-Abelian nature of gluons QCD possesses 
an interesting property called "asymptotic freedom" [1] according to which the interactions 
between quarks become weak at very high energy (or at short distance). Thus, QCD enjoys 
asymptotic freedom at short distances while shows confinement at large distances (of the order 
of hadronic length scale). One can address the scattering experiments at very high energy (such 
as Deeply Inelastic Scattering (DIS), Deeply Virtual Compton Scattering (DVCS), etc) through 
perturbative QCD. But, investigation of low energy (or long distance) physics (such as hadronic 
bound state problem, nucleon-nucleon interactions), where the coupling strength become large, 
is not possible by perturbative theory and one needs a non-perturbative tool to compute the QCD 
bound state spectrum. It is also desirable to have informations about the hadronic wavefunctions 
which are essential to calculate the QCD observables such as form factors, structure functions, 
distribution functions, decay rates etc. We do not know how to exactly solve QCD. So, various 
approximations go into QCD calculations. But, no approximation is valid for all length scales of 
QCD. For example, perturbation theory which is applicable in high energy domain is not at all 
applicable for low energy calculations and low energy models assume dominance of few particle 
states. Again, the issues of confinement and hadronic bound states are much more difficult to 
be addressed with proper confidence due to the lack of adequate nonperturbative methods. To 
have control over realistic QCD calculations it is desirable to assess different nonperturbative 
methods in terms of their strengths and weaknesses. In this thesis, we will mainly concentrate 
on this aspect for different nonperturbative methods. 

Till date, the most practiced non-perturbative method is the lattice gauge theory [2]. But 
lattice gauge theory has its own difficulties and limitations. In this method, without fixing any 
particular gauge, one calculates different n-point Green functions by path integral formalism in 
Euclidean space, but can not have any firsthand knowledge about the bound state wavefunctions. 
To explore the structure of hadrons ( baryons or mesons) in terms of their constituent degrees 
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of freedom, a straightforward way is to solve the Hamiltonian eigenvalue equation 

H\ X ¥)=E\W) (1.1) 

expanding the eigenstate into multiparticle Fock states. The conventional Fock state expansion 
becomes intractable in equal-time framework (by equal-time framework we refer to the instant 
form of Hamiltonian dynamics with quantization surface given by x° = [3]) due to the com- 
plicated vacuum of the relativistic quantum field theory. Again in the equal-time framework the 
square root operator in E = \fP 2 +M 2 brings in severe mathematical difficulties. Even if one 
can solve these problems once, the eigensolutions are found in the rest frame, finding the solu- 
tions in a moving frame is highly nontrivial. This is because boosts are interaction dependent 
(dynamical) operators and boosting the system is as complicated as solving the theory. 

An elegant way to avoid all the difficulties mentioned above is to choose the light-front (also 
called light-cone) framework [3] where the quantization surface is chosen to be the tangential 
plane to the light-cone. The light-cone or light-front coordinates are defined as 

and the quantization surface is now chosen to be light-front time x + = instead of x° = in 
equal-time. The operator conjugate to x + is the light-front Hamiltonian P and P + which is 
conjugate to x~ is the light-front longitudinal momentum. (The detail of light-front coordinate 
system will be discussed in the next chapter.) The eigenvalue E in Eq. (1.1) does not involve 

j_2 2 

any square root operator and takes the form E = ^ , the vacuum structure is relatively sim- 
ple and the boosts are interaction independent (kinematical) operators. Thus, unlike the usual 
equal-time Hamiltonian formalism, it is possible to have a frame independent (boost invariant) 
description of the bound state wavefunctions in the light-front Hamiltonian approach. Since 
the vacuum state of the free Hamiltonian is also an eigenstate of the full QCD Hamiltonian, 
the Fock space expansion on this vacuum provides a complete relativistic many-body basis for 
diagonalization of the full QCD Hamiltonian. Again, as the many-body states are also high 
energy states (see the discussion on light-front in Chapter 2), one can hope to have few body 
description of the bound states and can reconcile relativistic field theory like QCD in one hand 
with constituent quark model (CQM) in the other hand. [The CQM was motivated by the hadron 
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spectroscopy. According to CQM, mesons are made of a quark and an antiquark and baryons 
are made of three quarks or three antiquarks bounded by some empirical (phenomenological) 
potential.] 

Two light-front non-perturbative methods that one can take up to address the QCD bound 
state problem are Similarity Renormalization Group (SRG) approach [4, 5] to construct an 
effective theory with few constituents and the light-front transverse lattice (LFTL) approach [6] 
which is a clever combination of lattice gauge theory and light-front Hamiltonian formalism. 
We will devote this thesis to investigate these two non-perturbative methods in the context of 
meson bound state problem. 

In (3+1) dimensions, Euclidean lattice gauge theory predicts linear confinement for quarks 
while the lowest order calculation in SRG scheme produces logarithmic confinement which 
again violates rotational symmetry. Automatically, questions come. Is the confinement pro- 
duced by SRG an artifact of the scheme? What type of confinement does it produce in (2+1) 
dimensions? Does that also violate rotational symmetry? In light-front field theory ultraviolet 
(UV) divergence issue is complicated due to different power counting rule on the light-front. 
In (2+1) dimensions this issue becomes simplified due to the absence of ultraviolet divergences 
except in mass corrections. Since in (2+1) dimensions QCD is superrenormalizable, the cou- 
pling does not run and one can keep the coupling arbitrarily small and study the structure of the 
bound states in a weakly coupled theory. In (2+1) dimensions one component of the gauge field 
remains dynamical and one can systematically study the effects of dynamical gluons without 
additional complications of (3+1) dimensions. One can also hope to enlarge the Fock space 
sector and investigate their effect on restoring Lorentz invariance. It is expected that such in- 
vestigations are more viable in (2+1) dimensions compared to (3+1) dimensions due to less 
severe demand on computational resources. A major part of this work is devoted to understand 
these issues in (2+1) dimensional QCD using SRG scheme. SRG is a modification over Bloch 
effective theory [7]. A study of meson bound state problem using Bloch effective Hamiltonian 
can serve as benchmark to assess the strengths and weaknesses of SRG approach. 

As we have already mentioned, another nonperturbative approach in light-front framework 
is the light-front transverse lattice formalism. In this approach, one keeps the light-front time 
(x + ) and longitudinal direction (x~) continuous while the transverse plane (jr 1 = (x l ,x 2 )) is 
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discretized on a square lattice. In the light-front field theories, UV divergences come only from 
small transverse separations. Lattice provides the gauge invariant UV cutoff on the transverse 
space. Thus, in the light-cone gauge A + = 0, one can still preserve x~ independent residual 
gauge invariance on the transverse plane. In this approach one avails the advantages of gauge 
invariant UV cutoff coming from the lattice on the transverse plane and the beautiful features of 
light-front framework. But, it is well known that formulation of fermions on a lattice is compli- 
cated due to the notorious problem of species doubling. In the usual (Euclidean or Minkowski) 
lattice gauge theory extra species of fermions are generated from the corners of the Brillouin 
zone. Is there any fermion doubling on light-front transverse lattice? If yes, do they also come 
from the corners of the Brillouin zone as usual lattice ? What are the possible ways to remove 
the doublers on the transverse lattice? In usual lattice gauge theory there are rigorous the- 
orems and anomaly arguments regarding fermion doubling. In standard lattice gauge theory, 
some chiral symmetry needs to be broken in the kinetic part of the action to avoid the doublers. 
In the light-front chirality is the same as helicity even for a massive fermion. The constraint 
equation for fermion in the light-front field theory violates the usual chiral symmetry. Here one 
should ask the question, is it still possible to relate light-front chirality and fermion doubling on 
the transverse lattice? Since the notion of light-front chirality is different from usual chirality, 
is there any way of formulating fermions on a light-front transverse lattice without generating 
extra species of fermions? There are many such questions one has to answer in order to use 
this method as a practical tool for QCD calculations. Our work shows that one can exploit 
the constraint equation for the fermionic field to formulate fermions in two different ways on 
the transverse lattice. In one formulation where we use forward and backward lattice deriva- 
tives without spoiling the Hermiticity of the Hamiltonian doublers do not appear. In the other 
way of formulating fermions people use symmetric lattice derivatives and encounters doublers. 
We have investigated both ways of formulating fermions on the transverse lattice and tried to 
understand the origin of doublers, possible ways to remove them and the symmetry relevant 
for doubling. Only when one understands the properties of fermions on the transverse lattice, 
realistic QCD calculations become viable. 



5 



Since we have two two possible ways of formulating fermions on a light-front transverse 
lattice (a) with forward and backward lattice derivatives and (b) with symmetric lattice deriva- 
tive, we need a comparative study of these two approaches in order to decide which one is best 
suited for QCD calculations. We take meson bound state problem in (3+1) dimensions for this 
purpose. Again we focus on the comparative analysis of the strengths and weaknesses of the 
different fermion formulations to deal with different QCD interactions rather than fitting data. 

One important issue in nonperturbative analysis is the numerical procedure. Since analytic 
solutions of the nonpertubative bound state equations are not possible, one needs to solve them 
numerically. In Hamiltonian formalism, we need to diagonalize the Hamiltonian numerically 
in a suitable basis. There are several numerical procedures to diagonalize a matrix. When one 
deals with a theory like QCD, several complicated interactions come with different ultraviolet 
and infrared singularities. It is very important to know the efficiency of the numerical proce- 
dures in handling singular interactions. In the meson bound state problem using SRG approach 
we use Gauss Quadrature method to convert the integral equations into a matrix eigenvalue 
problem and discuss the efficiency of this method in handling different infrared singular terms. 
Due to longitudinal dynamics, in the bound state equation for transverse lattice Hamiltonian, 
linear light-front infrared divergences and logarithmic infrared divergences in self energy dia- 
grams arise and one need to add counterterms to cancel them on the computer. We also discuss 
the efficiency of counterterms to cancel the divergences on a computer. 

1.2 Organization of the thesis 

The main objective of our work is to assess different nonperturbative approaches in the 
light-front formalism. It is very important to know the strengths and weaknesses of different 
approaches when one wants to do any nonperturbative QCD calculation. We investigate the me- 
son bound state problem in light-front QCD using different nonperturbative methods. Whenever 
it is possible, we make comparative studies of different approaches. We should again empha- 
size that our main aim is not to fit data but to assess the different nonperturbative approaches in 
terms of their strengths and weaknesses to have control over the QCD calculations. 

In Chapter 2, we discuss the basic features of light-front field theory which will heavily 
be used in this thesis. Then, we embark on detailed investigations of the meson bound state 
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problem in light-front QCD. We can broadly categorize the thesis into two parts, namely, (1) 
SRG approach and (2) transverse lattice formalism. 

In the SRG scheme, starting from a bare cutoff Hamiltonian and by some similarity transfor- 
mations one arrives at a low energy effective Hamiltonian which is band diagonal in Fock basis. 
The effective Hamiltonian is then diagonalized nonperturbatively. We have already mentioned 
that SRG is a modification over Bloch effective theory and a study of Bloch theory can serve 
as a benchmark to assess the SRG approach. In Chapter 3, we investigate the meson bound 
state problem with Bloch effective Hamiltonian in (2+1) dimensions. We start with a brief dis- 
cussion on the basic ideas about the effective field theory (EFT). We observe that Bloch theory 
is infected with infrared divergences in (2+1) dimensions and as one cranks up the strength of 
the coupling constant eigenvalues diverge bringing instability in the system. Thus we define a 
reduced model which is free from the divergences but still has confinement in the lowest non- 
trivial order of the expansion in coupling constant. This allows us the opportunity to study the 
manifestation and possible violation of rotational symmetry in the context of light-front field 
theory. 

After having the knowledge of difficulties and shortcomings of Bloch effective theory in 
the context of (2+1) dimensional meson bound state problem, we investigate the same meson 
bound state problem using similarity renormalization approach in Chapter 4. After a brief re- 
view of renormalization approach of constructing a low energy effective theory, we move on to 
the detailed analysis of SRG scheme. We immediately see the improvements due to SRG over 
Bloch effective theory. The bound state integral equation is converted into a matrix diagonaliza- 
tion problem by Gauss Quadrature method. The Gauss Quadrature method is quite efficient in 
handling divergences. We find that SRG produces linear confinement in transverse direction but 
only square root potential in the longitudinal direction. Thus, it severely violates the rotational 
symmetry. Our results show that higher order calculations are essential to investigate whether 
rotational invariance can be restored or not. By tuning the similarity cutoff and the quark mass 
we also study the interplay between SRG generated confining interaction and the Coulomb 
interaction which in (2+1) dimensions gives rotationally invariant logarithmic confinement. 

In Chapter 5, we introduce the other promising nonperturbative approach in the light-front 
framework, namely, light- front transverse lattice (LFTL). We have already emphasized the 
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problem of fermion doubling in usual lattice gauge theory and the relevant questions on the 
transverse lattice one should ask. In this chapter, we propose a new method of putting fermions 
on a light-front transverse lattice which is free from fermion doubling. In this approach we use 
forward and backward lattice derivative in such a way that the Hermiticity of the action is not 
spoiled. In that case an irrelevant helicity flip interaction survives in the free field limit. We 
also discuss the violation of rotational symmetry on the transverse lattice. Then we discuss the 
transverse lattice Hamiltonian with symmetric lattice derivative where one encounters doublers. 
Our results show that the origin of doublers on the transverse lattice is not the same as usual 
lattice gauge theory. Here, the doublers appear due to decoupling of even and odd sub-lattices. 
The removal of doublers by two different ways, namely staggered fermion formulation and Wil- 
son fermion formulation are discussed. We identify a even-odd helicity flip symmetry on the 
transverse lattice relevant for fermion doubling. 

Once we understand the properties of fermions, in Chapter 6 we investigate the meson bound 
state problem in (3+1) dimensions with fermions formulated with forward and backward lattice 
derivatives and fermions formulated with symmetric lattice derivative. In the case of symmetric 
lattice derivative, we add a Wilson term to remove the doublers. In our investigation we use 
one link approximation, i.e., quark and antiquark at most can sit one lattice spacing apart on the 
transverse plane. We compare and contrast the two ways of fermion formulation in the context 
of meson bound state problem. The major difference between the two approaches is that, with 
forward and backward lattice derivatives, hopping of quark (or antiquark) in the transverse plane 
with helicity flip interferes with helicity non flip hopping, while there is no such interference 
with symmetric lattice derivatives. The consequence of this interference in the spectrum is also 
studied. 

Summary and conclusions are given in Chapter 7. To ease the reading, we provide refer- 
ences at the end of each chapter. Several appendices are provided to clarify the notations and 
formalisms used and also to elucidate different intermediate steps. 
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CHAPTER 2 



Some Basic Features of Light-Front Field Theory 



In this Chapter, we introduce the basic features of light-front field theory, in the context of 
light- front QCD (LFQCD). Since in the next few chapters we will heavily use the light-front 
coordinates and the features of field theory in this formalism, it is instructive to discuss its ad- 
vantages and the problems that one must understand for a successful practical application of 
the theory. As, in Chapter 1, we have already emphasized the advantages of using Hamiltonian 
formalism to investigate the hadronic bound state problem and as we will perform all our in- 
vestigations in the light-front Hamiltonian approach, here we shall only be concerned with the 
Hamiltonian formulation of LFQCD. 

In 1949, Dirac [1] showed that there are three independent parametrizations of the space 
and time that can not be mapped on each other by a Lorentz transformations and discussed 
three forms of Hamiltonian dynamics. In the equal-time Hamiltonian formulation of field the- 
ory, quantization conditions in the form of commutator (or anticommutaor) of dynamical fields 
and their conjugate momenta are specified on the space-like hypersurface x° = and the Hamil- 
tonian generates the time-evolution of the system (Dirac called it instant form as the kinematical 
part of the Lorentz group leaves the instant invariant). In the front form, the quantization con- 
ditions are specified on a light-like hypersurface x + = x° + x 3 = (called a light- front) and the 
light-front Hamiltonian generates the evolution for a new time (x + ). This formulation is known 
as the light-front Hamiltonian field theory. Another form that Dirac mentioned is the point 
form of Hamiltonian dynamics where the quantization hypersurface is given by the hyperboloid 
x^Xj_i = K 2 with x° > and K 2 > 0, and the Lorentz group leaves a point invariant. However, 
later on two more possibilities of parametrization of the space and time were found [2]. The 
quantization hypersurfaces for these two parametrizations are given by x\ — x\ — x\ = K 2 > 
with xo> and Xq — x 2 = K 2 > with xq > 0. Among all the parametrizations, the front form 
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has the largest stability group, the subgroup of the Poincare group that maps the quantization 
hypersurface onto itself. 

There is no well defined guideline to decide which parametrization one should use. High 
energy experiments (e.g., deep inelastic scattering) probe the hadrons near the light-cone. It 
motivates people to use light-front parametrization of space and time to explore the QCD ob- 
servables. One may hope that for highly relativistic systems in which cases the world-line lies 
very close to the light-cone, physics will be more transparent and it will be relatively easy to 
extract them if one uses light-front field theory. 

In the context of current algebra, Fubini and Furlan [3] introduced another notion of Lorentz 
frame known as Infinite-Momentum Frame (IMF) as a limit of a reference frame moving with 
almost the speed of light. Weinberg [4] using old-fashioned perturbation theory for scalar me- 
son showed that vacuum structures become simplified in the infinite-momentum limit. Later, 
Susskind [5] established that although the Lorentz transformation required to arrive at IMF is 
evidently singular (7 = 1 / \J 1 — ^ — > 00 as v — > c), the singularity cancels in the calculation of 
physical objects (like Poincare generators) and results in an effective coordinate change given 
by 

x ± =x°±x 3 , x ± = {x l ,x 2 }, (2.1) 

same as the light-front coordinate we defined in Chapter 1. Thus, one can see the fact that what 
one obtains after going through singular limiting procedure in IMF is built in quite naturally in 
the light front field theory. That is why, light-front field theories are also sometimes referred as 
field theories in the infinite-momentum frame. But, we should reemphasize that the formulation 
here is as prescribed by Dirac and has no connection with any singular limiting procedure. For 
a review and exhaustive list of references on light-front field theories see Ref. [6]. 

2.0.1 LF dispersion relation 

The inner product between two four- vectors is defined on the light front as 

1 1 

*-y= 2 x+y ~ + 2 x ~ y+ ~ x± ' y± - ( 1 - T> 

In analogy with the light-front space-time variables, the light-front four momenta are defined as 

k ± =k°±k\ k- L = {k\k 2 }, (2.3) 
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where kr being conjugate to x + is the light front energy and k + which is conjugate to x~ is the 
light-front longitudinal momentum. With the above definitions, the dispersion relation, i.e., the 
relation between light- front energy kr and the spatial components of momenta (k + ,k. ± ), for an 
on mass-shell particle of mass m, is given by, 

One of the remarkable features of this relativistic dispersion relation is that there is no square 
root involved in contrast to the relativistic equal-time dispersion relation E = \fP + m 2 . This 
provide great simplification when one tries to solve eigenvalue equation which we have already 
emphasized in chapter 1 . Secondly, the numerator in Eq. (2.4) being always positive implies that 
the particles with positive light-front energy (k~) always carry positive longitudinal momentum 
(k + > 0). As usual, the particles with negative k~ which must have negative k + are mapped 
to antiparticles with positive k~ and k + . As a consequence, we always have k + > for real 
particles. Thirdly, kr becomes large for the large value of k 1 - as well as very small values of 
k + . This makes light front renormalization aspects very different from the usual one. Lastly, the 
dependence on the transverse momenta is just like a nonrelativistic dispersion relation. We 
shall see later in this Chapter the crucial implications of these novel features of the light-front 
dispersion relation in the light-front field theory. 

2.0.2 The light-front vacuum 

The above dispersion relation has profound consequence in the vacuum structure of light- 
front field theory. Vacuum state is always an eigenstate of the longitudinal momentum P + \ 
0) = 0. The positivity condition of k + (k + > 0) implies that the vacuum | 0) is either a no 
particle state or, at most can have particles with longitudinal momenta exactly equal to zero. 
Now, if we consider a cut-off theory where longitudinal momentum is restricted to be k + > e, 
the vacuum state | 0) becomes completely devoid of any particle and therefore, an eigenstate 
of the full interacting Hamiltonian with zero eigenvalue. Thus, the light-front vacuum becomes 
trivial. It should be contrasted with equal-time case where the vacuum has highly complicated 
structure. In equal-time case, vacuum can contain infinite number of particles moving with 
positive and negative momenta adding up to zero. Another aspect of the cutoff k + > £ is that 
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it automatically puts a restriction on the number of constituent particles a state with finite P + 
can have. A composite state with total longitudinal momentum P + now can have at most P + /e 
constituents. This again simplifies the Fock space expansion for the hadronic bound states. 
Again, as small k + means high energy (large k~), one can hope to have a few body description 
for the low lying hadron states and reconcile QCD with CQM which is beyond hope in equal- 
time formalism. 

On the other hand, a complicated vacuum structure is supposed to be responsible for spon- 
taneous chiral symmetry breaking or confinement in QCD. It seems that with the trivial vacuum 
structure in light- front theory after removing the zero modes (k + = 0), we may lose these im- 
portant aspects in our theory. It should be emphasized that we have not simply removed the zero 
modes from our theory. The longitudinal momentum cut-off (g) should be removed from the 
theory at the end of any calculation by adding necessary counter terms in the effective Hamil- 
tonian to render the observables independent of £. Thus, we expect to get back all the effects of 
zero mode as an effective interaction in the Hamiltonian through renormalization. 

2.0.3 Poincare generators in light-front 

In equal-time theory, out of the ten Poincare generators (Hamiltonian (P°), three linear 
momenta (P), three angular momenta (/) and three boosts (K)) six are kinematical {P,/}, i.e., 
they do not depend on the dynamics (interactions) and other four are dynamical {P°,K}. 

In light-front, Poincare generators can be constructed in the same way as in equal-time 
case. Starting from Lagrangian density we construct the energy momentum stress tensor T^ v 
and from the stress tensor we construct the four-momentum P^ and the generalized angular 
momentum M^ v defined in the following way. 



In light-front dynamics P is the Hamiltonian and P + and P ! with (/ = 1 , 2) are the longitudinal 
and transverse momenta. M + ~ = 2K 3 and M +l = E l axe the boost operators and M 12 = J 3 and 
M~' = F' are generators for rotations. In light-front theory, boost operators (^ 3 and E l ) are 
kinematical. Longitudinal boost is like a scale transformation and the transverse boosts behave 




(2.6) 



(2.5) 
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like Gallilean boosts in the nonrelativistic theory. To elucidate it further let us consider the boost 
along the 3-axis (K 3 ) as an example. In equal-time, K 3 transforms the time and the 3-axis 
(x 3 ) but leaves the transverse space invariant. 

x° = y(x°-px 3 ), x 3 = 7(x 3 -j8x°), x l > 2 =x 1 > 2 , (2.7) 

where j8 = ^ and 7 = -j= ^ . From the above equations we see that K 3 changes the quantization 
surface x° = and hence, K 3 is a dynamical generator in equal-time theory. Introducing the 
parameter § such that 7 = cosh0 and Py= sinh0, we see that, in the light- front 

S*=&+Z i = e-+x + , r =x°-x 3 = e* x~, x l > 2 =x l > 2 . (2.8) 

It clearly shows that K 3 , which is known as generator of longitudinal boost in light-front, be- 
haves like a scale transformation. In particular, it keeps the quantization surface x + = in- 
variant. Therefore, it is a kinematical generator in light-front theory. On the other hand, two 
rotations about transverse axes (F 1 and F 2 ) which are kinematical in equal-time case become 
dynamical in light-front theory. Thus, in the light-front theory, we have seven kinemetcal ( 3 
boosts, 3 translations and rotation about 3-axis), and three dynamical (Hamiltonian and two 
rotations about the transverse axes) generators. 

Notice that the boost generators form a closed algebra among themselves which is similar 
to the generators of non-relativistic dynamics: 

[E\E 2 ]=0, [K 3 ,E i ]=iE\ (2.9) 

and 

[j 3 ,E t \=ie iJ E j , (2.10) 

where £ 12 = — e 21 = 1 and £ H = e 22 = 0. Also, F 1 , F 2 and J 3 form a closed algebra. 

[F 1 ,F 2 ]=0, [j 3 ,F i ]=ie ij Fi. (2.11) 

For more details of the Poincare algebra in light- front see Ref. [6]. 

Since the kinematical subgroup of the Poincare group enlarges and contains seven genera- 
tors in light-front theory, it is expected that defining a system will be easier in the light-front 
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theory as we can fix more variables of the system irrespective of its dynamics. Moreover, since 
different set of generators are kinematical in light- front compared to the equal-time theory, it is 
worth pursuing this theory, for certain things difficult to study in equal-time may just become 
simpler here. One such example is the feasibility of representing the QCD-bound states in terms 
of just a few boost invariant multi-particle wave-functions in the Fock-space expansion, which 
we discuss next. 

2.0.4 Basic strategy for bound state problem 

The starting point is the Hamiltonian eigenvalue equation 

p 2 i y/2 

P-\^)= p |y> (2-12) 

where is the invariant mass-squared of the state | As it is mentioned already, trivial 
structure of light-front vacuum makes it feasible to study the hadronic bound states in Fock 
language. Since the Fock-states form a complete basis, any state vector, in principle, can be 
expanded in terms of that basis introducing corresponding amplitude for each Fock-basis. The 
bound state of a hadron on light-front can be simply expanded in terms of the Fock states as 

I W ) = £ / dx i d2K -Li I n ^ Xi P ± + K U X) &n(Xi, KuM , (2-13) 
n.Xi 

where n represents n constituents contained in the Fock state | n,Xj,XiP± + Kj_j,A,-), is the 
helicity of the z'-th constituent, X[ is the fraction of the total longitudinal momentum carried by 
the z'-th constituent, and JCy is its relative transverse momentum with respect to the center of 
mass frame. 

p+ 

Xi = pT ' *i±=Pi±- x i p ±i (2.14) 
with pf, p t ± being the longitudinal and transverse momenta of the z'-th constituent 



J> = 1, and 2> ±J - = (2-15) 

and /' denotes the integral over the space. 4>„(^-, Kj_,-, Aj) is the amplitude of the Fock state 
| n,xi,XiP^ + Kj_,-, Xi), i.e., the multi-parton wave function, which is boost invariant and satisfies 
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the normalization condition: 



£ J dxid 2 K ±i \<f> n (xi,K ±i ,Xi)\ 2 = l. 

n.Xi 



(2.16) 



For example, if we consider the meson bound state problem, then after expanding the eigenstate 
in the basis of Fock states, the light- front bound state equation can be written as, 

<J>, 



v i=\ x i y 



4> 
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(qq\H int \qq) (qq\H int \qqg) 
(qqg\H int \qq) 



4>, 
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(2.17) 



Here H int is the interaction part of the light-front QCD Hamiltonian. 

The expansion has infinite number of terms and it is impossible to solve the bound state 
equation, Eq. (2.17), which is an infinite dimensional coupled equation. To make any practi- 
cal calculation viable using Fock-expansion, one needs to truncate the expansion at a suitable 
maximum particle number with the hope that a first few terms in the expansion may give useful 
information. We know two important informations from the light-front dispersion relation, one 
is that longitudinal momentum p + is always positive and secondly, states with small p + are 
high energy states (large p~). Since p + is always positive, constituents of a many particles state 
with a fixed longitudinal momentum P + carry only small amounts of p ; + and hence the state is 
of high energy. Since high energy states are weakly coupled one can hope that the dynamics of 
the bound states is dominated by few particle states and multiparticle states can be considered 
in a bound state perturbation theory in a consistent manner. 

The situation should again be contrasted with equal-time approach. As each Fock-state is 
obtained by operating various creation operator(s) on the vacuum of the theory, if the vacuum 
already has a complicated structure (as is the case in equal-time theory), which may contain 
arbitrary number of particles and thereby, the vacuum itself needs a Fock-expansion. This, in 
effect, render the Fock-expansion in equal-time theory meaningless for any practical applica- 
tion. This is not the case in light-front theory due to the simplicity of the vacuum. 

2.0.5 Renormalization aspects 

In light-front field theory in the Hamiltonian framework, the renormalization is a more com- 
plicated issue mainly due to the noncovariant structure of the theory and is quite different com- 
pared to the usual co variant one. This is due to the fact that the power counting in light front 
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is very different. For a detailed discussion on light- front power counting, see the Ref. [7]. 
Here we notice the fact that only transverse directions x 1 - carry the mass dimension, while the 
longitudinal direction x~ has no mass dimension. Thus, one has to treat transverse and longi- 
tudinal directions separately in determining the superficial degree of divergence of a divergent 
integral by power counting, in contrast to the covariant case where all the space-time direc- 
tions are treated democratically. This is also evident in the single particle dispersion relation 
k = ^ — |+ — , which shows that there are two sources of divergences: k + — > 0+ and k — > oo. 
The divergence coming from k + — > 0+ is referred as infrared (IR) divergence, whereas k 1 - — > °° 
is known as the ultraviolate divergence (UV) in light-front theory. 

For the above reason, dimensional regularization, which is so elegant and commonly used in 
covariant perturbation theory, is of very little importance in light- front theory. Only in the trans- 
verse direction, one may use dimensional regularization. Since we know that the lattice gauge 
theory provides gauge invariant UV cutoff, another way to regulate the UV divergences is to 
discretize the transverse plane on a square lattice [8]. IR divergences are also regularized by 
putting a small longitudinal momentum cut-off, which is equivalent to using principal value pre- 
scription for the integration over longitudinal momenta. Also the fact that the light-front theory 
being gauge fixed and noncovariant, leads to new type of divergences like quadratic divergences 
(if we are using cut-off instead of transverse dimensional regularization) in mass renormaliza- 
tion or mixed divergences involving both IR and UV ones. To remove these divergences one has 
to add counter terms to the canonical Hamiltonian, which are often nonlocal and help to restore 
the invariance of the theory that might be broken in the process of manipulation. For detailed 
discussion on this subject, see the Refs. [7, 9, 10, 11]. Another method specially designed to 
address the bound state problem in light-front, is that of similarity renormalization introduced 
by Glazek and Wilson [12] and Wegner [13], where first an effective Hamiltonian is obtained 
perturbatively, by performing a similarity transformation on the bare UV cutoff Hamiltonian. 
Similarity renormalization approach will be discussed in detail in Chapter 4 and in Appendix 
E. Notations and conventions are given in Appendix A. 



17 



BIBLIOGRAPHY 

[1] P. A.M. Dirac, Rev. Mod. Phys. 21, 392(1949). 
[2] H. Leutwyler and J. Stern, Ann. Phys. 112, 94 (1978). 
[3] S. FubiniandG. Furlan, Physics, 1, 229 (1965). 
[4] S. Weinberg, Phys. Rev. 150, 1313 (1966). 
[5] L. Susskind, Phys. Rev. 165, 1535 (1968). 

[6] A. Harindranath, An Introduction to Light-Front Dynamics for Pedestrians in Light-Front 
Quantization and Non-Perturbative QCD, edited by J.P Vary and F. Wolz (IITAP, Ames, 
Iowa, USA, 1997), hep-ph/96 12244. 

An exhaustive reference of light-front literature can also be found at 
http://theory.saha.ernet.in/~hari/light/light.html. 

[7] K.G. Wilson, T.S. Walhout, A. Harindranath, W.M. Zhang, R.J. Perry, and St.D. Glazek , 
Phys. Rev. D 49, 6720 (1994); 

[8] W. A. Bardeen and R. B. Pearson, Phys. Rev. D 14, 547 (1976); W. A. Bardeen, R. B. Pear- 
son and E. Rabinovici, Phys. Rev. D 21, 1037 (1980). 

[9] W.M. Zhang and A. Harindranath, Phys. Rev. D 48, 4881 (1993). 

[10] A. Harindranath and W.M. Zhang, Phys. Rev. D 48, 4903 (1993). 

[11] R.J. Perry, Ann. Phys., 232, 116 (1994). 

[12] S. D. Glazek and K.G. Wilson, Phys. Rev. D 48, 5863 (1993); D 49, 4214 (1994). 

[13] F. Wegner, Ann. Phys. (Leipzig) 3, 77 (1994). 



18 



CHAPTER 3 



Bloch Effective Hamiltonian and Bound State Problem in 
(2 + 1) -Dimensional Light-Front QCD 

3.1 Introduction 

One of the very well known techniques to extract relativistic bound state solutions is the 
Bethe-Salpeter formalism [1]. Though Bethe-Salpeter equation is formally an exact equation 
for bound state problem, provides a covariant formalism and successful in quantitative under- 
standing of bound states in different models [2] and positronium bound states in QED [3], the 
calculations in this approach are very complicated and almost out of control beyond the ladder 
approximation (for a review in the context of QCD see [4]). 

The straightforward way to extract the relativistic and nonperturbative wavefunctions is the 
Hamiltonian approach where one solves the eigenvalue equation //(^F) =£'| V F). But the straight- 
forward diagonalization of the Hamiltonian has two major problems namely (1 ) it involves in- 
finitely many energy scales and (2) the rapid growth of the dimension of the Hamiltonian matrix 
with particle number. In the spirit to diagonalize the Hamiltonian in a single step, one may 
implement Discretized Light Cone Quantization (DLCQ) [5]. DLCQ has been quite successful 
in two dimensional models, but for QCD this approach may be quite ambitious. Typical Hamil- 
tonians of interest couple low energy scales with high energy scales which results in ultraviolet 
divergences. Furthermore, Hamiltonian couples every particle number sector allowed by sym- 
metries and at strong coupling, brute force particle truncation can fail miserably. An alternative 
approach will be to use an effective Hamiltonian that operates in a few particle basis. 

Effective field theory (EFT) [6] relies on the assumption that physics at a low energy scale 
is insensitive to the microscopic details of the underlying physics at a high energy scale. EFT 
provides a powerful framework to study low-energy phenomena where one can replace the 
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microscopic degrees of freedom and their interactions by effective macroscopic degrees of free- 
dom and their effective interactions. The basic procedure of EFT is to separate out the important 
field components and redefine the theory with those fields within a certain range of energy and 
momentum, so that, acting on a limited Hilbert space it produces the same result as the original 
theory. In other words, to construct a low energy effective theory one needs to "integrate out" 
the high energy degrees of freedom (degrees of freedom above the scale considered) from the 
theory. The effective interactions are renormalized accordingly to incorporate the effects of the 
degrees of freedom above the cutoff. An effective theory describes the main features of the 
original theory below the scale one considers in a simpler way. Effective field theory has wide 
applications in different wings of physics such as condense matter physics, nuclear physics, 
high energy physics, etc. 

The strength or validity of an effective theory depends on how accurately the effective in- 
teractions mimic the effects of the degrees of freedom thrown out from the theory. There are 
many approaches to construct the effective Hamiltonian which acts on few particle states. It is 
also well known that there are some or the other drawbacks in all effective Hamiltonians. As 
the light-front framework is very much suitable for Hamiltonian formalism and due to triviality 
of the vacuum one can expand the bound states in Fock basis states, several attempts have been 
made for nonperturbative diagonalization of the light-front Hamiltonian for relativistic bound 
states (for a review see, Ref. [5]). One of the first attempt was to implement the Tamm-Dancoff 
truncation [7] or B loch-Horowitz effective Hamiltonian [8]. Though Tamm-Dancoff was suc- 
cessful in tackling (1+1) dimensional gauge theories, its deficiencies become apparent when 
attempts were made in (3+1) dimensions. First and foremost is the lack of confinement in the 
case of QCD in the first non-trivial order. Second is the appearance of the bound state eigen- 
value in the energy denominators. This has two undesirable consequences. Firstly, a light-front 
singularity of the type Ar, where k + is the light front longitudinal momentum of the exchanged 
gluon, remains in the bound state equation, which would have canceled if free energies appeared 
in the energy denominators. Secondly, for example, consider the meson bound state problem 
and truncate the Fock space with qq and qqg states. From the fermion self energy contribution 
(Fig. 3. 1(a)), in addition to the mass divergence another ultraviolet divergence appears (for an 
example in the context of (3+1) dimensional Yukawa model see Ref. [9]) which contributes 
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to the renormalization of the coupling. This contribution is also infrared divergent and can be 
identified as arising from fermion wave function renormalization (Z2). It is the Fock space trun- 
cation that has produced this unphysical divergence which would otherwise have been canceled 
by vertex renormalization (Zi) (Fig. 3. 1(b) and (c)) in a strict order by order perturbative calcu- 
lation. Thus, it severely violates the gauge invariance (gauge invariance demands Z\ = Zj) and 
one has to abandon the Tamm-Dancoff formalism in more than ( 1 + 1 ) dimensions. Another ap- 




(a) (b) (c) 

Figure 3.1: x + -ordered Hamiltonian diagrams for (a) self energy, (b) and (c) vertex correction. 
Diagrams (b) and (c) are not allowed by Tamm-Dancoff truncation as they involve two gluons 
in the intermediate state (shown by dashed lines). 

proach mentioned above is the Bloch-Horowitz effective Hamiltonian approach. But one of the 
major drawbacks of Bloch-Horowitz formula for effective Hamiltonian is that the transforma- 
tion rule does not preserve the ortho-normalization condition of the wavefunctions which is very 
much important for observable calculations. Bloch effective Hamiltonian [10] (also reinvented 
by Wilson [11] in the context of renormalization group) is a modification of Bloch-Horowitz 
Hamiltonian and some of the deficiencies of the Bloch-Horowitz formalism are absent in the 
Bloch effective Hamiltonian ( for a review see [12]). 

Bloch Hamiltonian has two desired properties, namely, the effective Hamiltonian is (1 ) Her- 
mitian and (2) involves only unperturbed energies in the energy denominator. The basic for- 
malism of calculating Bloch effective Hamiltonian is discussed in Appendix B. Use of Bloch 
effective Hamiltonian eliminates two major problems of the Tamm-Dancoff approach to gauge 
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theories mentioned above. However, Bloch effective Hamiltonian involves the undesirable van- 
ishing energy denominators. To the best of our knowledge, Bloch effective Hamiltonian was 
never assessed in terms of its strengths and weaknesses in the study of bound state problems in 
field theory. We first explore the Bloch effective Hamiltonian in the context of meson bound 
states in (2 + 1) dimensional light-front QCD. 

A major feature of gauge theories on the light-front is severe light-front infrared divergence 
of the type where k + is the exchanged gluon longitudinal momentum which appears 
in instantaneous four-fermion, two-fermion two-gluon, and four-gluon interactions. In old- 
fashioned perturbation theory these divergences are canceled by transverse gluon interactions. 
In similarity perturbation theory [13] which is considered in the next Chapter, the cancellation is 
only partial and singular interactions survive. Before embarking on a detailed study of effective 
Hamiltonian in the similarity renormalization approach which is a modification of the Bloch 
effective Hamiltonian, it is quite instructive to study the Bloch effective Hamiltonian itself. The 
result of such a study serves as benchmark against which one can evaluate the merits of simi- 
larity renormalization scheme. This also provides us quantitative measures on the strengths and 
weaknesses of numerical procedures in handling singular interactions (in the context of light- 
front field theory) on the computer. It is crucial to have such quantitative measures in order to 
study the effects of similarity cutoff factors on the nature of the spectrum. 

Just as the Tamm-Dancoff or the Bloch-Horowitz formalism, Bloch effective Hamiltonian of 
QCD in the first non-trivial order also does not exhibit confinement in (3+1) dimensions. Since 
one of our major concerns is the study of spectra for confining interactions, we go to (2+1) 
dimensions. In this case, in the limit of heavy fermion mass, a logarithmic confining potential 
emerges. There are several other reasons also to study light- front QCD in (2+1) dimensions. 
They arise from both theoretical and computational issues which we discuss next. 

First of all, issues related to ultraviolet divergence become more complicated in the light- 
front approach since power counting is different [14] on the light front. We get products of 
ultraviolet and infrared divergent factors which complicate the renormalization problem. Going 
to two space one time dimensions greatly simplifies this issue due to the absence of ultraviolet 
divergences except in mass corrections. An extra complication is that Fock space truncation in- 
troduces extra ultraviolet divergences which complicate the situation in non-perturbative bound 
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state computations [9]. Such special divergences do not occur in (2+1) dimensions. A third 
complication one faces in (3+1) dimensions is that on enlarging the Fock space in a bound state 
calculation, one soon faces the running of the coupling constant. At low energy scales, the 
effective coupling grows resulting in a strongly coupled theory [15] making the weak coupling 
approach with a perturbatively determined Hamiltonian unsuitable or making it mandatory to 
invent mechanisms like non-zero gluon mass to stop the drastic growth [14]. In (2+1) dimen- 
sional QCD we do not face this problem since the coupling constant is dimensionful in this 
superrenormalizable field theory and does not run due to ultraviolet divergence. We can keep 
the coupling arbitrarily small and study the structure of the bound states in a weakly coupled 
theory. 

Secondly, in (1+1) dimensions, in the gauge A + = 0, dynamical gluons are absent and their 
effect is felt only through instantaneous interactions between fermions. Further, recall that in 
light front theory, vacuum is trivial. As a result, the Fock space structure of the bound states 
are remarkably simple. For example, the ground state meson is just a qq pair both at weak and 
strong couplings. In contrast, in (2+1) dimensions, one component of the gauge field remains 
dynamical and one can systematically study the effects of dynamical gluons. Also note that 
(2+1) dimensions are the lowest dimensions where glueball states are possible and offer an 
opportunity to study their structure in the Fock space language without additional complications 
of (3+1) dimensions. 

A third reason deals with aspects of rotational symmetry. (2+1) dimensions offer the first op- 
portunity to investigate violations of Lorentz invariance introduced by various cutoffs (momenta 
and/or particle number) in the context of bound state calculations. This is to be contrasted with 
(1+1) dimensions where the sole Lorentz generator, namely boost, is kinematical in light-front 
field theory. Since in (2+1) dimensions we have a superrenormalizable field theory, violations 
introduced by transverse momentum cutoffs are minimal. Thus in contrast to (3+1) dimensions, 
one can study the violations caused by truncation of particle number alone and longitudinal 
momentum cutoffs. It is also conceivable that one can enlarge the Fock space sector and in- 
vestigate their effect on restoring Lorentz invariance. It is expected that such investigations are 
more viable in (2+1) dimensions compared to (3+1) dimensions due to less severe demand on 
computational resources. 
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A fourth reason concerns similarity renormalization approach. In (3+1) dimensions it has 
been shown that similarity renormalization group approach [13] to effective Hamiltonian in 
QCD leads to logarithmic confining interaction [16]. It is of interest to investigate corresponding 
effective Hamiltonian in (2+1) dimensions especially since the canonical Hamiltonian already 
leads to logarithmic confinement in the nonrelativistic limit in this case. It is also known that in 
(3+1) dimensions the confining part of the effective Hamiltonian violates rotational symmetry. 
Does the violation of rotational symmetry occur also in (2+1) dimensions? If so, how does it 
manifest itself? 

In [17] we initiated a systematic study of light-front QCD in (2+1) dimensions to investigate 
the various issues discussed above. 

3.2 Canonical Hamiltonian 

In this section we present the canonical light front Hamiltonian of (2+1) dimensional QCD. 
The Lagrangian density is given by 

^=[-\(hca) 2 + W(y X (idx+gAx)-m)v] (3.1) 

with 

pV-Xa _ ^^pXa _ fiX pjla _|_ gfabc j^Jlbj^Xc (?>2) 

We have the equations of motion, 



ifd^+gfA^-m 



y = 0, (3.3) 



d^ Va +gf abc A^ v +gfr v T a W = 0. (3.4) 

Because we are in (2+1) dimensions, we immediately face an ambiguity since there are no y 
matrices in (2+1) dimensions. In the literature both two component [18] and four component 
representation [19] have been in use. For simplicity, we pick the two component representation. 
Explicitly, 
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^ = r°±r 2 , r + = 



o o 

2i 



r = 



-2i 




a*.^. a + ^(; s).a-.(|» ;). 



Fermion field operator i/z^ = y. We have 



vr 



o 



(3.6) 



(3.7) 



(3.8) 



where ^ and 77 are one component fields. We choose the light front gauge A +a = 0. From the 
equation of motion, we get the equation of constraint for fermion 



id + y = a 1 (id 1 + gA l ) + fm 



Thus the fermion constrained field 



d+ 



(id 1 +gA l ) + im 



1- 



The equation of constraint for the gauge fields is 



-(d + ) 2 A- a = -d l d+A la -gf lbc A lb d + A lc -2grfT a r). 



(3.9) 



(3.10) 



(3.11) 



Using these equations of constraint, we eliminate i// and A in favor of dynamical field \\f + 
and A 1 and arrive at the canonical Hamiltonian given by 



H = H + H int = J dx-dx l (M + J^int). 

The free Hamiltonian density is given by 

^ = f -^f + m 2 + _Ul A l fl)9 l A l fl 
w + 2 

The interaction Hamiltonian density is given by 



(3.12) 



(3.13) 



(3.14) 



with 



g^A l —r]+g^^-(A l j] 



tl 1 

gmrfA l — T) +gmri i — (A l ri) 
2g-^ T (d l A la )ri J 'T a ri+gf abc d l A la -^ T (A lb d + A lc ) 
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(3.15) 



and 

+2g 2 f abc 4-(ri J( T a ri)4-(A lh d + A lc ) 

a + o + 

+^g 2 f abc f ade -^(A lb d+A lc )-^ : (A ld d + A le ). (3.16) 

The one component fermion field is given by 

r 7 (x + = 0,x-,x 1 ) = / f \f ■ \b{k)e-** + <ft{k)M ■ (3.17) 
7 2(2^) 2 VA: + L J 



2(2^) 2 v^+ 

The Fock operators obey the anti commutation relation 

{b(k), b\q)} = 2(2K) 2 k + 8 2 (k-q), {d(k), d\q)} = 2{2n) 2 k + 8 2 {k - q) , (3.18) 

other anti commutators being zero. Note that in two component representation, light front 
fermions do not carry helicity in (2+1) dimensions. 

In free field theory, the equation of motion of the dynamical field A 1 is the same as that of a 
free massless scalar field [20] and hence we can write 

AV = 0,x-,x l ) = J ^f k+ [a{k)e- ihx + a\k)e- ihx \ (3.19) 

The Fock operators obey the commutation relation 

[a{k),a\q)} = 2(2k) 2 k+ 8 2 (k - q), (3.20) 

other commutators being zero. 

We substitute the Fock expansions, Eqs. (3.17) and (3.19) into the Hamiltonian and treat all 
the terms to be normal ordered. Thus we arrive at the canonical Hamiltonian in the Fock basis. 

3.3 Bloch effective Hamiltonian in the meson sector and the bound state 
equation 

In this section we evaluate the Block effective Hamiltonian to the lowest non-trivial order 
for a meson state and derive the effective bound state equation. We define the P space to be qq 
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sector of the Fock space and Q space to be the rest of the space (see Appendix B for details). In 
the lowest non-trivial order, the Bloch effective Hamiltonian is given by 



{a\H eff \b) = (a\(H + H int )\b) 



+ -Y,(a\H int \k}(k\H int \b) 



1 



— + ■ 

£k 8b- £k 



(3.21) 



The states | a) and | b) are, explicitly, 



a) = b\p l ,a)d 1; (p2,a)\0), 

b) = b\ P3 ,P)d\p 4 ,P)\0), 



(3.22) 



where p\, p2 denote momenta and a, j5 denote color which is summed over. Explicitly, p\ = 
(pt'Pl) etc -' wnere P\ is the longitudinal component and p\ is the transverse component. 
For simplicity of notation, we will denote the transverse component of momenta without the 
superscript 1. 

The free part of the Hamiltonian leads to the matrix element 



(a\H\b) 



m 2 + p\ ^ m 2 + p\ 



Pi 



+ 



pi 



2(2n) 2 p+8 2 ( Pi - p 3 )2(2n) 2 p+8 2 (p 2 - P 4)8 a p. (3.23) 
From the four fermion interaction, we get the contribution 

-\g 2 {T a T a ) aa ^\^ 2 2(2Ti) 2 8 ap . (3.24) 

Next we evaluate the contribution from the second order term. The intermediate state | k) is 
any state in higher Fock space, e.g., qqg, qqgg, qqqq, and so on. In the lowest nontrivial order 
we take | k) a quark, anti-quark, gluon (qqg) state. This intermediate state gives rise to both self 
energy and gluon exchange contributions. 

The self energy contributions are 



g l C f 8 ap p12(2k) 2 8 2 ( Pi -p 3 ) p^2{2kY8\ P 2-pa) 
dk^dki 



2(2jr) 2 (tf-*? 



„ (pi — ki) k\ p\ . m . m I 1 



(Pt-K) ' *j pi 

„ (p\ —k\) ki p\ . m . m 
I (pi-ki) + kl + p+ +l k+ l p\ 



k+ p\ J <§\ 



+8 Z C f 8 aj5 pX2{2%) 2 8 2 { Pl -p 3 ) pX2{2kY8\ P 2-Pa) 
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with 



/: 



dkjdkj \ (p2 — £2) ^2 P2 . m . m 



2(2jf)2(p+-*J)1 *J p+ k+ pt 

~{pi — ki) k2 P2 .m .ml ,„ „^ 

-2 ? +^ + f T + l TT- 1 — ' (3 - 25) 

(J? 2 ^2 J K 2 ^2 ^2 ^2 J 



p 2 +m 2 m 2 + ^2 (pj-Jfej^ 



^ 2 = pj + m 2 m 2 + fc 2 ( P2 -fc 2 ) 2 

^2 ^2~ (P2 ~ ^2) 



and the color factor Cf = (T a T a ) = for N number of colors. The gluon exchange contri- 
butions are 



g Cfliln) 8 (pi+P2-P3-P4)ypiP2PiPj 

2 (P1-P 3) , P3 , Pi , m , . m ] f ^ (pi -p 3 ) _ p 2 , P4 , , m . m 
(Pl~P 

ie(p+- P +) 1 1 i 



+ 



2 (Pl-Pt) I m2 +P4 (P1-P3) 2 ^+P2 m2 +P\ (P1-P3) 2 ^+P3 

P4 (P1-P3) pi pi (pi-pi) pi 



g Cf 2(2%) 8 (pi+P2-p3-P4)yPiP2PiP4 

(j>3-Pi) , P3 , Pi .™ .m \ \ (ps-pi) p 2 t Pa .m . m 



(P3-P1) P3 Pi P3 Pi J I (Pi -Pi) Pi Pa P2 Pa 

i °ipj-Pi) { ! + 1 I (3 27) 

2 (pt-pt) I m2+P 2 (P3-Pl) 2 m2 +P4 m2 +P3 (P3-Pl) 2 ™ 2 +p\ f 

I pi (pi-pi) pi pi (pi-pi) pi I 

After the construction of H e ff in the two particle space, we proceed as follows. Consider 
the bound state equation 

M 2 + P 2 

H ef f I V> = I *> (3.28) 

where P + , P, and M are the longitudinal momentum, the transverse momentum and the invariant 
mass of the state respectively. The two particle (qq) bound state | is given by 
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/3 J y/2{2n)2p+ J y/2{2K)*pt 

^2(2k) 2 P+8 2 (P -p3- pa) (3.29) 

with the normalization 

(V(Q)\y(P)) = 2(2n) 2 P + 8 2 (P-Q) (3.30) 

provided 

J J dp\dp\ | (h(P;pi,P-pi) | 2 = l. (3.31) 

We symbolically represent the above state as 

l*>=X>yl7>- (3-32) 

j 

Taking projection with the state (/ |= (0 | d(p2,0c)b(p\,a), we get the effective bound state 
equation, 

M 2 + P 2 



P+ 



-hi = H 0i<hi + Y,(i I Hieff I j) <hj- (3-33) 

Introduce the internal momentum variables (x,k) and (y,q) via p\ = xP + , p\ = xP + k, p^ = 
(1 -x)P+, p 2 = (1 -x)P-k, p+ = yP+, P3 =yP + q, p+ = (1 -y)P+, p 4 = (1 -y)P-q and 
the amplitude fa(P; pi, Pi) = ^=W 2 (x,k). 

The fermion momentum fractions x and y range from to 1. To handle end point sin- 
gularities, we introduce the cutoff £ < x,y < 1. This does not prevent the gluon longitudinal 
momentum fraction (x — y) from becoming zero and we introduce the regulator 8 such that 
I x — y | > 8. To regulate ultraviolet divergences, we introduce the cutoff A on the relative trans- 
verse momenta k and q. We remind the reader that in the superrenormalizable field theory under 
study, only ultraviolet divergence is in the fermion self energy contribution which we remove 
by a counterterm before discretization. 

The bound state equation is 

x(l-x)J— - s ^ x ^- 4 ^Cfjdydq ¥ 2(y,q)j^ 2 
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■ , 2 m 2 + k 2 



2(2k) : 



C f / dydq \\fi{y,q) 



i y 

Ys' 



(3.34) 



The self energy contribution 



S = 



2{2Kf Cf 



' xy 



q , k 2(k-q) \ 2 , m 2 (x-y) 2 
y + x -J^yj) + *V 



j^dy J dq 

+ j; dy j dq(1 - x)(l . y) ^ItT- 

The boson exchange contribution 



(ky — qx) 2 + m 2 (x — y) 2 

q I _k_ I 2(g-fc) \ 2 m 2 (y-x) 2 

1-y ^ 1-x ^ (y-x) J + {\-x)i{\- y y 



[k(l -y) -q(l -x)] 2 + m 2 (x-y) 2 



y 



0(x-y) 

(x-y) 



1 



X 



m 2 +q 2 , (k-q) 2 

L y + {x-y) 
K(k,x,q,y) + iVi 



m 2 +k 2 m--ri\- j 

1— x ' x—y 



+ ■ 



1 



2 +fc 2 (fc_ g )2 m 2 +g 2 



B{y-x) 
(y-x) 



1 



+ 



1 



2 +k 2 (q-k) 2 q2+m 2 m 2 +q 2 (q-k) 2 m 2 +k 2 



(y-x) 



l—y y—x 



l-x 



X 



K(q,y,k,x) + iVj 



where 

K(k,x,q,y) = 



q ^k 2{k — q) 
y x (x-y) 



)(r 



2(k-q)\ m 2 (x-yf 



'\k-q) \ 
(x-y) J 



-y 1 —x (x — y) J xy(l— x)(l— y) 



(3.35) 



(3.36) 



(3.37) 



m 



xy(l-x)(l-y) 



[4(2-y-3x) + &(3y + x-2)]. 



3.4 Divergence Structure 



(3.38) 



Now, let us analyze the divergence structure of the effective bound state equation. We 
encounter both infrared and ultraviolet divergences. 

3.4.1 Ultraviolet Divergences 

First consider ultraviolet divergences. In the super renormalizable field theory under con- 
sideration, with the terms appearing in the canonical Hamiltonian as normal ordered, ultraviolet 
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divergence is encountered only in the self energy contributions. To isolate the ultraviolet diver- 
gence, we rewrite the self energy integrals as 

- 2 rx r+A r (x + y) 2 Am 2 



S = 



'2{2n) lCf 
c 

2{2n) 2 f 



x r+A 

dy / dq 

J-A 
1 r+A 

dy / dq 

J-A 
2 



xy(x — y) 2 (ky — qx) 2 + m 2 (x — y) 2 



(2-x-y) 



Am 1 



(3.39) 



{y-x)\\ -x){l -y) [k(l -y) -q{\ -x)] 2 + m 2 (x-y) 2 

The first term inside the square brackets in the above equation is ultraviolet divergent, which 
we cancel by adding an ultraviolet counterterm given by 



CT = + 



+A 



(x + y) 2 /•! 
o ay xy(x-y) 2+ Jx ay \y-x) 2 (\-x)(\-y) 



(2-x-yf 



(3.40) 



After the addition of this counterterm, the bound state equation is ultraviolet finite. 
3.4.2 Infrared Divergences 

The infrared divergences (IR) that appear in the bound state equation are of two types: 
(1) light front infrared divergences that arise from the gluon longitudinal momentum fraction 
x g = 0, (2) true infrared divergences that arise from gluon transverse momentum k g = and 
gluon longitudinal momentum fraction x g = 0. The IR divergences of type (1) are generated 
due to elimination of the constrained degrees of freedom. 

Cancellation of Light-front Infrared Divergences in the Effective Bound State Equation 

First consider light front infrared divergences. The effective bound state equation, Eq. 
(3.34), explicitly has a linear light front infrared divergent term , _ . 2 coming from instan- 

(x y) 

taneous gluon exchange. The most divergent part of the numerator of the transverse gluon 
exchange term in this equation is —4 • After combining the terms, the linear infrared di- 
vergent term is completely canceled and the resultant effective bound state equation takes the 
form 



M 



m 2 + k 2 



x(l —x) 



\l/ 2 (x,k) = Si \i/ 2 (x,k) - 



8 



2{2k) 2 



C f / dydq \\fi{y,q) 



1 

x - 
2 



Vi V 2 , 7 / 1 1 

— + — + iV/ I 1 

Ei E 2 \Ei E 2 



(3.41) 



31 



The self energy contribution, made ultraviolet finite by the addition of the counterterm is 



Si 



rx r+A 

Cf dy dq 

JO J -A 



2{2n) 2 

.2 



4m 2 



(ky — qx) 2 + m 2 (x — y) 2 

I r +A 4 m 2 



+ 2(L)2 C/ / dy L A dq [k(\-y)-q(\-x)] 2 + m 2 (x-y) 2 - ° A2) 



The energy denominator factors are 



1 



xy 



1 



(l-x)(l-y) 



E\ [ky — qx] 2 + m 2 (x — y) 2 ' E 2 [k(l — y) — q{\ — x)] 2 + m 2 (x — y) 2 ' 
The vertex terms are 

Vi = 6{x-y)U{k,x,q,y) + 9{y -x)U(q,y,k,x), 



(3.43) 



(3.44) 



V 2 = 0(x-y)U(k,l-x,q,l-y) + e(y-x)U(q,l-y,k,l-x), 



with 



U{k,x,q,y) 



. m 

4 



m 2 (x — y) z 



xy xy(l—x)(l—y) 



+ 



+ 



-2 



1 



+ 2- 



1 



(3.45) 



y(l-y) x(l-x) (x-y)x(l-x) (x-y)y(l-y) 



+ ■ 



kq 



+ 



kq 



+ 2- 



kq 



l-2y 



1 -2x 



(3.46) 



x(l-y) y(l-jc) (Jc-y) Ly(l-y) x(1-jc). 

In addition to the \ singularity which is canceled, transverse gluon exchange contributions also 
contain — singularity which is removed by the principal value prescription. Cancellation of this 

Xg 

singularity is an appealing feature of the Bloch effective Hamiltonian in contrast to the Tamm- 
Dancoff effective Hamiltonian where the singularity cancellation does not occur because of the 
presence of invariant mass in the energy denominator [21]. 

"True" infrared divergences 

Next we consider true infrared divergences. Consider the self energy integrals. The energy 
denominators in these expressions vanish when k = q and x = y which correspond to vanishing 
gluon momentum. By carrying out the integrals explicitly, in the limit A^<»we get, 

- 2 M x 1 I-* 



Si 



2k 



C 



f 



, x 
In - + 

x o l—x 



In 



(3.47) 
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Thus the singular part of self energy is 

Sl = - If °f xjhx) ln 5 - (3 - 48) 

The infrared divergent contribution from self energy gives a positive contribution to the fermion 
mass. It is important to note that the vanishing of energy denominator is possible also in (3+1) 
dimensions, but in that case we do not encounter any divergence. It is the peculiarity of (2+1) 
dimensions that the vanishing energy denominators cause a severe infrared divergence problem. 

The same vanishing energy denominators occur also in the one gluon exchange contribu- 
tions. Let us now consider various terms in the numerator separately. The terms proportional to 
4m 2 arose from the denominator of the transverse gluon exchange. A straightforward calcula- 
tion shows that this term leads to both finite and infrared divergent contributions. The infrared 
divergent contribution is given by 

^T- C f TT x ln 5 (3-49) 
2k 1 x(l-x) 

which exactly cancels the infrared divergent contribution from self energy. The finite part, in 
the nonrelativistic limit, can be shown to give rise to the logarithmically confining potential. 
Next we have to consider the remaining terms in the numerator. Rest of the terms proportional 
to m 2 are multiplied by (x — y) 2 so that they do not lead to an infrared divergence problem. The 
numerator of the imaginary part vanishes at k = q, and x = y and hence is also infrared finite. 
It is easy to verify that the rest of the (transverse momentum dependent) terms in the numerator 
does not vanish when the denominator vanishes and hence the resulting bound state equation is 
inflicted with infrared divergences arising from the vanishing energy denominator. This prob- 
lem was first noted in the context of QED in (2+1) dimensions by Tarn, Hamer, and Yung [22] 
but was not investigated by these authors. We remind the reader that this is a peculiarity of 
(2+1) dimensions which provides us a unique opportunity to explore the consequences of the 
vanishing energy denominator problem. 

3.5 Numerical study of the bound state equation 

Once we derive the effective Hamiltonian, we need to diagonalize it nonperturbatively. For 
that, we convert the integral equation into a matrix equation with the use of Gaussian Quadra- 
ture. (For details of the numerical procedure see Appendix C.) The color factor Cf is set to 1 for 
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all the numerical calculations presented. As mentioned before, an important feature of gauge 
theories on the light-front is the presence of linear infrared divergences. They appear in the 
canonical Hamiltonian in instantaneous four fermion interaction term. When the qqg states are 
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Figure 3.2: Cancellation of infrared divergence. Full line denotes the full Hamiltonian. (a) 
shows the cancellation of light-front infrared divergence by switching on and off the instanta- 
neous interaction. Filled circles - without instantaneous interaction, (b) shows the cancellation 
of logarithmic infrared divergence by switching on and off the self energy term. Filled circles - 
without self energy. The parameters are g = 0.2, e = 0.00001, m = 1, K = 20, ri\ = 40, ri2 = 50. 



integrated out completely in perturbation theory, they also appear in the effective four fermion 
interaction and cancel against each other. Non-cancellation of this divergence is a major feature 
of similarity renormalization approach. We first address the issue of how linear divergences 
manifest in the non-uniform grid of the Gaussian Quadrature and how well it can handle lin- 
ear light front infrared divergence. We have studied numerically discretized versions of Eq. 
(3.34) where the divergences are present separately in the discretized version together with the 
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counterterm given in Eq. (3.40). For g = 0.2, we have calculated the eigenvalues with and with- 
out the instantaneous interaction. The result presented in Fig. 3.2(a) for the lowest eigenvalue 
shows that the Gaussian Quadrature can handle the cancellation very efficiently. 

After the cancellation of linear light-front infrared divergence, a logarithmic infrared diver- 
gence which arises from the vanishing energy denominator survives in the bound state equation. 
Here we have to distinguish two types of terms. First type, where the coefficient of the loga- 
rithmic infrared divergence is independent of the fermion transverse momentum and the second 
type where the coefficient is dependent. Self energy and Coulomb interaction are of the first 
type. In the weak coupling limit, since the wavefunction is dominated by very low transverse 
momentum, we anticipate that contributions of the second type will be dynamically suppressed 
even though both are multiplied by the same coupling constant. This is especially true of any 
discrete grid which automatically imposes a lower limit on the smallest longitudinal momen- 
tum fraction allowed. Thus at weak coupling, even if there are uncanceled infrared divergences 
(divergences of the second type), they may not be significant numerically whereas divergences 
of the first type are significant. By switching the self energy contribution off and on, we have 
studied this interplay. The lowest eigenvalue with and without self energy contribution is plot- 
ted in Fig. 3.2(b). This shows the cancellation of the dominant logarithmic infrared divergence. 
Since there are still uncanceled infrared divergences in the bound state equation (with coeffi- 
cient proportional to fermion transverse momenta) these figures further illustrate the fact that 
such divergences are not numerically significant at weak coupling. 

As the strength of the interaction grows, wavefunction develops medium to large transverse 
momentum components and the infrared catastrophe triggered by the vanishing energy denom- 
inator becomes manifest numerically. This is illustrated in Table 3.1 where we present the 
variation with <5 of the first five eigenvalues for two different choices of the coupling g. The 
table clearly shows that on a discrete grid, the uncanceled infrared divergences due to the van- 
ishing energy denominator problem are not numerically significant at weak coupling but their 
effect is readily felt at a stronger coupling. 
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8 


5 


eigenvalues (M 2 ) 


0.2 


0.01 


4.0870 


4.0972 


4.0972 


4.0973 


4.0973 


0.005 


4.0901 


4.1066 


4.1099 


4.1100 


4.1112 


0.001 


4.0913 


4.1113 


4.1122 


4.1181 


4.1209 


0.0001 


4.0913 


4.1113 


4.1122 


4.1181 


4.1209 


0.00001 


4.0913 


4.1113 


4.1122 


4.1181 


4.1209 


0.6 


0.01 


4.5735 


4.7337 


4.7667 


4.7832 


4.8277 


0.005 


1.9094 


1.9415 


3.1393 


3.1399 


4.5697 


0.001 


-187230.4 


-187225.4 


-186664.9 


-186664.8 


-31506.9 


0.0001 


-187230.4 


-187225.4 


-186664.9 


-186664.8 


-31506.9 



Table 3.1: Variation with 5 of the full Hamiltonian. The parameters are «i=40, «2=50, 
£=0.00001, ?c=20.0 ink= ^tan(f) 



3.6 Reduced Model 

In this section we consider a model Hamiltonian free from infrared divergences constructed 
by dropping the transverse momentum dependent terms from the numerator of the effective 
Hamiltonian. For convenience, we further drop the terms proportional to (x — y) 2 and the imag- 
inary part. This defines our reduced model which is also ultraviolet finite. The equation gov- 
erning the model is given by 

\l/2(x,k) = Si \l/2(x,k) (3.50) 

The self energy contribution Si is the same as given in Eq. (3.42). The boson exchange contri- 
bution SS is given by 

g 2 f l f +A Am 2 
^ = ~ W C/ /o dy L dq (ky- qX ) 2 + l n 2 ( X -y) 2 
g 2 f 1 f +A Am 2 

Note that in the above approximations we dropped only the term sick with vanishing energy 
denominator and not so important imaginary terms and Eq. (3.50) still represents a relativistic 
bound state equation. Though the rotational symmetry is not manifest in this equation, Eq. 
(3.50) in the nonrelativistic limit reduces to a Schrodinger equation with explicit rotational 
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2 m z + k 2 



x(l—x) 



symmetry (see Appendix D). This model provides us an opportunity to study the simplest 
manifestation and possible violation of rotational symmetry in the context of light-front field 
theory. 

3.6.1 Numerical study of the reduced model 

Again we discretize the Eq. (3.50) by Gaussian Quadrature. The convergence of the eigen- 
values as a function of the number of grid points is presented in Table 3.2. In this table we also 
present the (in)dependence of eigenvalues on the momentum cutoff. 



n\ 


«2 


eigenvalues (lowest five) (k =10.0) 


20 


20 


4.08926 


4.10605 


4.10768 


4.11061 


4.11085 


30 


30 


4.09045 


4.10909 


4.11038 


4.11516 


4.11699 


40 


30 


4.09045 


4.10913 


4.11035 


4.11524 


4.11697 


40 


40 


4.09102 


4.11052 


4.11154 


4.11711 


4.11951 


40 


50 


4.09136 


4.11133 


4.11222 


4.11811 


4.12096 


50 


50 


4.09136 


4.11135 


4.11219 


4.11816 


4.12095 


50 


60 


4.09158 


4.11188 


4.11263 


4.12189 


4.12290 


46 


60 


4.09158 


4.11187 


4.11264 


4.11877 


4.12189 


46 


66 


4.09168 


4.11212 


4.11284 


4.11905 


4.12231 


46 


74 


4.09179 


4.11237 


4.11305 


4.11934 


4.12276 


n\ 




eigenvalues (lowest five) (k=20.0) 


46 


74 


4.09179 


4.11240 


4.11301 


4.11940 


4.12273 



Table 3.2: Convergence of eigenvalue with n\ and ri2 (reduced model). The parameters are 
m =1.0, g=0.2, e = 0.00001. 



(2+1) dimensions provide an opportunity to study the manifestation and violation of rota- 
tional symmetry in light front field theory in a simpler setting compared to (3+1) dimensions. 
The absence of spin further facilitates this study. Rotational symmetry in this case simply im- 
plies degeneracy with respect to the sign of the azimuthal quantum number / (see Appendix D). 
Thus we expect all / ^ states to be two fold degenerate. By a suitable change of variables, 
one can easily show that our reduced model, in the nonrelativistic limit reduces to Schrodinger 
equation in two space dimensions with a logarithmic confining potential. In the weak coupling 
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8 


eigenvalues 


0.2 


This 
work 


4.0918 
(4.1227,4.1235) 


(4.1124,4.1130) 
(4.1268,4.1273) 


4.1194 
(4.1298,4.1303) 


Koures 
(Ref. [23]) 


4.0925 (1 = 0) 
4.1260 (1 = 2) 


4.1144 (1 = 1) 
4.1303 (1 = 1) 


4.1214 (1 = 0) 
4.1340 (1 = 3) 


0.6 


This 
work 


4.5856 
(4.8767,4.8816) 


(4.7741,4.7821) 
(4.9094, 4.9184) 


4.8390 
(4.9458,4.9481) 


Koures 
(Ref. [23]) 


4.5806 (/ = 0) 
4.8827 (1 = 2) 


4.7777 (/ = 1) 
4.9205 (1 = 1) 


4.8409 (1 = 0) 
4.9545 (/ = 3) 



Table 3.3: Reduced model. The parameters are Hi =46, ^2=74, £ = 0.00001, m=1.0. k = 
tan(qn/2) /jc, K = 20.0. Eigenvalues within () are ±1 degenerate (broken) states. 





n\ n 2 


eigenvalues 




40 50 


18.217 


f 30.702 > 
V 33.499 J 


35.206 


( 39.955 > 
^41.159 ) 




(A\332\ 
^43.271 ) 




/44.134\ 
v 45.272 ) 


46 70 


18.276 


f 30.774 ^ 
V 33.616 j 


35.318 


( 40.106 ^ 
V41.331 J 




f 41-483^ 
V 43.477 ) 


( 


' 44.375 \ 
, 45.503) ) 




40 50 


18.980 


( 31.507 \ 
\ 34.219 ) 


35.826 


f 40.406 \ 
^41.888 J 


^ 41.921 \ 
V 43.788 J 


( 44.345 \ 
V 45. 163 J 


46 70 


19.008 


f 31-542^ 
^34.319 J 


35.935 


( 40.626 \ 
^42.031 J 


( 42.088 \ 
^44.010 J 


( 44.647 \ 
V 45.780j 



Table 3.4: First few eigenvalues in the reduced model. The parameters are g=5.0, m=1.0, 
£=0.00001. (I) for the parametrization k = wAra/((l - u 2 )A + m), A = 40.0. (II) for the 
parametrization k = tan(un /2)/?c, K = 10.0. Eigenvalues within ( ) are ±Z degenerate (bro- 
ken) states. 



limit, since C/ is set to 1, we can compare our results of the reduced model (where we do not 
make any nonrelativistic approximation) with the spectra obtained in nonrelativistic QED2+1. 
Tarn et al. [22] solved the radial Schrodinger equation in momentum space for / = states 
and Koures [23] solved the coordinate space radial Schrodinger equation for general /. Since 
we are solving the light front bound state equation, rotational symmetry is not at all manifest. 
However, at weak coupling we expect that the spectra exhibit rotational symmetry to a very 
good approximation. Our numerical results are compared with those of Koures in Table 3.3 for 
two values of the coupling. At g = 0.2 we find reasonable agreement with the degeneracy in 
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the spectrum. Even at g = 0.6 the violation of rotational symmetry is very small. Splitting of 
levels which are supposed to be degenerate become more visible at very strong coupling as can 
be seen from Table 3.4 for g = 5.0. 

Along with the eigenvalues, the diagonalization process also yields wavefunctions. We have 
plotted the wavefunctions corresponding to the first four eigenvalues in Fig. 3.3 as a function 
of x and k. All wavefunctions are normalized to be Jq dx J dk y 2 (x,k) = 1. The lowest state 




Figure 3.3: The wavefunctions corresponding to the lowest four eigenvalues of the reduced 
model as a function of x and k. The parameters are g = .2, £ = 0.00001, m = 1, K = 10, n\ = 46, 
n.2 = 74. (a) Lowest state, (b) first excited state, (c) second excited state, (d) third excited state. 
The first and second excited states should be degenerate in the absence of violation of rotational 
symmetry. 



is nodeless and corresponds to / = 0. The next two states correspond to I = 1 and have one 
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node. It is interesting to note the way the node appears in the wavefunctions which correspond 
to degenerate levels. Since the rotational symmetry cannot be manifest in the variables x and 
k, how can the wavefunctions still indicate this? From Fig. 3.3, it is clear that the way this 
problem is resolved is by one wavefunction having a node in k and the other wavefunction 
having a node in x. Thus even if we did not know about the underlying symmetry from other 
means, the light-front wavefunctions have a subtle way of indicating the symmetry. 

3.7 Summary 

The numerical solutions of the bound state equations are performed by using the Gaussian 
Quadrature (GQ) which is a straightforward procedure to solve the integral equation by con- 
verting it into a matrix equation. The investigation shows the efficiency of the GQ method 
in handling linear and logarithmic light-front infrared divergences. The manifestation of rota- 
tional invariance in light-front framework is demonstrated very clearly in the reduced model. 
But, our study of the Bloch effective Hamiltonian indicates that in the context of Fock space 
based effective Hamiltonian methods to tackle gauge theories in (2+1) dimensions, approaches 
like similarity renormalization method are mandatory due to uncanceled infrared divergences 
caused by the vanishing energy denominator problem. It is important to recall that Bloch effec- 
tive Hamiltonian is generated by completely integrating out the intermediate gluons irrespec- 
tive of whether they are low energy or high energy which has no clear justification specially 
in a confining theory. Once we have obtained quantitative measures of the vanishing energy 
denominator problem and the nature of the spectra at weak coupling of the Bloch effective 
Hamiltonian, the next step is to study QCD2+1 in the similarity renormalization approach which 
avoids the vanishing energy denominator problem. An important issue here is the nature of new 
effective interactions generated by the similarity approach. It has been shown that in (3+1) di- 
mensions, similarity approach generates logarithmic confining interactions [16] which however 
breaks rotational symmetry. It is interesting to investigate the corresponding situation in (2+1) 
dimensions. 
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CHAPTER 4 



Similarity Renormalization Group Approach to Meson Sector 

in (2 + 1 ) Dimensions 

4.1 A Brief Review of Renormalization Group Approach 

In the previous chapter we have studied the meson bound state problem in (2+1) dimensional 
QCD with Bloch effective Hamiltonian where we encountered the problem of vanishing energy 
denominators. We also concluded that Similarity Renormalization Group (SRG) approach is 
mandatory to get rid of the problem. In this chapter we discuss the same bound state problem 
in SRG scheme. But before that, we recapitulate the basic concept of Renormalization Group 
and explain why SRG is preferred over original Wilsonian Renormalization Group approach. 

It is well established that the most important tool to construct a low energy effective field 
theory is the Renormalization Group (RG). The concept of Renormalization Group was first 
introduced by Stueckelberg and Peterman [1] and Gell-Mann and Low [2] and further devel- 
oped by Bogoliubov and Shirkov [3]. RG as a practical tool to construct effective theory was 
developed by K.G. Wilson [4, 5, 6]. The aim of RG is to simplify the problem with many 
energy (length) scales involving many degrees of freedom which are coupled through the inter- 
actions. In the Wilsonian approach the cutoff on energy is lowered and the number of degrees 
of freedom is reduced in an iterative way and in each step one has to construct the effective 
interactions for the effective degrees of freedom. The simplification of Renormalization Group 
lies in the hope that the effective interactions are local interactions i.e., only nearby degrees 
of freedom are directly coupled by the interactions which holds true for original local theories 
we normally deal with [5]. Construction of nondiagrammatic RG transformations enables one 
to solve them numerically in a computer and hence the problems which cannot be solved by 
Feynman diagrams can be solved by using RG. 
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The starting point of the RG transformations is a bare Hamiltonian Hq with cutoff Ao- The 
transformation t converts Hq to Hi, Hi to H2 etc. as the cutoff is lowered Ao > Ai > A2 in 
each step and thins the degrees of freedom. This transformation is to be iterated until one gets 
the effective Hamiltonian at the desired low energy scale A. RG transformation is the evolution 
operator of the Hamiltonian as the cutoff changes. Here I should mention that the RG trans- 
formation is free of (UV or IR) divergences, since in each step, a momentum integral involves 
only a finite range of momentum. But the divergences occur as a result of many iterations of 
the RG transformation. The logic of RG transformation is best explained by the "triangle of 
Renormalization" [6]. Suppose one is interested to solve a theory at the energy scale 1 in some 
suitable unit. Consider, for example, that in each step we lower the cutoff by a factor of 1/2, 
i.e., we have a discrete set of cutoffs, A = 2 N for N = 1 , 2, 3,4, • • • , 00. For each N, we apply the 
RG transformation to produce a sequence of effective Hamiltonians Hq, Hf , Hj , ■ ■ ■ until we 
reach H$ with cutoff A = 1 . Thus the transformation produces a triangle. 

A = 4 

A = 2— Hi 
I T 

A =l- flj H\ 

The N — > 00 limit along any row produces the infinite cutoff limit or the renormalized Hamilto- 
nian. For example, the N — > °° limit of Ao = 1 row generates the Hamiltonian renormalized at 
the scale A = 1 . 

The fixed point of the transformation is defined by 

t(H*)=H*. (4.1) 

The fixed point of a transformation is a property of % and does not depend on the initial Hamil- 
tonian Hq. 

But, this formalism to construct an effective low energy Hamiltonian acting on a limited 
Fock space is again plagued with the problem of vanishing energy denominator. If we view the 
Hamiltonian as a matrix, the energy cutoff limits the size of the matrix and as one lowers the 
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Hi 



ttN 

H N 2 



It It 



ttN 



It It 



H 



N 



N 



cutoff, the matrix size is also reduced as shown in Fig. 4.1(a) [7]. The matrix elements near 
the diagonal region (gray region in Fig.4.1(a)) involve states with almost same energy (nearly 
degenerate) and involve degenerate perturbation theory for calculation of the effective interac- 
tions. In nonperturbative theory like QCD we do not even know how to do that. To overcome 










(a) (b) 
Figure 4.1: Effective Hamiltonian as the cutoff is lowered (solid to dashed line), (a) energy 
cutoff, (b) cutoff on energy difference. 



this problem Glazek and Wilson [8] and Wegner [9] independently developed the similarity 
renormalization group (SRG) technique to calculate low energy effective Hamiltonian. In this 
approach the cutoff is not on the energy of the states but on the energy difference of the states. 
In place of removing the states, off-diagonal matrix elements involving large energy transfer 
are removed. If the free energy difference between two states is larger than the cutoff, interac- 
tions between those two states are then removed from the effective Hamiltonian. The working 
prescription for SRG is as the following. Again, as in the standard Wilsonian RG, the starting 
point is a finite and bare cutoff Hamiltonian Hg at some ultra-violet cutoff A. Define a similar- 
ity transformation that converts Hg into a band-diagonal H 0l as the energy scale is lowered to 
<7i and removes the coupling between the states with energy difference greater than 0"i . This 
process should be repeated until one produces an effective Hamiltonian H a at a low energy 
scale a. It can again be viewed as a new "triangle of renormalization" [10] as discussed above 
in the context of Wilsonian RG transformations. The way it works can be compared with the 
standard numerical algorithm for matrix diagonalization [11] where to keep control over the 
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complexities, the matrix is first brought to a tri-diagonal form which is then diagonalized. In 
the SRG approach the Hamiltonian is brought perturbatively into a band-diagonal form (see 
Fig. 4.1(b)) which is then diagonalized nonperturbatively. Since it does not involve any energy 
jump below the energy scale o, the energy denominator cannot be smaller than o and is free 
from the problem of vanishing energy denominator. 

The diagonalization of the Hamiltonian using SRG is a two step process. In the first step, 
by removing any direct interactions between states with energy difference larger than the cutoff 
one arrives at an effective Hamiltonian which is in a band diagonal form. At this step, one can 
identify the ultraviolet divergent part of the counterterms needed to be added to the Hamiltonian 
to remove ultraviolet divergences. It is quite advantageous to treat the ultraviolet divergences 
perturbatively especially in gauge theories since one can avoid pitfalls of other effective Hamil- 
tonian approaches. For example, it is well known that a simple truncation of the Fock space 
(like a Tamm-Dancoff truncation) leads to uncanceled divergences as a result of violations of 
gauge symmetry. 

In the second step, the effective Hamiltonian is diagonalized exactly. It is important to note 
that in the process of removing the interactions with very large energy exchange we integrate out 
small x gluons, i.e., gluons having small longitudinal momentum fraction. Since the vacuum 
is trivial, it is hoped that, as a result of integrating out small x gluons which are sensitive to 
long distance physics on the light front, the effective Hamiltonian may contain interactions 
responsible for low energy properties of QCD. Indeed Perry [12] found a logarithmic confining 
interaction in the qq sector in the lowest order effective interaction. 

Initial bound state studies in the similarity renormalization approach worked in either the 
non-relativistic limit [13] or in the heavy quark effective theory formalism [14] to investigate 
heavy-quark systems. Only in last few years, some works have been done [15] in the context of 
glueball spectrum to address many practical problems, especially the numerical ones that one 
faces in this approach. Since the conceptual and technical problems one encounters in QCD are 
numerous, we initiated a study of bound state problems in QCD in (2+1) dimensions [16, 17]. 
Our main motivation is not the fitting of data but a critical evaluation of the strengths and 
weaknesses of the various assumptions and approximations made in the similarity approach. 
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In the previous chapter we have discussed the meson sector of (2+1) dimensional light front 
QCD using a Bloch effective Hamiltonian [18] in the first non-trivial order. The resulting two 
dimensional integral equation was converted into a matrix equation and solved numerically. We 
have already discussed in detail the problem of vanishing energy denominator which leads to se- 
vere infrared divergences in (2+1) dimensions and came to the conclusion that in the context of 
Fock space based effective Hamiltonian methods to tackle gauge theories in (2+1) dimensions, 
approaches like similarity renormalization method is mandatory due to uncanceled infrared di- 
vergences caused by the vanishing energy denominator problem. Now, we discuss the similarity 
renormalization approach in the first non-trivial order to the same problem. The detail of simi- 
larity renormalization theory for the Effective Hamiltonian in both Glazek- Wilson and Wegner 
approaches is discussed in Appendix E. 

4.2 Effective bound state equation in the qq sector in SRG scheme 

In similarity renormalization approach due to Glazek and Wilson, to second order, the inter- 
acting part of the effective Hamiltonian at a scale o is given by (see Appendix E for details) 



H Iaij - -Y^ HmikHmk j 

k 



gaijk 8 a jik 



Pk-Pj Pk-^ 



(4.2) 



where 

d_ 

ic 
d 



f°° d 
gaijk = faij J do' f a ' ik -^pf a ij k , 

r i d 

go jik = faij J a do fa'jkj^jfa'ik- ( 4 -3) 
. The similarity factor faij(x) is such that 

when o 2 » AM?-, f(x) = 1 (near diagonal region); 

when o 2 « AMfj, f(x) = (far off diagonal region); 

in between f(x) drops from 1 to (transition region) . (4.4) 

Here AM 2 j = (Mf—M 2 ) denotes the difference of invariant masses of states i and j. We restrict 
ourselves to the qq sector. Then the states involved in the matrix elements i and j refer to qq 
states and k refer to qqg states. 
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Following the steps similar to the ones outlined in Chapter 3, we arrive at the bound state 
equation 



M 



2 m 2 + k 2 



x(l — x). 



S \i/ 2 (x,k) -*2(h)2 C f J dy j dqfai J Tx~Z 



2(2^) 2 



{x-yf 

(4.5) 



Here x and y are the longitudinal momentum fractions and k and q are the relative transverse 
momenta. We introduce the cutoff e such that £ < x,y < 1 — £. We further introduce the 
regulator 8 such that | x — y |> <5. Ultraviolet divergences are regulated by the introduction of 
the cutoff A on the relative transverse momenta k and q. The self energy contribution 



q k 2(k-q) 



^2 C fJ Q d y Jdqe(x-y) [l-faik] xy 

ff [dy j dq9(y-x) [I - f aik ] (l-x)(l-y) 



+ x 2 y 2 



(ky — qx) 2 + m 2 (x — y) 2 



2(27T) 



x 



q | k | 2{g-k) 
l-y -T l- x "I" (y_ x ) 



+ 



m 2 (y-x) 2 
(l-x)^(l-y)^ 



[fc(l-y)-<7(l-*)] 2 + m 2 (*-v) 2 ' 
The boson exchange contribution 



(4.6) 



Ojx-y) 



L y 



Sojik , Saijk 

+q 2 , (fc-g) 2 m 2 +fc 2 m 2 +fc 2 , 
v ~ T " (x— y) x 1— x x— y 



x K{k,x,q,y) + iV/] 
0(y--y) gg/U 



n 2 +q 2 
l-y 



(y-x) 



Saijk 



X 



i 2 +£ 2 , (g-^) 2 _ g 2 +m 2 + m 2 + g 2 , (q-k) 2 _ m 2 + fc2 
x (y— x) y 1— y y— x 1— x 

K{q,y,k,x) + iV/] , 



(4.7) 



where 



K{k,x,q,y) = 



y x (x-y) 



( 2(*-g) 
1 — y 1 — jc (x — y) 



+ 



+ 



m 



xy(l -x)(l -y) 



m 2 (x — y) 2 
xy(l -x)(l -y)' 



[q(2-y-3x)+k(3y + x-2)}. 



(4.8) 



(4.9) 
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For all the / and g factors, 



t k 2 + m 2 t a 2 + m 2 

Mf = — ^- and M) = (4.10) 
x(l-x) J y(l-y) 



(k-q) q+m z k 2 + m 2 /J11S 
ForjOy, ^ ^ + ^ + (4.ll) 

* — .y j l 



and 



fory> ,,M? = « + 4±^ + ^!. (4 . l2 ) 

y—x l — j ^ 

Before proceeding further, we perform the ultraviolet renormalization. The only ultraviolet 
divergence arises from the term involving the factor l inside the square bracket in Eq. (4.6). We 
isolate the ultraviolet divergent term which is given by 



,2 



S divergent ~ 2(2k) 2 ^^ 



n dy r dq ^r 

Jo J -A xy(x-y) z 

-_y? ' 

y)(x-y) 2 



f 1 J f+ A J (2-x-y) 2 
+ dy dq \ 

Jx+8 J-A (1-X)(1- 



(4.13) 



which is canceled by adding a counterterm. 

4.3 Similarity factors 

Up to this point, we have not chosen any particular form of the similarity factor f(x) . Any 
function that satisfies the criterion (4.4) is a suitable candidate for the similarity factor. Here 
we consider three possible choices of similarity factor for comparative studies in the context of 
meson bound state problem. 

4.3.1 Parameterization I 

In the SRG study of glueball [15], the following form for the similarity factor has been 
chosen: 



,2)2 

faij = e (4.14) 



(AMf Jt 
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9 9 9 9 9 ( l^-hWI^ 

with AM,-,- = Mf — Mf where Mf denotes the invariant mass of the state i, i.e.,Mf = E/^-^- — '- 
Then 

d 



goijk = fcij J do' f a i ik -^pfa'jk 



(AM- 



2\2 



e ^ 



(AM 2 ,) 2 



(AM|)2 + (AM^ 



1 -e 



(4.15) 



For the self energy contribution, i = j and we get 



1 



Saijk — Sajik — Saiik — ^ 



l-e 



2(AM|)2 



(4.16) 



Due to the sharp fall of / with o, the effective Hamiltonian has a strong dependence on a. Note 
that this parameterization emerges naturally in the Wegner formalism (see Appendix E). 

4.3.2 Parameterization II 



St. Glazek has proposed the following form [19] for fcij- 

1 



fcij — 



1 + 



with 



Uqjjjl-UQ) 
«0(l-«<7i/) 



AM 2 . 



2"s 



(4.17) 



(4.18) 



uq a small parameter, and n g an integer. The mass sum ZMfj = Mf + Mj. The derivative 



d ha = 1 n g _ 2(7 

do 



2"s 



(1 - wo) 



2"s 



EMj + a 2 \u J (l-u oij ) 2 " g + l 



1 + 



U ai j(l-UQ) 

uo{l-Uoij) 



2"a 



2 ■ 



(4.19) 



Note that for small u, both 1 — f(u) and 4£ vanish like u 2 " s 



4.3.3 Parameterization III 



For analytical calculations it is convenient to choose [12] a step function cutoff for the 
similarity factor: 



faij = 0(o 2 — AM 2 -) . 



(4.20) 
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Then 



goijk = 9(o — AMfj) e(AM%-a z ) e(AM%-AMf k ). (4.21) 

It is the factor 6(AM 2 k — o 2 ) in g G ijk that prevents the energy denominator from becoming 
small. 

4.4 Analytical calculations with the step function similarity factor 

In this section we perform analytical calculations to understand the nature of the effec- 
tive interactions generated by the similarity factor. The differences in the divergences at some 
places with those discussed in Chapter 3 come due to the similarity factors associated with 
different terms in SRG scheme. As we have already emphasized that SRG is a modification 
over Bloch perturbation theory, we will see here how the divergences are made softer and their 
consequences in similarity group transformed effective bound state equation. Since there are no 
divergences associated with e and A, we suppress their presence in the limits of integration in 
the following equations. 

4.4.1 Self energy contributions 

Consider the self energy contributions to the bound state equation Eq. (4.6). Rewriting the 
energy denominators to expose the most singular terms, we have, 



k 2{k-q)\ 2 . m 2 (x-y) 2 



(x-y) 



g 2 r (' f. 0(x-S-y) 2 

f h dy J dq x-y 11 (k- q ? | M 



(x-y) + y 



.2 



Y q k 2(q-k) \ 2 m\y-x) 1 1 

.V-y + i-* + (y-x) ) + (i-x)Hi-y) 2 \ 



(q-k) 2 _|_ q 2 +m 2 _ k 2 +m 2 
y—x ' l—y 1— x 



(4.22) 



The terms associated with 1 in the curly brackets are the same as in Bloch effective Hamiltonian 
and lead to ultraviolet linear divergent terms which we cancel by counterterms. They also lead 
to an infrared divergent term [16] which remains uncanceled. Explicitly this contribution is 
given by 

g 2 m 2 rx-S r i 

4 2{2nj 2 Cf Jo dy J dq [ky-qx] 2 + m 2 (x-y) 2 
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g 2 m 2 f l f I 
+4 2(2^j2 Cf J x+S dy .1 dq [k{\-y)-q{\-x)] 2 + m 2 {x-y) 2 ' (4>23) 

This is simply indicative of the fact that terms associated with 1 in the curly bracket still has a 
vanishing energy denominator problem. We will address the resolution of this problem shortly. 

Let us next consider new infrared divergences that arise as a result of the modifications due 
to similarity factor. 



Leading singular terms 



Keeping only the most infrared singular terms in the numerators (i.e., for x > y, 4 ^_ q \ 2 and 



for y > x, 4 j^jr) and denominators (i.e., for x> y, ^ x -y) an( * ^ or y >x > )' we nave ' 
* = ^^CfJ^dyJdqeix-d-y)^-^ 



{x-yf 



2{2%f 



8 



2{2k) 2 

The integral with 9 -function similarity factor is given by 

J Jo {x-yr V x-y J Jx+d 



CfJ^dy J dqe(y-x-8)ti ik 4j-±-y. (4.24) 



1 



x+ s {y-xf 



9 a z - 



y — x 



(4.25) 



We change the transverse momentum variable, p = k — q. For x — 8 >y, we set x — y = z and 
for y > x + 8 we set y — x = z. Then, we have, 



I 6 f-i 

2{2k) 2 f 



I6g 2 



2(2n) 



1 



1 



Vx y/l-x. 



(4.26) 



Sub-leading singular terms 



Next we study sub-leading singular terms containing ^ in self energy generated by the 
similarity transformation. They are given by 



Si 



-4 



2(2k) 



2 C f 



Jo J \ x — y J x — y V x y / (k — 



Jx+8 J V y~ x J y~ x ^~x \-y) {q- 



kf 



ci) 2 
(4.27) 
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where we have kept only (k — q) 2 term in the denominator since the rest vanish in the limit 
x — > y. As before, for x — 8 > y, we put x — y = z, k — q = p. With the symmetric integration in 
p, terms linear in p do not contribute. Only potential source of 8 divergence is the p 2 term in the 



singularity. Same situation occurs for y > x. Thus there are no terms divergent in 8 coming 
from sub-leading singular terms. 

4.4.2 Gluon exchange contributions 

Let us next consider the effect of similarity factors on gluon exchange terms. 
Instantaneous gluon exchange 

From instantaneous interaction we have, 



integrand. Since p max = 



G^/z, after p integration - is converted into -4= which is an integrable 




(4.28) 



For the sake of clarity, it is convenient to rewrite this as 




+ 



6(y-x-8) 



^!_(£!_^) +m 2 ( l_l) 

y— x v y x ' v x y ' 



y — x 



(q-k) 2 -(£-^)(y-x)+m 2 (\-\)(y-x) 




(4.29) 



We have to seperately analyze the three types of terms in the numerator. 

First consider terms proportional to m 2 in the numerator. They are given by 




1 1 



x 



[ky — qx] 2 + m 2 (x — y) 2 [k(l —y) —q(l —x)] 2 + m 2 (x — y) 2 



(4.30) 
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which leads to the logarithmic confining interaction in the nonrelativistic limit. Note, however, 
that Eq. (4.30) is affected by a logarithmic infrared divergence arising from the vanishing 
energy denominator problem. The logarithmic infrared divergence is canceled by the self energy 
contribution, Eq. (4.23). Thus it explicitly shows that the logarithmically confining Coulomb 
interaction survives similarity transformation but the associated infrared divergence is canceled 
by self energy contribution. 

Next we look at the most singular term in the numerator in the limit x — > y. In this limit we 
keep only the leading term in the denominator and we get 



-4- 



° f I dy I dq ^° 



e(x-y) + 6(y-x) 



(4.31) 



'2(2^) 2 ~ 7 J " J J *** rzv "^ JaiJ \ (x-y) 2 ' (y-x) 2 
Lastly we look at the rest of the terms in the instantaneous exchange. Since we are interested 
only in the singularity structure, we keep only the leading term in the denominator and we get 



-2- 



2(2n)' 



C f J dy J dq faij 



(k-qf 



x 



6(x-y) \q 2 (\-2y) k 2 (l-2x)i 0(y-x) [k 2 (l - 2x) q 2 (\-2y) 
x — y L y(l-y) x(l—x)l y—x i x(l—x) y(l — y) 



(4.32) 



Transverse gluon exchange 

First, consider the most singular terms. 

Keeping only the most singular terms, the gluon exchange contribution is 

„2 



2(2^) 2 
0(x-8-y) 



Cf J dy J dq y 2 (y,q) x 



x-y 

„2 



gojik + gcijk , Ak-q) 2 
1 \x-y) 2 

{x-y) 



+ 



0(y-x-S) 



y-x 



ga jik + goijk ^_ 4 ^ (q-k) 2 



(q-k) 2 

(y-x) 



(y-xY 



2(2n) 2 

8(x-8-y) 
(x-y) 2 

Explicitly, for x > y, 



° f I dy I dq X 

9(y-x-8) 



8 a jik gaijk 



+ ■ 



(y-xf 



8a jik goijk 



(4.33) 



goijk = e(o 2 -M 2 )e(M 2 k -Ml)e(M 2 k -o 2 ), 
sajik = e(o 2 -M 2 ) e(M 2 k -M 2 k ) e(M? k -o 2 ). 
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(4.34) 



We are interested in the situation x near y and i near j. Then 9 (M 2 k — M 2 k ) = | = 9 (M 2 k — M 2 k ) 
and 9{a 2 —M 2 ) = 1. Then the gluon exchange contribution is 



(x-y) 2 ) V 



(4.35) 



where we have used 0(x) = 1 — 0(— x). Combining with the most singular part of the instanta- 
neous contribution given in Eq. (4.31) we arrive at 

.2 



4 2(ijr)2 Cf J dy J dq 



2{2k) 
0(x-8-y) 

(x-y) 2 



X 



e a 2 



x — v 



+ 



{y-xf 



9 a 2 



y — x 



. (4.36) 



+ . 

For convenience we change variables. For x > y, we put x — y = pf and k — q = p and for 
y > x, we put j — x = and q — k = p l where P + is the total longitudinal momentum. Thus 
we arrive at 



g 



2 C f P+ 



2(2n) 

+ [dp f 
J Jp 



{ P X ) 2 P+ 



p+s 

Consider the Fourier transform 

,2 



P+(l-x) n + 1 / (n l \ 2 P^ 



V(x-,x L ) = -4 



g 



2{2k 



a c f p+ 



P+x 



p+ 



dp 



(P + ) 



(4.37) 



P+S (p 



+ \2 



/~^~Pmax , £ _i_ — -11 { ^ 

dp 1 ei p x -' px 9 I a 2 
~P max V 



dp+ ^ dplMx -_ ip <^ W 



p+s (p + ) 2 - 1 -" 1 



(p 1 ) 2 p + \ 
p + J 



(4.38) 



where p] nax = G\J p%. We are interested in the behavior of V (x , x l ) for large x ,x l . For large 
x~, nonnegligible contribution to the integral comes from the region q + < rpj. For large x l , we 



need pln ax x l to be small, i.e., {p l max ) 2 < i.e., p + < ^2^2 • Thus we have the requirements 
p + < p+ < . 2 - We make the approximations 

\X I (X J CT 



~^Pmax . 11 



Pmax 



2 pL 



(4.39) 
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and e^ q+x 



For large x , we have p + < -ppy < , > the upper limit of p + integral is cut off by pp 
Adding the contributions from both the integrals (which are equal), for large x~ , we have 



V{x~,x l ) « 32 



2(2tt) 



1 



(4.40) 



Thus for large the similarity factors have produced a square root potential but it is also 
infrared singular. 



For large x l the upper limit of p + integral is cut off by , z and we get, 

(X ) (7 



V(x~,x l ) « 32 



2(2^) 



1 



(4.41) 



For large x 1 , similarity factors have produced a linear confining potential which is also infrared 
singular. We note that the rotational symmetry is violated in the finite part of the potential. In 

„2 < 

both cases, however, the infrared singular part is — 32 2 (2iz) 2 a ~^rs wn i cn * s exactly canceled 
by the infrared contribution generated by similarity transformation from self energy, Eq. (4.26). 
Lastly, we consider the terms that go like -j— . Keeping only the leading term in the energy 

x y 

denominator, we have, 

-°^-^i c f j d y jdqwiiy.q) 



x 



2(27T) : 

0(x-y) ggjik + gaijk [ k z (l-2x) _ q 2 (l-2y) 
x-y (k-q) 2 I x(l-x) y(l-y) 



y-x (q-k) 2 
With the step function cut off we have 



0(y-x) gajik + gaijk \ q 2 {\-2y) _ k 2 (l-2x) 

y(l-y) x(l-x) J 



(4.42) 



gajik + goijk « faij (AM| - O 2 ) . 



(4.43) 



Then, combining Eq. (4.32) and Eq. (4.42) for the sub-leading divergences, we get, 
„2 



-2 



2(27t) : 



Cf J dy J dq faij \l/ 2 (y,q) 
'6(x-y) 1 



x 



x-y {k-qf 



l a 2_( k -^ 2 \ [k 2 (l-2x) q 2 {\-2y) 



1 



Bjy-x) 
y-x (q-k) : 



6 a 



x-y J I x(l-x) y(l-y) J 



2 {q-k) 2 \\q 2 {\-2y) k 2 (l-2x) 



y-x 



y(l-y) x(l-x) 



(4.44) 
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Taking the Fourier transform of this interaction, a straightforward calculation shows that no 
log 8 divergence arise from this term. 

Summary of divergence analysis 

The logarithmic confining Coulomb interaction of (2+1) dimensions is unaffected by simi- 
larity transformation and is still affected by a logarithmic divergence which is however canceled 
by a logarithmic divergence from self energy contribution. Similarity transformation leads to a 
non-cancellation of the most singular ( ,_ y p ) term between instantaneous and transverse gluon 
interaction terms. This leads to a linear confining interaction for large transverse seperations 
and a square root confining interaction for large longitudinal separations. The confining inter- 
actions generated by similarity transformations violates rotational symmetry in lowest order of 
perturbation theory and needs higher order calculations to see if restoration of the symmety oc- 
curs. However, non-cancellation also leads to -j= divergences where 8 is the cutoff on \ x — y \. 
This divergence is cancelled by new contributions from self energy generated by similarity 
transformation. The subleading singular -3- terms do not lead to any divergence in 8. 

x y 

4.5 Numerical studies 

The integral equation is converted in to a matrix equation using Gaussian Quadrature. The 
matrix is numerically diagonalized using standard LAPACK routines [20]. We follow the same 
procedure as what we adopted for study with Bloch effective theory (see Chapter 3) and the 
details of numerical procedure are discussed in Appendix C. With the exponential form and 
the step function form of the similarity factor, the integral over the scale in the definition of g a 
factors in Eq. (4.3) can be performed analytically as shown in Sec. 4.4. For parametrization II, 
we perform the integration numerically using n s quadrature points. 

The first question we address is the cancellation of divergences which are of two types: 
(1) the In 8 divergence in the self energy and Coulomb interaction which has its source in the 
vanishing energy denominator problem that survives the similarity transformation and (2) -7= 
divergences in the self energy and gluon exchange generated by the similarity transformation. 
Here I should remind the reader once again that the main motivation of our work is not to 
fit data but to assess the strengths and weaknesses of SRG scheme over Bloch perturbation 
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theory and comparative study of different choices for similarity factor. In Table 4.1 we present 



8 


Parametrization I 


0.1 


4.89535 


4.90359 


4.90420 


4.90420 


4.90482 


0.01 


5.62612 


6.38083 


6.82963 


7.20037 


7.36414 


0.001 


5.68417 


6.42147 


6.90879 


7.30609 


7.64650 


0.0001 


5.68432 


6.42148 


6.90909 


7.30611 


7.64677 


0.00001 


5.68432 


6.42148 


6.90909 


7.30611 


7.64677 


8 


Parametrization II 


0.1 


4.55364 


4.55668 


4.55668 


4.55668 


4.55669 


0.01 


4.86066 


5.33491 


5.49693 


5.59838 


5.79111 


0.001 


4.87607 


5.35671 


5.59226 


5.64476 


5.88613 


0.0001 


4.87604 


5.35671 


5.59236 


5.64477 


5.88615 


0.00001 


4.87604 


5.35671 


5.59236 


5.64477 


5.88615 


8 


Parametrization III 


0.1 


5.12410 


5.13039 


5.13101 


5.13101 


5.13754 


0.01 


6.02600 


6.94326 


7.50445 


7.98054 


8.39927 


0.001 


6.00968 


6.97160 


7.55376 


8.07749 


8.49199 


0.0001 


5.96636 


6.97160 


7.51814 


8.07751 


8.46524 


0.00001 


5.96636 


6.97160 


7.51814 


8.07751 


8.46524 



Table 4.1: Variation with 8 of the first five eigenvalues of the full Hamiltonian (excluding 
the less significant imaginary term). The parameters are m = 1.0, g = 0.6, ri\ = 58, ri2 = 
58, e = 0.00001, A = 20.0, o = 4.0, (w = 0.1, n g = 2 and n s = 500 (for a integration) in 
parametrization II) 

the 8 independence of the first five eigenvalues of the Hamiltonian for g=0.6. Results are 
presented for three parametrizations of the similarity factor, namely, the exponential form, the 
form proposed by St. Glazek and the step function form used in our analytical studies. It is clear 
that the Gaussian Quadrature effectively achieves the cancellation of 8 divergences. Recall that 
in the study of the same problem using Bloch approach in Chapter 3, negative eigenvalues 
appeared for g=0.6 when 8 was sufficiently small (for example, 0.001) which was caused by 
the vanishing energy denominator problem. Our results in the similarity approach for the same 
coupling shows that this problem is absent in the latter approach. 

Next we study the convergence of eigenvalues with quadrature points. In Table 4.2 we 
present the results for all three parametrizations of the similarity factor for the coupling g=0.2 
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with the transverse space discretized using k = ^tan^ where w's are the quadrature points (see 
Appendix C). The table show that for m = 1, convergence is rather slow for all three choices 
of the similarity factor compared to the results in Bloch approach. Among the three choices, 
parameterization II shows better convergence. 



(A) 




Figure 4.2: The ground state wavefunction for different choices of the similarity factor using 
the parametrization k = ^tan(q7t/2) for transverse momentum grid and for n\ = 40, n,2 = 80, 
m=1.0, g = 0.2, o = 4.0 and e = 8 = 10~ 5 , k = 10.0 as a function of x and k. (A) Parametriza- 
tion I, (B) Parametrization II with n g = 2, uq = 0.1, n s = 500, (C) Parametrization III. 

Let us now discuss the nature of low lying levels and wavefunctions. First we show the 
ground state wavefunctions for all three similarity factors for a given choice of parameters in 
Fig. 4.2. As is anticipated, step function choice produces a non-smooth wavefunction. For 
parameterizations I and II, the wavefunctions show some structure near x = 0.5. From our 
previous experience with calculations in the Bloch formalism, we believe that the structures 
indicate poor convergence with the number of grid points. 
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n\ 


«2 


Parametrization I 


10 


10 


4.320 


4.353 


4.357 


4.357 


4.361 


20 


20 


4.375 


4.442 


4.484 


4.484 


4.485 


20 


30 


4.398 


4.482 


4.546 


4.583 


4.610 


20 


40 


4.411 


4.502 


4.570 


4.615 


4.656 


30 


40 


4.412 


4.503 


4.570 


4.615 


4.655 


40 


50 


4.420 


4.515 


4.585 


4.634 


4.678 


40 


60 


4.426 


4.524 


4.594 


4.645 


4.692 


40 


80 


4.434 


4.535 


4.607 


4.661 


4.709 


n\ 


n 2 


Parametrization II 


10 


10 


4.163 


4.194 


4.194 


4.194 


4.203 


20 


20 


4.186 


4.244 


4.276 


4.276 


4.277 


20 


30 


4.192 


4.256 


4.296 


4.323 


4.329 


20 


40 


4.195 


4.262 


4.304 


4.335 


4.344 


30 


40 


4.195 


4.262 


4.304 


4.335 


4.344 


40 


50 


4.197 


4.266 


4.308 


4.341 


4.353 


40 


60 


4.199 


4.268 


4.311 


4.345 


4.359 


40 


80 


4.201 


4.272 


4.315 


4.350 


4.367 


n\ 


n 2 


Parametrization III 


10 


10 


4.360 


4.379 


4.381 


4.381 


4.391 


20 


20 


4.469 


4.533 


4.572 


4.572 


4.572 


20 


30 


4.503 


4.604 


4.674 


4.721 


4.749 


20 


40 


4.520 


4.632 


4.703 


4.768 


4.810 


30 


40 


4.528 


4.636 


4.714 


4.768 


4.811 


40 


50 


4.542 


4.657 


4.736 


4.797 


4.848 


40 


60 


4.548 


4.668 


4.748 


4.813 


4.866 


40 


80 


4.556 


4.683 


4.764 


4.833 


4.889 



Table 4.2: Convergence of eigenvalues with n\ and n 2 for the parametrization k = ^tan(qn/2). 
The parameters are m=1.0, g=0.2, £ = 0.00001, 8 = 0.00001, k = 10.0, a = 4.0, (w =0.1, n g =2 
and the number of quadrature points n s = 500 for o integration for parametrization II). 
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Now consider the structure of low lying levels. Recall that in the Bloch formalism discussed 
in the previous chapter, the ordering of levels was / = 0, 1,0, ... corresponding to logarithmic 
potential in the nonrelativistic limit (see Appendix D). In the presence of effective interactions 




Figure 4.3: The wavef unctions corresponding to the lowest four eigenvalues as a function of 
x and k for parametrization II. The parameters are as in FIG. 4.2. (A) Lowest state, (B) first 
excited state, (C) second excited state, (D) third excited state. 



generated by the similarity transformation, obviously the level ordering changes. Now we have 
additional confining interactions which, however, act differently in longitudinal and transverse 
directions. From our analytic calculation in Sec. 4.4 we know that the confining potential 
generated by similarity transformation is linear in x l and square root in x~ . Thus a node in x l 
costs more energy than a node in x~ and thus the states with nodes in x l (i.e., in k) will be of 
higher energy compared to states with nodes in x~ . This feature is well manifested in Fig. 4.3 
where the wavefunctions for the first four low lying levels are presented for parametrization II. 
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There is extra freedom in parametrization II due to the presence of £M^- in the definition of 
Uaij, Eq. (4.18). For zero transverse momentum of constituents, EM? has the minimum value 
8m 2 . Thus relative insensitivity of parametrization to o in parametrization II for small values 
of (7 may be due to this factor. When we consider the heavy fermion mass limit, presence of 
8m 2 in u G ij enhances the effect of similarity factor. In Fig. 4.4 we present the wavefunctions 
corresponding to first four levels for parametrization II with 8m 2 subtracted from £M ; 2 in u G ij- 
From Figs. 4.3 and 4.4, note that the fourth level is different for parametrization II with and 
without 8m 2 in u a ij. 




By suitable choice of parameters we can study the interplay of rotationally symmetric loga- 
rithmically confining interaction and effective interactions generated by similarity transforma- 
tion. Since for a given coupling constant g, strength of the logarithmic interaction and similarity 
generated interactions are determined by m and o respectively, for m > > a we should recover 
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the Bloch spectrum (Fig. 3.3). Upto what levels the recovery occurs, of course depends on the 
exact value of m and the energy scale a at which the effective Hamiltonian is constructed. As 
we have already observed, for parametrization II this will happen only if 8m 2 is subtracted from 
EM 2 . For this case, we present the first four levels for m = 10.0 and a = 4.0 in Fig. 4.5 which 
clearly shows the level spacing corresponding to the Bloch spectrum presented in Fig. 3.3. 




Figure 4.5: The wavef unctions corresponding to the lowest four eigenvalues as a function of 
x and k for large fermion mass with parametrization II with 8m 2 subtracted from EM ; 2 . The 
parameters are m = 10.0 g = 0.2, e = 8 = 0.00001, k = 10.0, n g = 2, uq = 0.1, n\ = 40, 
«2 = 80. (A) Lowest state, (B) first excited state, (C) second excited state, (D) third excited 
state. 

Finally, we discuss the sensitivity of the spectra to the similarity scale o. Ideally the low ly- 
ing energy levels should be insensitive to a. However we have calculated the effective Hamilto- 
nian to only order g 2 and we expect significant sensitivity to o. In Tables 4.3 and 4.4 we present 
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Eigenvalues (M 2 ) 


a 


Parametrization I 


2.0 


4.214 


4.285 


4.329 


4.365 


4.393 


4.0 


4.434 


4.535 


4.607 


4.661 


4.709 


6.0 


4.701 


4.821 


4.914 


4.980 


5.043 


a 


Parametrization II 


2.0 


4.157 


4.214 


4.248 


4.260 


4.276 


4.0 


4.201 


4.272 


4.315 


4.350 


4.367 


6.0 


4.266 


4.350 


4.404 


4.446 


4.483 


a 


Parametrization III 


2.0 


4.254 


4.346 


4.395 


4.443 


4.477 


4.0 


4.556 


4.683 


4.764 


4.833 


4.889 


6.0 


4.927 


5.075 


5.179 


5.262 


5.333 



Table 4.3: Variation with o of the full Hamiltonian (excluding the imaginary term). The pa- 
rameters are m = 1.0, g = 0.2, nl = 40, nl = 80, £ = 0.00001, k = 10.0, (k = ±tan(qn /!)), 
8 = 0.00001, (wo = 0.1, n g = 2, and n s = 500 for o integration for parameterization II). 



the lowest five eigenvalues for all three parametrizations of the similarity factor for g = 0.2 and 
g = 0.6 respectively. As expected o dependence is greater for larger value of g. 

Among the three parametrizations, the paramterization II is least sensitive to o. In order to 
check whether this behaviour is due to the presence of 8m 2 in EM ; 2 in the definition of u a ij we 
present the results in Table 4.5 for parametrization II with 8m 2 subtracted from EM ; 2 . It is clear 
that sensitivity to o is still considerably less compared to the other two parametrizations. This 
may be due to the fact that EM?- is added to a 2 in the definition of o. 

Note that in parametrization II sensitivity to a is controlled also by additional parameters uq 
and n g . The sharpness of the cutoff in this paramatrization depends on n g . The cutoff becomes 
sharper as one increases the value of n g . o dependence of the eigenvalues for two different 
n g are prestented in Tables 4.5 and 4.6. The other adjustable parameter uq is a samll number 
(uq « 1) and the sensitivity to uq of the lowest five eigenvalue for the coupling g = 0.6 is 
presented in Table 4.7. 
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Eigenvalues (M 2 ) 


a 


Parametrization I 


2.0 


4.941 


5.418 


5.708 


5.713 


5.950 


4.0 


5.684 


6.421 


6.909 


7.306 


7.647 


6.0 


6.678 


7.617 


8.274 


8.801 


9.267 


a 


Parametrization II 


2.0 


4.769 


5.152 


5.226 


5.372 


5.485 


4.0 


4.876 


5.357 


5.592 


5.645 


5.886 


6.0 


5.071 


5.655 


6.016 


6.159 


6.316 


a 


Parametrization III 


2.0 


4.888 


5.589 


5.882 


6.120 


6.263 


4.0 


5.966 


6.972 


7.518 


8.077 


8.465 


6.0 


7.359 


8.603 


9.360 


10.083 


10.621 



Table 4.4: Variation with o of the full hamiltonian (excluding the imaginary term). The pa- 
rameters are m = 1.0, g = 0.6, nl = 58, nl = 58, £ = 0.00001, A = 20.0, (k = {l _ q q i™ A+m ), 
8 = 0.00001, (wo = 0.1, n g = 2,and n s = 500 for o integration for parameterization II). 

4.6 Summary and Discussion 

From our results of Chapter 3, we know that the attempt to solve (2+1) dimensional gauge 
theories using the Bloch effective Hamiltonian is unseccessful due to problem of uncancelled 
infrared divergences. They arise out of vanishing energy denominators and a more sophisticated 
tool is necessary to handle the bound state problems in QCD. Similarity renormalization formal- 
ism attempts to solve the bound state problem in a two step process. At the first step, coupling 
between low and high energy degrees of freedom are integrated out and ultraviolet renormal- 
ization carried out perturbatively. In the second step, the effective Hamiltonian is diagonalized 
non-perturbatively. 

Here I briefly summarize the main points addressed in this chapter. In order to have a bet- 
ter understanding of the numerical results, we have performed analytical calculations with step 
function form for the similarity factor. Many interesting results emerge from our analytical 
calculations. First of all, it is shown that due to the presence of instantaneous interactions in 
gauge theories on the light front, the logarithmic infrared divergence that appeared in the Bloch 
formalism persists in two places, namely a part of the self energy contribution and the Coloumb 
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Eigenvalues (M 2 ) 


8 


a 


Parametrization II 




2.0 


4.126 


4.166 


4.176 


4.189 


4.201 


U.z 


4.U 


4.172 


4.235 


4.273 


4.298 


4.304 




6.0 


4.241 


4.321 


4.371 


4.411 


4.445 




2.0 


4.739 


5.003 


5.020 


5.134 


5.193 


0.6 


4.0 


4.802 


5.222 


5.350 


5.470 


5.626 




6.0 


4.993 


5.541 


5.876 


5.940 


6.156 



Table 4.5: Variation with o of the full Hamiltonian (excluding the imaginary term) after sub- 
tracting 8m 2 from £M? in the definition of u C ij. The parameters are m = 1.0, £ = 0.00001, 
8 = 0.00001, u = 0.1, n g = 2, and n s = 500 

1) for g = 0.2, k = \tan{q%jl) with k = 10.0 and nl = 40, nl = 80 

2) for g = 0.6, k = ffi? with A = 20.0, and nl = nl = 58. 

' ° (1— q L A+m 



8 


a 


M 2 (Parameterization II) 




2.0 


4.781 


5.049 


5.051 


5.161 


5.221 


0.6 


4.0 


4.843 


5.235 


5.391 


5.464 


5.639 




6.0 


5.030 


5.538 


5.848 


5.956 


6.106 



Table 4.6: Variation with o of the full Hamiltonian(excluding the imaginary term) after subtract- 
ing the 8m 2 term from £M 2 - in the definition of u G ij- The parameters are m = 1.0, £ = 0.00001, 

8 = 0.00001, u = 0.1, n g = 1, and n s = 500 g = 0.6, k = (1 _^ +m with A = 20.0, and 

n \ = n 2 = 58. 



interaction that gives rise to the logarithmically confining potential in the nonrelativistic limit. 
However the terms that persist are precisely those that produce a cancellation of resulting in- 
frared divergences in the bound state equation. The rest of the infrared problem that appeared in 
the Bloch formalism due to the vanishing energy denomnator problem is absent in the similarity 
formalism. 

Similarity transformation however prevents the cancellation of the most severe j^zyp. sm " 
gularity between instantaneous gluon exchange and transverse gluon exchange interactions and 
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8 


U 


M 2 (Parameterization II) 




0.2 


4.990 


5.548 


5.887 


5.969 


6.171 




0.1 


4.802 


5.222 


5.350 


5.469 


5.626 


0.6 


0.05 


4.746 


5.076 


5.091 


5.253 


5.313 




0.01 


4.813 


5.003 


5.029 


5.069 


5.130 



Table 4.7: Variation with uq of the full Hamiltonian(excluding the imaginary term) after sub- 
tracting the 8m 2 term from EM?- in the definition of u C ij. The parameters are m = 1.0, 

£ = 0.00001, 8 = 0.00001, o = 4.0, n g =2, and n s = 500 g = 0.6, k = (1 _^ +>H with A = 20.0, 
and nl = n2 = 58. 

produces divergences in the self energy and gluon exchange contributions which cancel be- 
tween the two in the bound state equation. The resulting effective interaction between the quark 
and antiquark grows linearly with large transverse separation but grows only with the square 
root of the longitudinal separation. This produces severe violations of rotational symmetry in 
the bound state spectrum. We have also verified that no In 8 divergence results from the -3- 

x y 

singularity in the self energy and gluon exchange contributions. 

In the Glazek- Wilson formalism the exact form of the similarity factor f a is left unspecified. 
In the literature an exponential form has been used in numerical calculations [15]. For analyt- 
ical calculations it is convenient to choose a step function even though it is well known that 
it is not suitable for quantitative calculations [5]. There is also a proposal due to Stan Glazek 
which has two extra free parameters. We have tested all three parametrizations. Our numeri- 
cal results indeed show that step function choice always produces non-smooth wavefunctions. 
Parameterization II costs us an extra integration to be performed numerically but convergence 
is slightly better for small g compared to exponential form. All three parametrizations produce 
violations of rotational symmetry even for small g. When an exponential form is used in the 
Glazek- Wilson formalism, the resulting effective Hamiltonian differs from the Wegner form 
only by an overall factor that restricts large energy diffrences between initial and final states. 
Numerically we have found this factor to be insignificant. 

We have studied the sensitivity of the low lying eigenvalues to the similarity scale a. Since 
the effective Hamiltonian is calculated only to order g 2 results do show sensitivity to o. Among 
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the three parametrizations the form II is least sensitive to o due to the functional form chosen. 
We have also studied the sensitivity of eigenvalues to the parameters w and n g . 

The bound state equation has three parameters m, g 2 and o with dimension of mass. The 
strength of the logarithmically confining interaction is determined by m and the strength of the 
rotational symmetry violating effective interactions generated by similarity transformation is 
determined by o. For a given g we expect the former to dominate over the latter for m >> a. 
An examination of low lying eigenvalues and corresponding wavefunctions show that this is 
borne out by our numerical calculations. 

A major problem in the calculations is the slow convergence. Compared to the Bloch for- 
malism, in calculations with the similarity formalism, various factors may contribute to this 
problem with the Gauss quadrature points. One important factor is the presence of linear and 
square root confining interactions generated by the similarity transformation. It is well known 
that such interactions are highly singular in momentum space. Another factor is the presence of 

divergences, the cancellation of which is achieved numerically. It is of interest to carry out 
the same calculations with numerical procedures other than the Gauss quadrature. However, 
one should note that calculations in (3+1) dimensions employing basis functions and splines 
have also yielded [15] wavefunctions which show non-smooth structures. 

An undesirable result of the similarity transformation carried out in perturbation theory is 
the violation of rotational symmetry. Our results show that this violation persists at all values 
of g for m = 1. Such a violation was also observed in (3+1) dimensions. In that case the 
functional form of the logarithmic potential generated by similarity transformation is the same 
in longitudinal and transverse directions but the coefficients differ by a factor of two. Same 
mechanism in (2+1) dimensions makes even the functional forms different. The important 
questions are whether the confining interactions generated by the similarity transformation are 
an artifact of the lowest order approximation and if they are not, then, whether the violation of 
rotational symmetry will diminish with higher order corrections to the effective Hamiltonian. 
Recall that matrix element between low and high energy degrees of freedom has been integrated 
out and the effective low energy Hamiltonian determined only to order g 2 . A clear answer will 
emerge only after the determination of the effective Hamiltonian to fourth order in the coupling. 
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CHAPTER 5 



Fermion Formulation on a Light-Front Transverse Lattice 

5.1 Introduction 

In the previous two chapters we have discussed the meson bound state problem in (2+1) 
dimensions in the light-front framework using Bloch and similarity effective Hamiltonians and 
have seen that light-front framework provides the opportunity for non-perturbative studies in 
Hamiltonian formalism. In this chapter we introduce another approach in light-front Hamilto- 
nian framework. 

Till date, the most practiced non-perturbative technique is the lattice gauge theory [1]. Us- 
ing path integral formalism in Euclidean space, with no gauge fixing one calculates the n-point 
Green functions, but no direct informations about the bound state wavef unctions are accessible 
in this approach. Hamiltonian formalism provides bound state wavefunctions in a straight- 
forward way. Since boosts are kinemetical, one can have frame independent description of the 
bound state wavefunctions in light-front framework. One extra advantage of lattice gauge the- 
ory over light-front formalism is full gauge invariance. Light front Hamiltonian formulation of 
transverse lattice QCD [2, 3] is an optimum combination of lattice gauge theory and light-front 
QCD. It uses the power of light-front Hamiltonian formalism to produce the boost invariant 
wavefunctions and the advantage of having gauge invariant ultraviolet cutoff from lattice gauge 
theory. With the gauge choice A + = A + A 3 = and the elimination of the constrained vari- 
able A~ =A°-A 3 , it uses minimal gauge degrees of freedom in a manifestly gauge invariant 
formulation exploiting the residual gauge symmetry in this gauge. 

For Hamiltonian formalism light-front time (x + ) is kept continuous. Due to the constraint 
equations in fermion and gauge degrees of freedom in the light-front, non-locality comes in the 
longitudinal direction. Also there is no ultraviolet divergences coming from small x~ and hence 
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x~ is not latticized. Ultraviolet divergences come only from small transverse separations which 
we want to regulate with gauge invariant cutoff. Thus, the transverse plane (jr 1 = (x l ,x 2 )) is 
discretized on a square lattice. This defines the light-front transverse lattice (LFTL) in (3+1) 
dimensions. It is a promising and developing tool for non-perturbative investigations of QCD. 
So far encouraging results have been obtained in the pure gauge sector [4] and in the meson 
sector with particle number truncation (for a recent review see, Ref. [5]). 

It is well known that fermions on the lattice pose challenging problems due to the dou- 
bling phenomenon. Light-front formulation of field theory has its own peculiarities concerning 
fermions because of the presence of a constraint equation. As an example, the usual chiral trans- 
formation on the four component fermion field is incompatible with the constraint equation for 
nonzero fermion mass [6]. There have been previous studies of fermions on the transverse lat- 
tice [7, 8, 9, 10] in different contexts. But properties and origin of species doubling of fermion 
on light-front transverse lattice were not studied with proper care and desired details. So, before 
embarking on any QCD calculation on transverse lattice, it is better to understand the fermions 
on it in detail. I devote this chapter to discuss different ways of formulating fermion on a light- 
front transverse lattice and related issues such as absence or origin of doubers, different ways 
of removing doubers, relevant symmetry on light-front transverse lattice and so on [11]. 

As we shall see later in this chapter, the presence of the constraint equation in light front field 
theory allows different methods to put fermions on a transverse lattice. It is worthwhile to study 
all the different methods in order to examine their strengths and weaknesses. Here I should 
also mention two important points. One, our ultimate aim is to calculate QCD observables on 
LFTL where we deal with a Hamiltonian acting on a Fock space. For a reasonable size of Fock 
space, computing limitations will force us to be in a reasonably small lattice volume when we 
deal with realistic problems. The second point is that the currently practiced version of the 
transverse lattice gauge theory uses linear link variables [2, 3] justified on a coarse lattice and 
recovering continuum physics is nontrivial. Hence it is very important to carry out detailed 
numerical investigation of all the possible ways of formulating fermions on LFTL with proper 
attention to the finite volume effects. 

In one of the approaches of treating fermions on the light front transverse lattice, we main- 
tain as much transverse locality as possible on the lattice by using forward and backward lattice 
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derivatives without spoiling the hermiticity of the Hamiltonian. In this case doublers are not 
present and the helicity flip term proportional to the fermion mass in the full light front QCD 
becomes an irrelevant term in the free field limit. Thus in finite volume, depending on the 
boundary condition used, the two helicity states of the fermion may not be degenerate in the 
free field limit. However, we find that in the infinite volume limit the degeneracy is restored 
irrespective of the boundary condition. 

In the second approach [8], symmetric derivatives are used which results in a Hamiltonian 
with only next to nearest neighbor interaction when we take the free field limit. As a conse- 
quence even and odd lattice sites decouple and the fermions live independently of each other on 
the two sets of sites. As a result we get four species of fermions on a two dimensional lattice as 
excitations around zero transverse momentum. Note that this is quite different from what one 
gets in the conventional Euclidean lattice theory when one uses symmetric derivatives. In that 
case, doublers have at least one momentum component near the edge of the Brillouin zone. The 
doublers can be removed in more than one way. We also study the staggered fermion formula- 
tion on the light front transverse lattice to eliminate two doublers and reinterpret the remaining 
two as two flavors. In this light front staggered fermion formulation, there is no flavor mixing 
in free field limit. But, in QCD, we get irrelevant flavor mixing terms. An alternative which 
removes doubling completely is to add the conventional Wilson term which generates many ir- 
relevant interactions on the transverse lattice. Among them, the helicity flip interactions vanish 
but the helicity non flip interactions survive in the free field limit. 

5.2 Hamiltonian with forward and backward derivatives 
5.2.1 Construction 

In this section we propose to use different lattice derivatives for dynamical and constrained 
fermion field components on the transverse lattice in such a way that sacred Hermiticity of the 
Hamiltonian is preserved. The starting point of our discussion is the fermionic part of the QCD 
Lagrangian density 

^ f = xj/(ifD^-m)xjf (5.1) 
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with z'Z)^ = id^ —gA^ 1 . Since fermions on LFTL are our main concern in this chapter, we omit 
the pure gauge part from QCD Lagrangian. 

Moving to the light front coordinates we impose the light-cone gauge A+ = and introduce 
the transverse lattice by discretizing the transverse plane on a square lattice with lattice spacing 
a. Now, S£f on the LFTL can be written as 

& f = \i/ +Jf (id--gA-)\i/ + + \i/--id + y- 

-my'V^-m^Vv"- (5-2) 
Here r = 1 , 2 and Dp b is the forward/backward covariant lattice derivative defined as 

D(t] (x) = i [U r (x) T) (x + ar) - T) (x)] (5.3) 

and 

D b r r] (x) = - [tj (x) - u}(x - ar) r\ (x - af )] , (5.4) 

ft 

where a is the lattice constant and r is unit vector in the direction r = 1,2 and D{ = —D b r . 
U r (x) is the group valued lattice gauge field with the property C//(x) = U- r (x + ar). In the 
weak coupling limit 

U r (x)^e l8aAr{x+at/2 \ (5.5) 

For notational convenience we suppress x~ in the arguments of the fields. 

Our goal here is to write the most local lattice derivative. That is why, instead of using 
the symmetric lattice derivative, in the above we have used the forward and backward lattice 
derivatives. However, the Hermiticity of the Lagrangian (Hamiltonian) requires that if one of 
the covariant lattice derivatives appearing in Eq. (5.2) is the forward derivative, the other has to 
be the backward derivative or vice versa. 

The constraint equation is 

id + \j/- = (ia r D( + fm)\if + . (5.6) 
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Eliminating the constrained field component y in terms of dynamical field we obtain 



& f = \if + \id- -gA~)y + -m^ ++ yV~ -iy^ a r D h r x\f~ 
= X]/ + \id--gA-) ¥ + 

- y + \ia r D b r + 7°m] [ia s D{ + fm] y + . 



(5.7) 



The dynamical field can essentially be represented by two components [12] such that 



^) 



r\{x ,x L ) 




(5.8) 



where r\ is a two component field. Finally going over to the two component fields r\, the 
Lagrangian density can be written as 



se f = 7]\id-- g A-)j] 



- 7] t [ic r D b r - im] [id s D{, + im] T\ 



(5.9) 



&\ = 02 and 62 = — 0\ where a ; are Pauli spin matrices. Writing explicitly in terms of the link 
variables, the Lagrangian density is 

&f = lH^)(id--gA-)ri(x)-m 2 ri\x)^ T ri(x) 

- mT) ^ (x) V O r - -^rr \u r (x) T) (x + at) - 77 (x) 



m£ ^(x + ar)U?(x)-ri\x) a r -— tj(x) 



1 1 



-l£[r 7 t (x + ar)^(x)-r 7 t (x)]d r ^d,[f/. v (x)r 7 (x + as)-r 7 (x)]. (5.10) 



a id + 
1 

In the free limit the fermionic part of the Hamiltonian becomes 
Pj b = I dx-a 2 ^^ 3 

- ^2 ^ f ( x ) E T^T ( x + af ) - 2r 7 ( x ) + ^ ( x - flf ) 

I 1 
+ ^2 V f ( x ) X>™<^) [tj (x + ar) - 2 77 (x) + 77 (x - ar)} 
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(5.11) 



In order to get Eq. (5. 1 1), we have assumed infinite transverse lattice and accordingly have used 
shifting of lattice points which is equivalent of neglecting surface terms. The positive sign in 
front of the last term would change if we had switched forward and backward derivatives. One 
should note that the cross term coming from the last line of Eq. (5.10) survives in the free field 
limit. Explicitly, the term is 

l£[r 7 t (x + ar)-r 7 (x)]d r ^ T d s .[r 7 (x + as)-r 7 (x)]. (5.12) 

CI i 10 

It produces extra helicity nonflip hoppings in the transverse plane. In the continuum limit {a — > 
0) this term does not survive and one can recover the right continuum limit without this cross 
term. If one demands hypercubic (square) symmetry of the transverse lattice, then this cross 
term (Eq. (5.12)) vanishes. More detailed discussion about this term is provided in Appendix 
F. 

Because of the presence of 6 r , the last term of Eq. (5.11) couples fermions of opposite 
helicities. Note that it is also linear in mass. Such a helicity flip linear mass term is typical 
in continuum light-front QCD. Here in free transverse lattice theory this term arises from the 
interference of the first order derivative term and the mass term, due to the constraint equation. 
This is in contrast to the conventional lattice (see Appendix G) where no helicity flip or chirality- 
mixing term arises in the free theory if we use forward and backward lattice derivatives. 

5.2.2 Absence of doubling 

Consider the Fourier transform in transverse space 

!?(*-,*) = / (02^^) < 5 - 13 ) 

where — ^ < &i,&2 < +f • Then the helicity nonflip part of Eq. (5.11) becomes 




(5.14) 



Using 

a 2 jV k -P)' x = {2k) 1 S 2 (k-p) (5.15) 

X 
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we get, 

P «/ = / dx "/(^ 0kt(x_) ^ 0k(x_) 



(5.16) 



7* r \ & r a/2 

where we have defined k r a = k ■ ra. Note that the sine function vanishes at the origin k\ , &2 = 
but does not vanish at the edges of the Brillouin zone k\ ,&2 = ±f • 

~ sin & ~ 

Define fc r = fc r sl " ^ . In the naive continuum limit fc r — ► £ r . 

Now, let us consider the full Hamiltonian (Eq. 5.1 1) including the helicity flip term. In the 
helicity space we have the following matrix structure for P + P (since P is inversely propor- 
tional to the total longitudinal momentum P + , we study the operator P + P ) 

,2 i 4 v _4jn/;-„ 2M , _■ 2 



m 2 + JE, sin 2 k -f - ^ (i sin 2 ^ + sin 2 ^ ) 

^(/ sin 2^_ sin 2^) m 2 + ^ L( _ sin 2^ 

which leads to the eigenvalue equation 



(5.17) 



Jl 2 = m 2 + l£sin 2 ^ ± - j£sin 4 ^. (5.18) 

Third term in the above equation comes from the linear mass helicity flip term. If the mass 
m = 0, then it is obvious from Eq. (5.18) that M 2 = if and only if k\=k2 = 0. For nonzero 
m, one can also in general conclude that M 2 = m 2 only for the case k\=k% = 0. Thus there 
are no fermion doublers in this case (for physical masses am < 1). In the following for specific 
choices of momenta we elaborate on this further. 

If one component of the momentum vanishes, then 

^2 = m 2 + 4 ( l ±m)s . n 2^ (519) 

a a 2 

where k is the non- vanishing momentum component. Thus for am = 1, irrespective of the value 
of k we get ^ 2 = m 2 which is unwanted. In general, for am > 1, M 2 can become negative. It 
is important to recall that physical particles have m < ^ (the lattice cutoff) and hence are free 
from the species doubling on the lattice. With periodic boundary condition (discussed in the 
next subsection), allowed k values are k q a = ±2^+r> with q = 1, 2,3, for 2n + 1 lattice sites 
in each direction. Let k\ = 0. For ma = 1.0, Eq. (5.19) with the minus sign within the bracket 
gives M 2 = m 2 for all values of ki and we get 2(2n+ l)-fold degenerate ground state with 
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eigenvalue m . 

The two spin states (spin up and down) are degenerate for k\ = ki = 0. But if any one (or both) 
of the two transverse momenta is (are) nonzero then the degeneracy is broken on the lattice by 
the spin flip term proportional to m. So the total degeneracy of the lowest states for ma — 1.0 
can be calculated in the following way: (a) k\ = &2 = : Number of states =2 (spin up and 
spin down), (b) k\ = 0, &2 7^ : Number of states =2n and (c) k\ ^ 0, &2 = : Number 
of states =2n. Note that k[ can have 2n nonzero values and there is no spin degeneracy for any 
nonzero k[. So, the total number of degenerate states = 2 + 2n + 2n = 2(2n +1). But if ma ^ 1 
we cannot have m 2 eigenvalue for nonzero k t and we have only two (spin) degenerate states with 
eigenvalue m 2 . Again we see from Eq. (5.19) that if ma > 1, the kinetic energy term becomes 
negative and the eigenvalues go below m 2 . But ma > 1 means m > ^ (ultra violet lattice cutoff) 
and hence unphysical. 

5.2.3 Numerical Investigation 

For numerical investigations we use Discretized Light Cone Quantization (DLCQ) [13] for 
the longitudinal direction (— L <x< +L) and implement antiperiodic boundary condition to 
avoid zero modes. Then, 



T?(*-,x) = 4fI> I [b{l^X)e- M ^l^+d\l^-X)e Mx -l^] (5.20) 



In DLCQ with antiperiodic boundary condition, it is usual to multiply the Hamiltonian P by 
so that H = jP has the dimension of mass squared. 
In DLCQ the Hamiltonian (Eq. 5.1 1) can be written as, 




/=1,3,5,... 



with 



{b(l,x,X),b\l',x',X')} = {d(l,x,X),d\l',x',X')} = SiyS^f d X:X , 



(5.21) 



Hfb =Ho+H h f 



(5.22) 



where, the helicity nonflip part 



#o = t t (/,z,A)M/,z,A)+J t (/,z,A)J(/,z,A)j 
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-EEEEEyZW \b%x,k')b(l,* + at,k) -2b\l,z,X')b(l,z,X) 

*w v 1 111 



z r A X' / 



+^ t (/,z,A / )^(/,z-ar,A)+J t (/,z,A / )^(/,z + ar,A)-2J t (/,z,A / )^(/,z,A) 



z, X')d(l, z — or, <j) 



(5.23) 



and helicity flip term 



H hf 



z r X X' I 



-2b f (l,z,X')b(l,z,X)+b\l,z,X')b(l,z-ar,X)+d t (l,z,X')d(l,z + ar,X) 




(5.24) 



Let us now investigate the Hamiltonian (Eq. 5.22) with two types of boundary conditions 
on the transverse lattice: (1) fixed boundary condition and (2) periodic boundary condition. 
For each transverse direction, we choose 2n + 1 lattice points ranging from — n to +n where 
fermions are allowed to hop. For the study of the fermion spectra on the transverse lattice, the 
longitudinal momentum plays a passive role and for the numerical studies we choose the dimen- 
sionless longitudinal momentum (I) to be unity which is kept fixed. For a given set of lattice 
points in the transverse space we diagonalize the Hamiltonian and compute both eigenvalues 
and eigenfunctions. 

Boundary conditions 

Fixed boundary condition: To implement fixed boundary condition we add two more 
points at the two ends and demand that the fermion remains fixed at these lattice points. Thus 
we consider 2n + 3 lattice points. Let us denote the fermion wavefunction at the location s 
by u(s). We have u(s) ~ sin(s — \)ka with w(l) = u(2n + 3) = 0. Allowed values of k are 
(2n + 2)k p a = p% with p = 1, 2, 3, 2n + 1 and k p = ( 2 n+2) a P- Thus the minimum k p allowed 
is ^ (2n + 2) an d maximum k p allowed is f (2«+2) • ^ or exam pl e > for n = 1 we have k\ = ^, k 2 = 



Periodic boundary condition: Periodic boundary condition identifies the (2n + 2) lattice 
point with the first lattice point. In this case we have the fermion wavefunction u(s) ~ e lska 
with the condition u{s) = u(s + L) where L = 2n + 1. Thus [2n + \)k p a = ±2np so that k p = 



IE ka — IE etc 
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~^ {2n+X)a P' P = 0) 1)2, ...,n. Thus the minimum k p allowed is and the maximum k p allowed 
is (2n+i) a - Forn= l,we have, fc = 0,ki = ±§f, etc. 

Numerical results 

First we discuss the results for Ho given in Eq. (5.23). We diagonalize the Hamiltonian using 
basis states defined at each lattice point in a finite region in the transverse plane. Let us denote 
a general lattice point in the transverse plane by (x^y;). For each choice of n (measure of the 
linear lattice size), we have — n < X[,y[ < +n. Thus for a given n, we have a {In + 1) x {In + 1) 
dimensional matrix for the Hamiltonian. The boundary conditions do have significant effects at 



1.6 r 



1.4 - 



1.2 



2 



6 

n 



10 



Figure 5.1: Convergence of ground state eigenvalue (of forward-backward Hamiltonian) versus 
n(m = 1.0). 

small volumes. For example, a zero transverse momentum fermion at finite n is not allowed with 
fixed boundary condition. But with periodic boundary condition, we can have zero transverse 
momentum fermion for any finite n. With fixed boundary condition, in the infinite volume 
limit, we expect the lowest eigenstate to be the zero transverse momentum fermion with the 
eigenvalue m 2 . In Fig. 5.1 we show the convergence of the lowest eigenvalue as a function of n. 

For a zero transverse momentum fermion, the probability amplitude to be at any transverse 
location should be independent of the transverse location. Thus we expect the eigenfunction 
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Figure 5.2: Eigenfunctions of first three states of the Hamiltonian Hq with lattice points n = 5. 

for such a particle to be a constant. At finite volume, with fixed boundary condition, we do 
get a nodeless wave function which nevertheless is not a constant since it carries some non- 
zero transverse momentum. All the excited states carry non-zero transverse momentum in the 
infinite volume limit. All of them have nodes characteristic of sine waves. The eigenfunctions 
corresponding to the first three eigenvalues are shown in Fig. 5.2 for the case of fixed boundary 
condition. With periodic boundary condition, for any n, we get a zero transverse momentum 
fermion with a fiat wave function. 

Now, we consider the effect of helicity flip term. With fixed boundary condition the lowest 
eigenstate has non vanishing transverse momentum in finite volume. In the absence of helicity 
flip term positive and negative helicity fermions are degenerate. The helicity flip term lifts the 
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Figure 5.3: Spin splitting of the ground state caused by the spin dependent interaction as a 
function of n. Circles represent the eigenvalues without splitting and the vertical bars show the 
magnitude of splitting due to spin flip interaction. 

degeneracy. The splitting is larger for larger transverse momentum. In Fig. 5.3 we present the 
level splitting for the helicity up and down fermions as a function of n. As expected, the level 
splitting vanishes and we get exact degeneracy in the infinite volume limit. For the periodic 
boundary condition, the lowest state has exactly zero transverse momentum and we get two 
degenerate fermions for all n. 

5.3 Hamiltonian with symmetric derivative 
5.3.1 Construction 

Now, let us discuss the other way of putting fermions on a transverse lattice. Here we use 
the symmetric lattice derivative defined by 



D r ^ ± (x) 



1 



[U r (x) y ± (x + of) - U- r (x) y± (x - at )] . 



(5.25) 



2a 
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In place of using forward and backward derivatives in Eq. (5.2), we use the above symmetric 
derivative for all lattice derivatives. Proceeding as in Sec. 5.2, we arrive at the fermionic part of 
the QCD Hamiltonian 



sd 



= Jdx fl 2 £m 2 T7 + (x)^ T T7(x) 

- J ^-a 2 £jm^r7 t (x)£d r ^[t/ r (x)T7(x + OT)-t/_,.(x)T7(x-OT) 



— m 



2a ■ 



Y ri^(x-ar)d r U r (x-ar) - rf(x + ar)o r U- r (x + ar) — — T7(x) 



dx^a 2 ^^2^ *7 + ( x -ar)t/ r (x-af) - i] + (x + af)C/_ r (x + af) 



id+ 



U r (x)ri(x + ar) - U- r (x)ri(x - ar) 



(5.26) 



In the free limit, the above Hamiltonian becomes 



sd 



= J dx a 2 £|m 2 r7 t (x)^ T r](x) 

+^^[lH^ + at)-r 1 \x-ar)]^[r 1 (x + ar)-r 1 (x-ar)]^. (5.27) 

In the free field limit the two linear mass terms cancel with each other. 

When we implement the constraint equation on the lattice and use symmetric definition of 
the lattice derivative, it is important to keep in mind that we have only next to nearest neighbor 
interactions which can easily be seen from Eq. (5.27). Thus in each transverse direction even 
and odd lattice points are decoupled and as a result we have four independent sub-lattices in 
two dimensional transverse plane connecting (x 1 = even, x 2 = even), (x 1 = even, x 2 = odd), (x l 
= odd, x 2 = even) and (x 1 = odd, x 2 = odd) lattice points. 

Let us now address the nature of the spectrum and the presence and origin of doublers. 

5.3.2 Fermion doubling 

For clarity, the Hamiltonian (Eq. 5.27) can be rewritten as 

Kd = jdx-a 2 { 



£jm 2 T 7 t(x)^T 7 (x)-^ 2 [T 7 t (x)^£[T 7 (x + 2af)] + T 7 (x-2af)-2T 7 (x)] 
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+ £ [m 2 r 7 t (x)^r7(x)-^ 2 [r7 t (x)^ T £[r7(x + 2ar)] + r7(x-2ar)-2r 7 (x 

X e ,X ^ 

+ £ [ m 2^t (x) ^^ (x) _^ a 2^t (x) _i_^ [T?(x + 2af)] + 7?(x _ 2af) _ 2T?(x)] ] 

^mV^^x)-^ 2 ^)^ 



+ 



x ,x 



(5.28) 



where x\ stands for jc ! = even and x l () stands for x 1 = odd. Clearly the Hamiltonian is divided into 
four independent sub-lattices each with lattice constant 2a. As a result, a momentum component 
in each sub-lattice is bounded by ^ in magnitude. Again, going through the Fourier transform 
in each sub-lattice of the transverse space, we arrive at the free particle dispersion relation for 
the light front energy in each sector 

K = t+ t™ 2 + 4 E sin2 M • ( 5 - 29 ) 

ft CI y. 

For fixed k r , in the limit a — > , \ sin 2 k r a — > fc 2 and we get the continuum dispersion relation 



._ m 2 + k 2 
~k+ 



K = rz • (5-30) 



Because of the momentum bound of # doublers cannot arise from ka = it in sharp contrast 
with Euclidean lattice gauge theory where doubers come from ka = it [14]. However, because 
of the decoupling of odd and even lattices, one can get zero transverse momentum fermions one 
each from the four sub-lattices for a two dimensional transverse lattice. Thus we expect a four 
fold degeneracy of zero transverse momentum fermions. 

5.3.3 Numerical Investigation 

Using DLCQ for the longitudinal direction, Eq. (5.27) can be written as 

p- = k Hsd =L[ Hm + Hk ] (5.31) 

Jv Jv 



where 



H m = aV£££j 

l a z 1 

[b\l,z, a)b(l,z, o) + d\l,z, a)d(l,z, a)} (5.32) 
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and 



= EEEE) 



/ a z r ^ 

^(Z, z + of, o)b(l,z + ar, a) + z — ar, o)b(l, z — af , a) 

— b\l,z + ar, <j)b(l, z — af , a) — /?' (7,z — ar, o)b{l,z + ar, a) 
+ d^(l,z + at, a)d(l,z + ar, a) + d\l,z — ar, a)d(l,z — ar, a) 

— d/ (I, z + ar, a)d(l, z — ar, a) — (I, z — ax, a)d(l, z + ar, a) 



(5.33) 



For each transverse direction, we have 2n+l lattice points where the fermions are allowed 
to hop. Since even and odd lattice points are decoupled we need to fix the boundary conditions 
separately for even and odd sub-lattices in each transverse direction. 

Boundary conditions 

Fixed boundary condition: To implement the fixed boundary condition, we need to con- 
sider 2n + 5 lattice points. For one sub-lattice we have to fix particles at s = 1 and s = 2n + 5. 
We have, the wavefunction at location s, u s ~ sin (s — \ )ka which gives u s = for s = 1. We 
also need u s = for s = 2n + 5. Thus (2n + 4)k p a = pn, with p = 1,2,3, ...,«+ 1. For n=\, 
allowed values of k p are k p = ^, ||. 

For the other sub-lattice, we fix the particles at s = 2 and s = 2n + 4. The wavefunction at 
location s, u s ~ sin (s — 2)ka. u s = for s = 2 and s = 2n + 4. Thus (2n + 2)k p a = pit with 
p = 1,2,3, ....,n. For n = 1, only allowed value of k is k = ^. 

Combining the two sub-lattices, for n — 1, the allowed values of £ are ^, ^, and ||. 

Periodic boundary condition: To implement the periodic boundary condition, we need to 
consider 2n + 3 lattice points when fermions can hop in 2n + 1 lattice points. For one sub-lattice 
(2n + 3) rd lattice point is identified with the lattice point 1. For the other sub-lattice (2n + 2) nd 
lattice point is identified with the lattice point 2. Wavefunction at point s, u s ~ e lska . We require 
e ika = e i(2n+3)ka_ Thus k p a = ±j^^, p = 0,1 ,2, . For n = 1, we have, &o = 0, k\ = ±£. 

For the other sub-lattice we require e 2lka = e A 2n + 2 ) ka . Thus the allowed values of momentum 
k p are k p a = ±f p, p = 0, 1,2, For n=l, allowed value of k = 0. Thus for n = 1, taking 

the two sub-lattices together, the allowed values of k are 0, 0, ^. 
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Numerical results 

The results of matrix diagonalization in the case of the symmetric derivative with fixed 
boundary condition are presented in Figs. 5.4, 5.5 and 5.6. In Fig. 5.4 we present the lowest 




Figure 5.4: First four eigenvalues of Hamiltonian with symmetric derivative as a function of n. 

four eigenvalues as a function of n. At finite volume, the four states do not appear exactly 
degenerate even though the even-odd and odd-even states are always degenerate because of the 
hypercubic (square) symmetry in the transverse plane. The four states become degenerate in 
the infinite volume limit. The eigenfunctions of the lowest four states are presented in Fig. 5.5 
for n = 5 and fixed boundary condition. As they carry small nonzero transverse momenta they 
are not flat in the transverse plane. But, as they correspond to particle states, they are nodeless 
and same in shape and magnitude. All other states in the spectrum have one or more nodes. For 
example, in Fig. 5.6 we show the eigenf unction corresponding to the fifth eigenvalue for n = 5 
which clearly exhibits the node structure. 

With periodic boundary condition we can achieve exactly zero transverse momenta for any 
n and hence get four fold degenerate lowest eigenvalue corresponding to four zero transverse 
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Figure 5.5: Eigenfunctions of first four (degenerate) states for the case of fermion doubling with 
symmetric lattice derivative. 

momentum fermions. Corresponding wavefunctions are flat in transverse coordinate space and 
excited states (nonzero k states) show the expected node structures. 

5.4 Staggered fermion on the light-front transverse lattice 

As we have seen in the previous section that the method of symmetric derivatives results in 
fermion doublers, we now consider two approaches to remove the doublers. In this section we 
study an approach similar to the staggered fermions in conventional lattice gauge theory. In the 
next section we will take up the case of Wilson fermions. 
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Figure 5.6: Eigenfunction corresponding to the fifth state. 



In analogy with the Euclidean staggered formulation, define the spin diagonalization trans- 
formation 



We see from the QCD Hamiltonian given in Eq. (5.26) with symmetric derivative that in the 
interacting theory (except for the linear mass term) and also in the free fermion limit, even 
and odd lattice sites are decoupled and the Hamiltonian is already spin diagonal. So, it is very 
natural to try staggered fermion formulation on the light front transverse lattice. In this section 
we shall follow the Kogut-Susskind formulation [15] and present an elementary configuration 
space analysis for two flavor interpretation. After the spin transformation the linear mass term 
in the Hamiltonian (Eq. 5.26) becomes: 



rj(x h x 2 ) = (a 



(5.34) 




(5.35) 
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where, 0(x, r) = 1 for r = 1 and 0(x, r) = (— l) Xl for r = 2. After spin diagonalization, the full 
Hamiltonian in the free field limit becomes 

Psf = |rf*-a 2 £ jmV(x)^*(x) 

+^Y,[xH^ + at)-xH^-ar)]^ T [x(x + ar)-x(x- at)} 

- f (x) £ (x, r) \% (x + af ) - % (x - a?) ] 

-^m£[z t (x + af)-^( x _ fl r)]^(x,r)^ TZ (x)|. (5.36) 

The two linear mass terms cancel with each other in the free theory, but since they are present 
in the interacting theory we keep them to investigate the staggered fermions. 

Since all the terms in Eq. (5.36) are spin diagonal, we can put only a single component field 
at each transverse site. From now on, all the #'s and x^' s appearing in Eq. (5.36) can be taken 
as single component fermion fields. Thus we have thinned the fermionic degrees of freedom by 
half. Without loss of generality, we keep the helicity up component of % at each lattice point. 
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Figure 5.7: Staggered lattice. 
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Apart from the linear mass term in Eq. (5.36), all the other terms have the feature that 
fermion fields on the even and odd lattices do not mix. Let us denote (see Fig. 5.7) the even- 
even lattice points by 1, odd-odd lattice points by 1', odd-even lattice points by 2 and even-odd 
lattice points by 2', and the corresponding fields by %\, %v, X2 and %y. Then the first of the 
linear mass terms 

£z t (x)^ T £0(x,r)[ Z (x + ar)- Z (x-ar)] (5.37) 

x W r 

can be rewritten as (suppressing factors of a from now on), 

+ 1 +1 
Zi T ^(Vi#2 + V2#2') + #2 T ^(Vi£i-V 2 ;^) 

+Zi' t ^ T (ViZ 2 ' - V2Z2) +Z2' t ^ T (ViZi' + V2Z1) +B (5.38) 

where V] and V 2 are the symmetric derivatives in the respective directions and B represents 
the contribution from other blocks. Looking at Fig. 5.7 it is apparent that these Vi and V 2 can 
also be interpreted as a block derivative, i.e., finite differences between block variables. For 
example, V\Xi = Zl(l>0) — #i(0,0) where (1,0) and (0,0) are the block indices as shown in 
Fig. 5.7. 

Using Eq. (5.34), in terms of the nonvanishing components of 7], we have 

Vi=Xu 12 = 1X2, T]v=iXv, l2' = -X2'- (5.39) 

An interesting feature of lattice points 1 and V is that fermion fields r\\ and r\i< have positive 
helicity. r/2 and r\ 2 i have negative helicity. In terms of r\ fields the expression given in Eq. 
(5.38) can be written as 

n^-^xC-iVi^- V 2 T? 2 ') +/T7 2 t T^-(Vir]i +iV 2 TfiO 

+/T7i' t ^ T (-Vi77 2 / + /V 2 7] 2 )-T] 2 / t — (-/V lT j 1 , + V 2 i7i)+B. (5.40) 

Now, 

1?(1)-TJ(0) = I(i 7 (l)- n (-l)) + i(i 7 (l) + i 7 (-l)-2i 7 (0)) 

= Vl(0) + \v 2 ri(0), (5-41) 
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77(0)-r7(-l) = i(T7(l)-T7(-l))-i(77(l) + T7(-l)-277(0)) 
= VT7(0)-iv 2 r7(0) 



(5.42) 



where V and V z are respectively first order and second order block derivatives. So, we can write 
the expression (5.40) as 

1, 



■4-1 C \ — 1 — O "1 

tf+ I" /(Vl7?2 " 2 V ^ 2) ~~ (V27?2 ' " 2 V 2 T?2 7 
+/r 72 t ^ T {(V 1 r 7l + ^Vf7 ?1 ) + /(V 2 i7i' - ^Hi')} 

+ I ' 1 ii' t ^+ { ~ (V1172' + ^V?TJ2') + i(V 2 T? 2 + ^ife) } 

{ -i(VlT?l' - ^?7]l0 + (V 2 T?1 + l^hl)} ■ 



Let us introduce the fields 



Ml = -J=(l7l + 77l') 



M 2 

</ 2 



1 

7! 
1 

7! 
i 
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(772 + 772') 
(772-772')- 



(5.43) 



(5.44) 



Then, the first order derivative term in Eq. (5.43) can be written as 

w + w + w + 

where, d = & l d and the flavor isospin doublet 

/ = 



U 
d 



Similarly, we can write the second order block derivative term in expression (5.43) as 

2 id+ r 
where, T r s are the matrices in the flavor space defined as 

T l = -ia 2 , T 2 = -io l . 



(5.45) 



(5.46) 



(5.47) 



(5.48) 
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Similarly, the second term in Eq. (5.36) 

Y [tj + (x + at ) - ri f (x - at) ] -Jr- [tj (x + at) - tj (x - at)] (5.49) 

reads as 

*rf 1 j F *rf+ V^/t J_v2 /+ L { v r f_L r T r V 2 r f + V 2 f-L a r T r V r f] . (5.50) 

The full Hamiltonian given in Eq. (5.36) can now be written in two flavor notation (restoring 
factors of a) as 

P; f = Jdx~a^ { ™ 2 f^f+\lVrf^Vrf + a 2 nf^V 2 r f 

+ l ^rf^o r rv 2 j + v 2 r f^c r rv r f)\ 

-\m(f^a r V r f+^^rPjf+h.c)} . (5.51) 

The above simple exercise shows that applying the spin diagonalization on the symmetric 
derivative method, the number of doublers on the transverse lattice can be reduced from four to 
two which can be reinterpreted as two flavors. Although in the free case given by Eq. (5.51) 
the second and third lines are separately zero identically, we have kept these terms because in 
QCD similar terms will survive. These terms exhibit flavor mixing and also helicity flipping. 
The flavor mixing terms are always irrelevant. 

5.5 Wilson fermion on the light-front transverse lattice 

Since doublers in the light front transverse lattice arise from the decoupling of even and 
odd lattice sites, a term that will couple these sites will remove the zero momentum doublers. 
However, conventional doublers now may arise from the edges of the Brillouin zone. A second 
derivative term couples the even and odd lattice sites and also removes the conventional dou- 
blers. Thus, the term originally proposed by Wilson to remove the doublers arising from ka = % 
in the conventional lattice theory will do the job [8]. 

To remove doublers, add an irrelevant term to the Lagrangian density 

5jS?(x) = ^ E VK x ) t^( x ) V( x + a *) - 2 + ^M x ) V(* - at)] (5.52) 
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where K is the Wilson parameter. The constraint equation for i// in the presence of the Wilson 
term becomes 



id + Y (x) = mf\lf + (x) 

M Ya ¥+ ^ + ^ ~ U - r< ^ ¥+ ^ X ~ at ^ 

- - f [U r {x) y + (x + ar) -2\\r + (x) + C/_ r (x) y + (x - ar )] . (5.53) 
a 

The Wilson term generates the following additional terms in the Hamiltonian (5.26): 

P~ = -/^-a^j^l^x)^^^ 

k 1 r i 1 

-4— — 52 *7 (x-ar)6>t/ r (x-af) - r7 t (x + ar)d r i7_ r (x + ar) ^^r](x) 



a la 
+ J dx~a 2 ^ 



K 1 

a la 



£ £ \r\ f (x - af ) U r (x - af ) + r\ f (x + af ) C/_ r (x + ar) 



1 



K 1 



L r . v (x)r](x + as) -£/_ s (x)?7(x-as) 



— — ^52 T7 t (x-ar)d r t/ r (x-ar) - r] t (x + ar)d , r 6 r - r (x + ar)j 



[/$ (x) T] (x + as) + U- s (x) 77 (x - as) 



~y a'^a 2 £| A i^r7 t (x)^£[t/ r (x)r7(x + ar)+t/_ r (x)r7(x-ar) 



jc r 1 1 

— V r] t (x-af)i7 r (x-af) + r] t (x + af)L r _ r (x + af) — — t](x) 

- I ' dx~ a 2 Yr2Y*lL i]\x-ar)U r (x-ar) + ri\x + ar)U- r (x + ar) 

J x a r s 

1 



t/j (x) r\ (x + as) + C/_ s (x) 7] (x - as) 

In addition, the factor m 2 in the free term in Eq. (5.26) gets replaced by ji 2 = (m + 4-|) 2 
In the free limit the resulting Hamiltonian goes over to 



(5.54) 



/ a 2 52 

+^E[T7 + (x + «r)-r 7 t (x-af)] T i T l[r 7 (x + af)-r 7 (x-af)] 



?'<9+ 2a 
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K 2 1 

+^Y,[ r lHx + ar)-2ri jf (x) + ri" i (x-ar)] — [tj (x + of) - 2tj (x) + 77 (x - a?)] 



r 

-2^l,r 1 \x)^ T [r 1 (x + ar)-2r 1 (x)+r 1 (x-ar)] 



We rewrite the free Hamiltonian (5.55) as 



The diagonal terms are 

= |jx-fl 2 ^T]t(x)^ T T](x) 

The nearest neighbor interaction is 

P 0D1 = - / <frTa 2 ££ 

^ X f 



jii 2 + ^ + 8jiiK- + 12jr- y 



or 



The next to nearest neighbor interaction is 



OD2 



77 f (x + at) t^t t? (x - or) + 77 + (x - at) rf(x + at) 



Using the Fourier transform in the transverse space, we get 

1 



2 [ 

= I dx ~ ! "rl^^/F^ 



(2*)*) 



sinL-a 



' 1 k r a 



(5.55) 



(5.56) 



(5.57) 



(5.58) 



(5.59) 



(5.60) 



Note that, as anticipated, Wilson term removes the doublers. Since the Wilson term intro- 
duces nearest neighbor interactions, the sub-lattices are now coupled to each other and we have 
only one transverse lattice. The lowest eigenvalue in Eq. (5.60) occurs only if all the k r 's are 
zero and there are no more doublers in the theory. 



94 



5.5.1 Numerical Investigation 

For our numerical investigation, we write the Hamiltonian (Eq. 5.56) in DLCQ as 



H w — Hd + H d\ + Hqd2 



(5.61) 



with 



H D = 



[a V + 1 + 8«AiK + 12k 2 ] £££ - 

1 (J z ' 

[b\l,z,a)b(l,z,o) +J + (Z,z, o)d(l,z, a)], 



H 0Dl = -[2k^ + 4k 2 ]££££) 



/ a z r I 

b\l,z, o)b(l,z + ar, a) +l? t (/,z, o)b(l,z — ar, a) 
+ d\l,z, a)d(l,z + ar, a) +d jf (l,z,a)d(l,z — ar, a) 



and 



Horn = -Ij-^EEEE) 



/ a z r 

b\l,z + ar, c)b{l, z — ar, a) + b\l,z — ar, a)b(l,z + ax, a) 
+ d^ (I, z + ar, a)d(l,z — ar, a) + d\l, z — ar, o)d (I, z + ar, a) 



Boundary condition 



(5.62) 



(5.63) 



(5.64) 



With the Wilson term added, we do not have decoupled sub-lattices. We have both nearest 
neighbor and next- to-nearest neighbor interactions. Since with fixed boundary condition, the 
lowest four eigenvalues are not exactly degenerate in finite volume, it is difficult to investigate 
the removal of degeneracy by the addition of Wilson term. With periodic boundary condition, 
for a lattice with 2n + \ lattice points in each transverse direction, we identify the (2n + 2) th 
lattice site with the first lattice site. Then for the Hamiltonian matrix we get the following 
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additional contributions. 



H = 



NN N \ 
NN 



NN 
\N NN 



(5.65) 



/ 



The matrix elements AW = — \ + K 2 and N = —lajiK — 4k 2 . For a given n, the allowed values 

of k are £ p a = ±J^y, p = 0, 1 , 2, Thus for n = 3, we expect multiples of ^ apart from 0. 

For n = 5, apart from 0, allowed values of k are multiples of j{. 

Numerical results 

Since the Wilson term connects even and odd lattices, the extra fermions that appear at zero 
transverse momentum are removed once Wilson term is added as we now have nearest and 
next to nearest neighbor interactions. For large n, we get the expected spectra but, numerical 
results suggest that the finite volume effect is larger for small K which is obvious because K is a 
mass-like parameter. For small K the wavefunctions become fat and requires larger lattice to fit 
into. For example, with periodic boundary condition, for n = 3, for K = 1 .0, 0.5, 0.4, we get the 
expected harmonics but not for K = 0. 1 . The situation is similar for n = 5. For n = 10, expected 
harmonics emerge even for k = 0.1 but not for K = 0.01. Since all realistic calculations are 
done in a small lattice, it is desirable to have the Wilson parameter K not too small. 

5.6 Doubling and symmetries on the light front transverse lattice 

Let us now try to understand the fermion doubling in terms of the symmetries of the trans- 
verse lattice Hamiltonians. We are aware that there are rigorous theorems and anomaly argu- 
ments in the conventional lattice gauge theories [16] regarding presence of fermion doublers. 
In standard lattice gauge theory, some chiral symmetry needs to be broken in the kinetic part 
of the action to avoid the doublers. On the light-front, chirality means helicity. For example, 
a standard Wilson term which is not invariant under chiral transformations in the conventional 
lattice gauge theory, is chirally invariant on the light-front in the free field limit. The question 
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is then why the Wilson term removes the doublers on the light-front transverse lattice. The 
argument that there is nonlocality in the longitudinal direction cannot hold because, in the first 
place, having nonlocality is not a guarantee for removing doublers and secondly, there is no 
nonlocality on the transverse lattice. One, therefore needs to find a reasoning that involves the 
helicity in some way. 

Because of the constraint equation which is inconsistent with the equal time chiral trans- 
formation in the presence of massive fermions, we should distinguish between chiral symmetry 
in the equal time formalism and in the light-front formalism. For example, the free massive 
light-front Lagrangian involving only the dynamical degrees of freedom is invariant under y 5 
transformation. On the light-front, helicity takes over the notion of chirality even in presence of 
fermion mass which can be understood in the following way. 

In the two component representation [12] in the light-front formalism, let us look at the 
objects and y/jj". We have 



with 



T){X) 



7]{x) 




m(*) 



(5.66) 



(5.67) 



The projection operators are P R = 1(1 + y 5 ) and P L = Ul — y 5 ) with 



2 V ' 1 I ) 1 L — 2 < 

3 



r 5 



o 



(5.68) 



Then 



Yr=PrY + = 



+ _ 



and 



o 



V°/ 
/o\ 





(5.69) 



(5.70) 



Thus = PR\y + represents a positive helicity fermion and = PlW + represents a negative 
helicity fermion, even when the fermion is massive. This makes sense since chirality is helicity 
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even for a massive fermion in front form. This is again to be contrasted with the instant form. 
In that case the right handed and left handed fields defined by \j/ R = P R \j/ = ^(1 + y 5 ) V and 

= P^Xj/ = i(l — y 5 ) \]/ contain both positive helicity and negative helicity states. Only in the 
massless limit or in the infinite momentum limit, ty/^ becomes the positive helicity state and ty/^ 
becomes the negative helicity state. 

As a passing remark, we would like to mention that in continuum light front QCD there is a 
linear mass term that allows for helicity flip interaction. 

In lattice gauge theory in the Euclidean or equal time formalism, because of reasons con- 
nected to anomalies (the standard ABJ anomaly in vector-like gauge theories), there has to be 
explicit chiral symmetry breaking in the kinetic part of the action or Hamiltonian. Translated 
to the light front transverse lattice formalism, this would then require helicity flip in the kinetic 
part. A careful observation of all the above methods that get rid of fermion doublers on the light 
front transverse lattice reveals that this is indeed true. 

In particular, we draw attention to the even-odd helicity flip transformation 

ri(x h x 2 ) - {d^{d 2 ) x ^{x u x 2 ) (5.71) 

that was used in Sec. 5.4 for spin diagonalization. It should also be clear that the form of the 
above transformation is not unique in the sense that one could exchange C\ and c 2 and their 
exponents x\ and x 2 could be changed by ±1 . 

Note that the Hamiltonians Pj b given in Eq. (5.11) and P~ given in Eq. (5.55) that do not 
exhibit fermion doubling are not invariant under the transformation Eq. (5.71). On the other 
hand the Hamiltonian P~ d given by Eq. (5.27) that exhibits fermion doubling is invariant under 
this transformation. 

5.7 Summary 

The presence of the constraint equation for fermions on the light front gives rise to interest- 
ing possibilities of formulating fermions on a transverse lattice. We have discussed in detail the 
transverse lattice Hamiltonians resulting from different approaches in two different boundary 
conditions. 
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In the first approach, we have proposed to use forward and backward lattice derivatives 
respectively for y + and (or vice versa) so that the resulting Hamiltonian is Hermitian. 
There is no fermion doubling. The helicity flip (chiral symmetry breaking) term proportional to 
the fermion mass in the full light front QCD becomes an irrelevant term in the free field limit. 
With periodic boundary condition one can get the helicity up and helicity down fermions to be 
degenerate for any transverse lattice size n. With fixed boundary condition, there is a splitting 
between the two states at any n but the splitting vanishes in the large volume limit. 

In the second approach, symmetric derivatives are used for both and . This results 
in four fermion species. This is a consequence of the fact that the resulting free Hamiltonian 
has only next to nearest neighbor interactions and as a result even and odd lattice sites get de- 
coupled. One way to remove doublers is to reinterpret them as flavors using staggered fermion 
formulation on the light front. In QCD Hamiltonian, it generates irrelevant flavor mixing in- 
teractions. However, in the free field limit, there is no flavor mixing. Another way to remove 
the doublers is to add a Wilson term which generates many extra terms in the Hamiltonian. In 
the free field limit, only the helicity nonflip terms survive. The Wilson term couples even and 
odd sites and removes the doublers. Numerically, we found that in small lattice volumes it is 
preferable to have not too small values of the Wilson mass k/ci. 

Chiral symmetry in light-front is different from the chiral symmetry in equal time formal- 
ism. So, fermion doubling on LFTL should be related with the chiral symmetry on the light- 
front. Light- front chirality is equivalent to helicity even for massive fermions. We have identi- 
fied an even-odd helicity flip symmetry of the light front transverse lattice Hamiltonian, absence 
of which means removal of doublers in all the cases we have studied. 
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CHAPTER 6 



Meson Bound States in Transverse Lattice QCD 
6.1 Introduction 

In the previous chapter we have discussed different ways of formulating fermions on a 
light-front transverse lattice. We have also discussed the advantages and disadvantages of all 
the methods in free field limit. But, the true testing ground for the strengths and weaknesses of 
different methods is the full QCD. To complete this comparative study, in this chapter we con- 
sider the meson bound state problem in (3+1) dimensional light-front QCD with two transverse 
directions (x 1 ,* 2 ) discretized on a square lattice. 

In this chapter we make a detailed comparison of two different light-front QCD Hamilto- 
nians. One is the Hamiltonian with fermions formulated using forward and backward lattice 
derivatives and the other Hamiltonian with fermions formulated with symmetric lattice deriva- 
tive with Wilson term [1] to get rid of doublers on the transverse lattice. Light-front staggered 
lattice formalism [2, 3] to remove doublers in case of symmetric lattice derivative that we have 
discussed in the previous Chapter is a different game altogether and will not be further inves- 
tigated in this work. For our calculation, we adopt the one link approximation in the meson 
sector which has been widely used in the literature [4, 5, 6]. (Only very recently, the effect of 
additional links in the meson sector has been investigated [7]). One link approximation is too 
crude to reproduce physical observables. So, rather than fitting the parameters to reproduce any 
physical observable we concentrate here to investigate the effects of various coupling strengths 
on the low-lying spectra and wave functions and compare two different formulations. 

We use Discretized Light Cone Quantization (DLCQ) [8] to address longitudinal dynamics. 
Because of the presence of severe light-front infrared divergences, a major concern here is 
the reliability of DLCQ results when calculations are done at finite resolution K and results are 
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extrapolated to the continuum (K — > oo). In meson calculations so far, K < 20 have been chosen. 
We perform a detailed study of the continuum limit of DLCQ by performing calculations at 
larger values of K [9]. 

In the meson sector, in the zero link approximation, at each transverse location we have a 
two-dimensional field theory which in the large N c limit (where N c is the number of colors) is 
nothing but the 't Hooft model. In this well-studied model, excited states are simply excitations 
of the qq pair, which contain nodes in the wavefunctions. The picture changes when one link 
is included thereby allowing fermions to hop. The admixture of qq link states with qq states is 
controlled by the strengths of the particle number changing interactions and the mass of the link 
field. One link approximation is a priori justified for very massive links and/or weak particle 
changing interaction since in this case low lying excited states are also qq excitations. Likewise, 
for large particle changing interaction strength and/or light link mass, low lying excited states 
are qq link states. We explore the spectra and wavefunctions resulting from the choice of various 
regions of parameter space. 

Details of the derivation of the fermionic part of the Hamiltonian are already discussed in 
the previous chapter. Here we give the details of the gauge field part of the QCD Hamilto- 
nian. Non-linear constraints on the unitary link variables make it difficult to perform canonical 
quantization. We also present the effective Hamiltonian when non-linear unitary variables are 
replaced by linear variables. 

6.2 Gauge field part of the Lagrangian density 

The gauge field part of the Lagrangian density in the continuum is 

^G = ^TrF pa FP a (6.1) 

where F po = dP A a -d a AP + [A p ,A a ] with A p = igAP a T a . Here p, a = 0, 1,2,3 and a = 
1, 2, . . . , 8. For ease of notation we suppress the dependence of field variables on the longitu- 
dinal coordinate in this section. With the gauge choice A + = 0, the Lagrangian density can be 
separated into three parts, 

J? G = ^T + ^L + ^LT. (6.2) 
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Here ££ T depends entirely on the lattice gauge field U r (x). 



-TjYX Tr\u r (x)U s (x + ar)U- r (x + ar + as)U- s (x + as)-l] I, (6.3) 

8 a r^s { J 



r, s = 1,2. The purely longitudinal part Jzf^ depends on the constrained gauge field A , 

J? L = l -{d + A- a f (6.4) 
and the mixed part J£ LT depends both on lattice gauge field and the constrained gauge field. 
^lt = ^Tr[F^ = ^Tr[D^)(D^U r (^] 



1 



g L a L 



Tr 



d^U r (x)d^Uj(x) +A M (x)[C/ r (x) d» t/t (x) ] 



+A M (x + ar)[t/ r t (x) <9 M U r (x) 



(6.5) 



2n = ^Tr[d^U r (x)d^Uj(x)] + gA - a J+ a NK . (6.7) 



where fj, = +, — only and the mixed co variant derivatives of the link variables are defined by 
(see Appendix H for the derivation) 

D M C/ r (x) = dpU r (x) +A^U r (x) -C/ r (x)A M (x + af). (6.6) 

In the A + = gauge, we can write J^lt as 

l^rr[^C/ r (x)^C/;(x)] + -i 

Here the link current 

Jli%k(*) = £^rjr"[[/ r (x)/i+ C/;(x) + C/;(x-af)/i+ C/ r (x-«f)]|. (6.8) 

Substituting back the expression for A a from the constraint equation 

(<?+) 2 A-« = J (J+% K -/+«)) (6.9) 

with 

7+ a (x) = 277 t (x)r a r7(x) (6.10) 

where 7] is the dimensionless two-component lattice fermion field, in the A~ a dependent terms 
in the Lagrangian density, namely, 

-^A-«/+« + I(a + A-«) 2 +I|A-«^ (6.11) 
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we generate the terms 

+1 / 1 \ 2 



ytS*{pl ] u^b^{lF) ^r'n- s -A{^) n<r«n. (6.12, 

Collecting all the terms, the canonical Lagrangian density for transverse lattice QCD is 



JSf = J? f + -^Tr[d^U r (x)d^(x)] 

a, g 

1 



7> U r (x)U s (x + ar)U- r (x + ar + as)U- s (x + as)-l 

a s \% 1 



^Z„a J LINK [ 2+ ) J L1NK + r ) „A8 Z 'Jq U [ 3 + ) ^ 



1 



1 



2 T +a 



2a* 6 q 



8 T+a 



a 



4 LINK 



1 



T+a 



(6.13) 



where =2/ is the fermionic part of the QCD Lagrangian density with forward and backward 
lattice derivatives given by 



f 



a 



1 ri\x)id r ](x)-^r7 t (x)^ T r7(x) 



a 



1 



+?'m^r7 t (x)d v -4r U s (x)r](x + as) - r](x) 



1 r 

+im—~ 

a L 



1 1 



j] f (x + at ) U} (x) - j] t (x) j a r - — 7] (x) 



^[T ? t (x + af)^(x)-r 7 t (x)]d r ^ T d v [t/. v (x)r 7 (x + as)-r 7 (x)], (6.14) 



or with symmetric lattice derivative with Wilson term given by 
&f = \riHx)id-ri(x)-\(m + 4-) 2 ri^x)^-ri(x) 

+^(^ + 4^)^|77 t (x)£a r ^ T [t/,(x)r 7 (x + af)-t/_ r (x)r 7 (x- 

E 77 t (x-af)d r L r r (x-ar) - r] t (x + flf)a , r 6 r - r (x + OT) -^17 (x) 



ar 



1 ?c 1 



— y< -^EE T7 t (x-ar)[/ r (x-OT) + r] t (x + ar)t/_ r (x + OT) 
a L \ a2a'—"—' 



1 



?c 1 

a 2a 



U s (x)ri(x + as)-U- s (x)ri(x-as) 



EE r 7 t ( x -OT)d r t/ r (x-ar) - (x + ar)d r U- r (x + ar) 

r s 
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1 



id+ 
1 1 



a 2 4a 2 
1 



U s (x)ri(x + as) + U- S (x)j](x - as) 



£ £ \r) 1 (x - at) a r U r (x - ar ) - r\ f (x + at ) 6 r U r (x + at) 



id+ 



U s (x)vj (x + as) — U— s (x)t] (x — as) 



+ J_ ( m + 4 -) - J T]t(x)^ T £[C/ r (x)T](x + flf)+C/_ r (x)T](x-flf)] 

+ T7 t (x-ar)6 r r (x-ar) +r] t (x + ar)L r - r (x + ar) t^t](x) 



1 K 2 



a <ar 



H — 2~im T7 + (x-af)/y r (x-af) + T7 t (x + af)/7_ r (x + af) 



C/j(x)7](x + fls) + £/_ iS -(x)T7(x-as) 



(6.15) 



Here we use the two-component representation [10] for the dynamical fermion field 







(6.16) 



where r\ is the dimensionless two component lattice fermion field and &\ = (72, <h = —Oi. 
Linearization of the link fields 

Because of the nonlinear constraints U^U = 1, det U = 1, it is highly nontrivial to quantize 
the system. Hence Bardeen and Pearson [11] and Bardeen, Pearson, and Rabinovici [12] pro- 
posed to replace the nonlinear variables U by linear variables M where M belongs to GL(N 7 c €), 
i.e., we replace ^U r (x) — > M r (x). This linearized approximation is somewhat meaningful only 
on a coarse lattice. Justification may come from the fact that particle structure of a hadron is a 
long-distance property which can be described by some effective variables without any explicit 
details of the microscopic variables [12]. The x~ independent residual gauge invariance on the 
transverse plane is still preserved with the linear gauge fields (see Appendix I for details). Once 
we replace U by M, many more terms are allowed in the Hamiltonian. Thus one needs to add 
an effective potential V e /f to the Lagrangian density 



V eff = -IL Tr(M^M) + Ai Tr[(M^M) 2 ] + X 2 [detM + H.c] + .... 



(6.17) 
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6.3 Effective Hamiltonian 



6.3.1 Hamiltonian with forward and backward derivatives 



Once we replace the nonlinear link fields by linear link field by including the effective poten- 
tial, we perform canonical quantization and construct the effective Hamiltonian for transverse 
lattice QCD. 

Thus, the effective Hamiltonian for QCD on the transverse lattice, when fermions are put in 
with forward and backward lattice derivatives, becomes 



fb ~ l ffree 



P f free + P V + P fhf + P hf + P chnf 



The free fermion part is 



/ free 



The effective potential part is 



+P +P +P +P 

1 qqc 1 ggc 1 1 qgc 1 1 p 



/2 1 
dx~ £ (m 2 + — ) 7] f (x) 7] (x) . 



(6.18) 



(6.19) 



Jdx-a 2 Y,(tL 7>(M f M)-Ai Tr^M) 1 }- X 2 [det M + H.c] + ... 



The free helicity-flip part is 

/I \ 
x s a ° 

Helicity flip associated with the fermion hop is 

f t 1 1 r 

P hf = ~ im 8 / ^"EEn T ( x )*j-^+ M s (x)J](x + as) 
J x s ad 

/r ill 
dx'YJL T? f (x + af)M r + (x) & r -—ri(x). 

Canonical helicity non-flip terms are 

p ; hnf = 4/^ 2 EDV(x+rfR t (x)]4^['iW] 

a j x rs to 

"4 / ^^ 2 H[r7 t (x)]^^ T d,[M,(x)T 7 (x + as)] 
a j x rs la 

-4 / ^« 2 ££[r7 t (x + af)M r t (x)]d r ^ T d,[M,(x)r 7 (x + as)]. 
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(6.20) 



(6.21) 



(6.22) 



(6.23) 



The four-fermion instantaneous term is 

P~ c = -2^|^-£7 7 t(x)rr 7 (x)^r 7 t (x)rr 7 (x). (6.24) 
The four link instantaneous term is 

Pggc = ~2a^ I ^ X E^£/wa:( x ) (d+yJ~UNK.( x )- (6.25) 



The fermion - link instantaneous term is 

P «gc = 2 ^J dx ~L J UNK(*)j^2ilHx)T a Tl(x). (6.26) 
The plaquette term is 

P~ = [dx-a 2 ^Y, | Tr\M r (x)M s (x + ar)M_ r (x + ar + as)M_ s (x + as)-l]\. 

■a J x r ^ s [ J 

(6.27) 

Here 

J+% K (x) =Y t Tr^T a [M r {x)i d+ Mj(x) + M}{x-ar)i d+ M r (x-ar)} J. (6.28) 

6.3.2 Hamiltonian with the Wilson term 

When one uses symmetric derivatives for the fermion fields, doublers arise as a result of the 
decoupling of even and odd lattice sites. Here we use the Wilson term to remove the doublers. 
In this subsection, the details of the structure of the Hamiltonian resulting with the modification 
of the Wilson term are presented. 

The effective Hamiltonian for this case can be written as 

P = P f free+ P V + P hf + P whf 
JrP chnf + P wnf\ + P wnf2 

+P qqc + P ggc + P qg c + P P ■ (6-29) 

The free fermion part is 

Pi free = / ^E? (m + 4^) V (x) ^ (x) . (6.30) 
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The helicity flip part is 



P hf = -gjdx-^l ( m + 4 ^)lr 7 t (x)£d^[M r (x)r 7 (x + ar)-M_ ; -(x)r 7 (x-ar)] 

X I T 



f K*\ 1 r i i 

fm + 4— j— J2 ^ t (x-ar)a , r M r (x-ar) -r] + (x + ar)a , r M_ r (x + ar) t^T7(x 



(6.31) 



The Wilson term induced helicity flip part 

P~ hf = /|^-£|^££[r 7 t (x-af)M r (x-ar)+77 t (x + ar)M_ r (x + ar) 



1 



M s (x)r](x + as) -M- s (x)ri(x-as) 



id+ 
k 1 r 

— 2~EE ^ (x-flf)d>M r (x-af) - r7 t (x + OT)d r M_ r (x + OT) 



Af s (x) r] (x + as) + M_ s (x) 7] (x - as) 



(6.32) 



The canonical helicity non-flip term arising from fermion constraint is 
P~ hnf = V/^~E^2LE[^ + (x-a?)a r M r (x-af)-77t(x + af)d r M_ r (x + af) 



1 



Mj(x)T](x + as) -M_,(x)T](x-as) 



The Wilson term induced helicity non flip terms are 



(6.33) 



P~ nfl = -^|^-£|( m + 4 ^)^r 7 t (x)^£[M r (x)r 7 (x + ar)+M_ r (x)r ? (x-af)] 

/ k* \ /c r ~i 1 ! 

+ (m + 4-J-£ Tj t (x-af)M r (x-a?) + Ti t (x + af)M_ r (x + af) ^77(x) >. 



(6.34) 



and 



/k^ r 
£ — £ £ *1 f ( x ~ at ) M r (x - ar) + r] + (x + ar)M_ r (x + af ) 

M s .(x)r](x + as) +M_. s .(x)r](x-as) 



1 



id+ 



(6.35) 



Comparing the Hamiltonians with a) forward-backward derivative and b) symmetric deriva- 
tive with the Wilson term we notice that the only differences are in the particle number changing 
interactions, namely, helicity flip and helicity non-flip terms. 
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6.4 Meson bound state in one link approximation 

6.4.1 Relevant interactions 

In one link approximation, for either Hamiltonian, the four link instantaneous term and the 
plaquette term do not contribute and only the link mass term of the effective potential con- 
tributes. Further, in the case of the forward-backward Hamiltonian, the helicity non-flip part 
proportional to g 2 does not contribute. For the Wilson term modified Hamiltonian, the Wil- 
son term induced helicity flip part P~ h p the canonical helicity non-flip term P~ nh j- and the term 
proportional to K 2 in the Wilson term induced helicity non-flip part do not contribute. Thus 
in the case of the Wilson term modified Hamiltonian the entire fermion hopping with no he- 
licity flip arises from the Wilson term. In the case of forward and backward derivatives, terms 
are also present in one link approximation which violate hypercubic symmetry on the trans- 
verse lattice (see Appendix F). They become irrelevant in the continuum limit when the linear 
variables M are replaced by non-linear variables U. We have removed them entirely from the 
Hamiltonian in the present investigation. The Hamiltonian matrix elements in DLCQ for both 
forward-backward and symmetric derivative with Wilson term are explicitly given in Appendix 
J. 

6.4.2 Comparison with one gluon exchange in the continuum 

It is interesting to compare the one link approximation on the transverse lattice with the one 
gluon exchange approximation in the continuum. In the latter, a major source of singularity 
is the pr term in the quark - gluon vertex where k 1 - (k + ) is the gluon transverse (longitudinal) 
momentum. This originates from the A 7+ interaction term in the Hamiltonian via j+d^ -A 1 - 
contribution to the constrained field A . This term gives rise to quadratic ultraviolet divergence 
in the transverse plane accompanied by linear divergence in the longitudinal direction in fermion 
self energy. On the transverse lattice, d + A~ oc ^ T j+ NK so that A~/+ — > J^ INK ^+^ Jq '• Thus a 
term which gives rise to severe divergence structure in the continuum gets buried in the fermion- 
link instantaneous interaction term which gives rise to a term in the gauge boson fermion vertex 
in the continuum in Abelian theory. In the non-Abelian gauge theory this gives rise to a term in 
the quark-gluon vertex and also to the instantaneous quark-gluon interaction in the continuum. 
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The transfer of the troublesome term from quark-gluon vertex in the continuum theory to 
quark - link instantaneous interaction term in the lattice theory has an interesting consequence. 
In the continuum theory, the addition of a gluon mass term by hand spoils the cancellation of 
the light front singularity between one gluon exchange and the instantaneous four - fermion 
interaction. On the transverse lattice, this cancellation is absent anyway with or without a link 
mass term. 

6.4.3 Longitudinal dynamics and effects of transverse hopping 

We first consider the dynamics in the absence of any link. In this case, fermions cannot hop, 
and at each transverse location we have (1+1) dimensional light front QCD which reduces to 
the 't Hooft model in the large N c limit. In this case quark and antiquark at the same transverse 
position interact via the spin independent instantaneous interaction which, in the non-relativistic 
limit reduces to the linear potential in the longitudinal direction. The only parameters in the 
theory are the dimensionless fermion mass nif = am and the gauge coupling g. The spectrum 
consists of a ground state and a tower of excited states corresponding to the excitations of the 
qq pair. 

Next consider the inclusion of the qq link states. There are four independent amplitudes 
corresponding to whether the quark is on the left or right of the antiquark or, above or below the 
antiquark. With non-zero mass of the link, these states lie above the ground state of pure quark 
- antiquark system. Further the q, q and link (which are frozen at their transverse positions) 
undergo fermion - link instantaneous interactions in the longitudinal direction which further 
increases the mass of qq link states. Now the quark or antiquark can hop via helicity flip or 
helicity non-flip. Here we find a major difference between the Hamiltonians resulting from 
forward-backward derivative and symmetric derivative. Let us first consider the helicity flip 
hopping term in the forward-backward case 

P~ = -imgJdx-Y,Y< lH*)6 r ^^ri(x + ar)+ri\x + ar)Ml(x)6 r l--^ T ri(x) . 

(6.36) 

If we consider transition from two particle to three particle state by a quark hop, then the first 
term in Eq. (6.36) corresponds to | 2) — >| 3a) and the second term corresponds to | 2) — >| 3b). 
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The helicity flip term in symmetric derivative case, after making some shifts in lattice points, 
can be written as 

K\ 1 



77 f (x) c r j^M r (x) T) (x + at) - 7] f (x) c r M r (x) 7] (x + at) 

Gr-^+M-rix) 7] (x - at) - 7] f (x) c r M r (x) ^ 7] (x - at) 



(6.37) 



For the Hamiltonian with symmetric derivative, a quark or antiquark hopping accompanied 
by helicity flip has opposite signs for forward and backward hops. On the other hand, hopping 
accompanied by helicity non-flip have the same signs. As a result, there is no interference 
between helicity flip and helicity non-flip interactions [1]. In the case of the Hamiltonian with 
forward-backward derivative, quark or antiquark hopping accompanied by helicity flip has the 
same sign for forward and backward hops. As a consequence the helicity non-flip hop can 
interfere with the helicity flip hop. This has immediate consequences for the spectrum. In the 
case with symmetric derivative, in lowest order perturbation theory, the helicity zero states mix 
with each other which causes a splitting in their eigenvalues resulting in the singlet state lower 
than the triplet state. On the other hand, helicity plus or minus one states do not mix with each 
other or with helicity zero states resulting in a two fold degeneracy. In the case with forward and 
backward derivatives all helicity states mix with each other resulting in the complete absence 
of degeneracy. Obviously, one has the freedom to tune the free parameters to minimize the 
splitting. 

6.4.4 Singularities, divergence and counterterms 

Since the transverse lattice serves as an ultraviolet regulator, we need to worry about only 
light front longitudinal momentum singularities. 

Tree level 

We take all the terms in the Hamiltonian to be normal ordered. At tree level this leaves 
us with singular factors of the form t-U in the normal ordered four fermion and fermion link 
instantaneous interactions. The singularities are removed by adding the counterterms used in 

111 



the previous work [4] on transverse lattice. The explicit forms of the counterterms are given in 
Appendix J. 3 in the appropriate places. 

Self energy corrections 

In the one link approximation, a quark can make a forward (backward) hop followed by a 
backward (forward) hop resulting in self energy corrections. In a single hop, helicity flip or non- 
flip can occur. In the case of symmetric derivatives, helicity flip cannot interfere with helicity 
non-flip, and as a consequence, self energy corrections are diagonal in helicity space. In the 
case of forward and backward derivatives, the interference is nonzero resulting in self energy 
corrections, both diagonal and off-diagonal in the helicity space. Similar self energy corrections 
are generated for an antiquark also. These self energy corrections contain a logarithmic light 
front infrared divergence which must be removed by counterterms. In Appendix J. 5 we present 
the explicit form of counterterms in the two cases separately. In previous works on one link 
approximation [1, 4, 5], these counterterms were not implemented. For low K values one may 
not feel the divergence, but as one increases K the need of self energy counterterms are readily 
felt. 

6.5 Numerical Results 

We diagonalize the dimensionless matrix a 2 P~. We further divide the matrix elements by 
g 2 Cf which is the strength of the matrix elements for four fermion and fermion - link instanta- 
neous interactions. Now, define the constant G with dimension of mass by G 2 = %C f . DLCQ 
yields the eigenvalue M 2 = 

The dimensionless couplings are introduced [4] as follows. Fermion mass nif = m/G, 
link mass \x b = jl/G, particle number conserving helicity flip coupling mf/{aG) = m/Ci, 
particle number non-conserving helicity flip y/Ngmf/ (aG) = rtifCi, and particle number non- 
conserving helicity non- flip VNg/{a 2 G 2 ) = C 3 . In the case of the Wilson term modified Hamil- 
tonian, we have fermion mass term my = (m + 4?c/a) / G, helicity-flip coupling \/Ngmf/ (2aG) = 
m,fC2, and helicity non-flip coupling y/NgnifK/ (aG) = nifC^. 
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Figure 6.1: Effect of counterterm on the ground state eigenvalue, (a) With and without the 
counterterm in the qq sector for rtif = 0.3. (b) With and without the counterterm in the qq link 
sector for nif = 0.3 and \L\, = 0.2. 



All the results presented here were obtained on a small cluster of computers using the Many 
Fermion Dynamics (MFD) code [13] implementing the Lanczos diagonalization method in par- 
allel environment. For low K values, the results were checked against an independent code 
running on a single processor. 

Cancellation of divergences 

As we already mentioned, we encounter jp-^ singularities with instantaneous four fermion 
and instantaneous fermion - link interactions which give rise to linear divergences. We remove 
the divergences by adding appropriately chosen counterterms. We have numerically checked 
the removal of linear divergence by counterterms in DLCQ. First we consider only qq states 
with instantaneous interaction. We study the ground state eigenvalue as a function of K with 
and without the counterterm. Results are presented in Fig. 6.1 (a). Next we consider only qq 
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link states with fermion-link instantaneous interaction with and without the counterterms. The 
behavior of ground state eigenvalue as a function of K is presented in Fig. 6. 1 (b). In both cases, 
it is evident that the counterterms are efficient in removing the divergence. 

6.5.1 qq at the same transverse location 

Next we study the spectrum of the Hamiltonian in the absence of any links. Since, in this 
case, the Hamiltonian depends only on the dimensionless ratio ^ we fix g = 1 and vary nif 
to study the spectra. The Hamiltonian matrix is diagonalized for various values of K. The 
convergence of the ground state eigenvalue as a function of K is presented in Table 6.1. 



K 


Eigenvalue (J- 


r 1 ) 


nif — 0.3 


ntf — 0.9 


nif = 3.0 


10 


0.620 


4.547 


39.233 


18 


0.693 


4.664 


39.861 


30 


0.745 


4.724 


40.053 


50 


0.788 


4.762 


40.163 


78 


0.819 


4.783 


40.220 


98 


0.832 


4.791 


40.241 




0.869 


4.820 


40.285 



Table 6. 1 : Ground state eigenvalue (in units of G 2 ) for qq sitting at the same transverse location. 



The ground state wavefunction squared as a function of the longitudinal momentum fraction 
x is plotted in Fig. 6.2. The convergence of the wavefunction has a very different behavior as 
a function of fermion mass nif. As can be seen from this figure, the convergence in K is from 
above for heavy nif and from below for light nif. As a consequence the wavefunction is almost 
independent of K when nif is of order g. 
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Figure 6.2: Quark distribution function | \j/(x) | 2 of the ground state in the qq approximation for 
three choices of quark masses with coupling constant g = 1.0. 



6.5.2 Results of the one link approximation 

We encountered logarithmic infrared divergences due to self energy corrections and, in Ap- 
pendix J. 5, we discuss the associated counterterms. In Fig. 6.3 we show the effect of self energy 
counterterms on the ground state energy in the two Hamiltonian cases we studied. 

The convergence of lowest four eigenvalues with K for the Hamiltonian with forward- 
backward and symmetric lattice derivatives is shown in Table 6.2 for nif = 0.3, = 0.2. We 
also show the results extrapolated to K — > oo. The convergence of the eigenvalues in K is very 
slow and one really needs to go for large K. 

The quark distribution function for the ground state and the fifth state for the set of param- 
eters ntf = 0.3, jl b = 0.2, C 2 = 0.4, C 3 = 0.01 and K = 30 is presented in Fig. 6.4. In this 
figure we also present separately the contribution from two particle and three particle states. As 
expected, the contribution from the three particle state peaks at smaller x compared to the two 
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Figure 6.3: Effect of self energy counterterms on the ground state eigenvalue in the case of 
(a) symmetric derivative with £2 = 0.4, C3 = 0.1 and (b) forward-backward derivative with 
C 2 = 0.4, C 3 = 0.01. m f = 0.3, \x h = 0.2 for both cases. 

particle state. The exact location of this peak depends on the link mass. The ground state is 
dominated by two particle sector while the fifth state is dominated by three particle sector. We 
have shown only one of the first four states since they are similar looking states with different 
spin contents. Fig. 6.4 warns us that one link approximation is inadequate to study the excited 
states. 

It is interesting to see the effect of fermion - link instantaneous interaction on the low lying 
eigenvalues. In its absence, there is no confining interaction in the longitudinal direction in 
the qq link sector. Furthermore, the mass of the lowest state in this sector corresponds to the 
threshold mass in this sector. Since its mass is lowered, it mixes more strongly with the qq 
sector in the ground state and from Fig. 6.5 we see that the ground state gets comparable 
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Forward-backward 






Symmetric 






(C 2 = 0.01, C 3 = 0.4) 


(C 2 = 0.1, C 3 = 0.4) 


K 




jet 














10 


0.38041 


0.4800 


0.4899 


0.5996 


0.3486 


0.4507 


0.4507 


0.5980 


18 


0.3722 


0.4968 


0.5110 


0.6447 


0.3402 


0.4673 


0.4673 


0.6409 


30 


0.3606 


0.5027 


0.5210 


0.6680 


0.3288 


0.4702 


0.4702 


0.6620 


42 


0.3511 


0.5029 


0.5240 


0.6765 


0.3189 


0.4677 


0.4677 


0.6682 


50 


0.3457 


0.5019 


0.5246 


0.6790 


0.3130 


0.4651 


0.4651 


0.6693 




0.3243 


0.5022 


0.5313 


0.6979 


0.2913 


0.4589 


0.4589 


0.6837 



Table 6.2: Lowest four eigenvalues (in units of G 2 ) in one link approximation. 

contribution from both sectors. The fifth state now corresponds to an almost free qq link state 
with infinitesimal qq component as shown in Fig. 6.5. 

6.6 Summary and Discussion 

We have performed an investigation of qq states using two different light front Hamiltonians 
in the one link approximation. The Hamiltonians correspond to two different ways of formu- 
lating fermions on the transverse lattice, namely, (a) forward and backward derivatives for 
and respectively or vice versa and (b) symmetric derivatives for both \]/ + and . In the 
latter, fermion doubling is present which is removed by an addition of the Wilson term. In this 
case there is no interference between helicity flip hop and helicity non-flip hop and, as a result, 
the qq component of the ground state wavefunction which has helicity plus or minus one are 
degenerate. In the former case, interference between helicity flip and helicity non-flip leads to 
the absence of degeneracy in the low lying spectra. One can recover approximate degeneracy 
of helicity plus or minus one components only by keeping the strength of the helicity non-flip 
hopping very small. 

Since the one link approximation is very crude and our motivation was to study and compare 
different fermions on the transverse lattice other than the assessment of the transverse lattice 
approach itself, we have not attempted a detailed fit to low lying states in the meson sector. 
Instead, we have explored the effects of various coupling strengths on the low lying spectra and 
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Figure 6.4: (a) Quark distribution function | y(x) | 2 of the ground state in the one link approx- 
imation, (b) qq contribution to the ground state, (c) qq link contribution to the ground state. 

(d) Quark distribution function | \]/(x) | 2 of the fifth eigenstate in the one link approximation, 

(e) qq contribution to the fifth eigenstate, (f) qq link contribution to the fifth eigenstate. The 
parameters are m f = 0.3, jl b = 0.2, C 2 = 0.4, C 3 = 0.01 and K = 30. 



associated wavef unctions. In this work, longitudinal dynamics is handled by DLCQ. We have 
performed a detailed study of various convergence issues in DLCQ using a wide range of K 
values. 

We summarize our results as follows. We have shown the effectiveness of appropriate coun- 
terterms in the qq and qq link sector to regulate the instantaneous fermion and fermion - link 
interactions respectively. We have also checked the cancellation of logarithmic divergences due 
to self energy effects. In the limit where fermions are frozen on the transverse lattice but un- 
dergo instantaneous longitudinal interaction, we have studied the convergence of ground state 
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Figure 6.5: Without the fermion - link instantaneous interaction: (a) Quark distribution function 
| | 2 of the ground state in the one link approximation, (b) qq contribution to the ground 
state, (c) qq link contribution to the ground state, (d) Quark distribution function | \\f(x) | 2 
of the fifth eigenstate in the one link approximation, (e) qq contribution to the fifth eigenstate 
multiplied by 10 4 , (f) qq link contribution to the fifth eigenstate. Parameters are the same as in 
Fig. 6.4. 



wavefunction with respect to K for three typical values of the fermion mass. We have stud- 
ied how the presence or absence of fermion - link instantaneous interaction in the qq link sector 
affects the wavefunction of low lying states. We have also studied the consequences of the inter- 
ference of helicity flip and helicity non-flip hopping in the Hamiltonian with forward-backward 
derivatives. This interference is absent in the symmetric derivative case. 
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CHAPTER 7 



Summary, Conclusions and Future Outlook 



Hadronic bound state problem is one of the most challenging tasks in nonperturbative QCD. 
Since it is very important to have direct access to the bound state wavefunctions which is essen- 
tial to calculate QCD observables, Hamiltonian approach is the most suitable candidate to ad- 
dress the problem of bound states. The simple structure of the vacuum in light-front QCD makes 
it possible to carry out a Hamiltonian analysis of the bound states in the Fock space language. 
In this thesis, we have investigated two nonpeturbative techniques in light-front QCD, namely, 
similarity renormalization group (SRG) approach and light-front transverse lattice (LFTL) ap- 
proach in the context of meson bound state problem. Unless and until we have complete control 
over the intricacies of the techniques in use, the dream to explore the nonperturbative QCD with 
full confidence will remain unrealized. In this work, we have performed a critical evaluation 
of the two nonperturbative approaches mentioned above and assessed them in terms of their 
strengths and weaknesses. 

We have first investigated Bloch effective Hamiltonian in the context of meson bound states 
in (2+1) dimensional QCD. There we encounter infrared divergences due to vanishing energy 
denominators in the bound state equation and SRG becomes mandatory to get rid of the diver- 
gences. The investigation of the qq bound states with Bloch effective Hamiltonian serves as a 
benchmark for comparative study of the same problem with SRG generated effective Hamilto- 
nian which is a modification over the Bloch Hamiltonian. Bound state equation in SRG scheme 
is free from the problem of uncanceled divergences coming from vanishing energy denomi- 
nators. We have also compared three different choices for similarity factor. To have better 
understanding, we have performed analytic calculation in the lowest order with step function 
similarity factor which in (2+1) dimensions generates linear confinement along the transverse 
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direction (for large x , V(x~,x l ) ~ x l ) while only square root confinement along the longitu- 
dinal direction (for large x~, V (x~,x l ) ~ Vx~ ) and thus breaks the rotational symmetry. Here 
one should recall that in (3+1) dimensions, in the lowest order, SRG generates logarithmic 
confining potential which also violates rotational symmetry. If the confinement generated by 
similarity transformation of the Hamiltonian is not an artifact of the lowest order approxima- 
tion, one might hope that the violation of the rotational symmetry will diminish with higher 
order corrections to the effective Hamiltonian. Higher order calculations are thus very impor- 
tant and illuminative in this context. It is also expected that higher order calculations in (2+1) 
dimensions will be much easier than in (3+1) dimensions due to simpler divergence structures 
and less demand of computing resources. 

We have studied another nonperturbative approach, the light-front transverse lattice formu- 
lation in this thesis. It is still a developing subject and is a very potential tool for nonperturbative 
investigations. Only very recently efforts are being given to formulate fermions and to study 
meson bound states on the LFTL. Lattice formulation of fermions is complicated due to gener- 
ation of extra species and needs special treatments to have a meaningful description of fermions 
in the discrete world. Thus it is very much important to know their origin and way(s) to over- 
come them before attempting any realistic QCD calculation on a lattice. We have shown that 
the origin of doublers on LFTL is completely different from usual lattice gauge theory. In the 
usual lattice gauge theory doublers come from the end of the Brillouin zone. But, when one 
uses symmetric lattice derivative to formulate fermions on a LFTL, doublers prop up due to 
decoupling of odd and even lattice points. We have also studied two different ways of removing 
doublers, namely, Wilson fermion and staggered fermion on LFTL. We have proposed another 
way of formulating fermions on LFTL by using forward and backward lattice derivatives in 
such a way that the Hermiticity of the Hamiltonian is preserved. In our method, there is no 
generation of extra fermion species. In this case, the helicity flip term proportional to fermion 
mass in full QCD becomes an irrelevant term in the free field limit. 

To assess which one of these two fermion formulations is better than the other we have 
compared them in the context of the meson bound state problem in (3+1) dimensional QCD. 
In the zero link approximation i.e., when q and q can sit only at the same transverse location, 
in the limit of large number of colors, it reduces to the 't Hooft model with linearly confining 
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instantaneous interaction along the longitudinal direction. In the one link approximation, q and 
q can be separated at most by one lattice point. With this approximation, the major difference 
between these two methods (in the case of symmetric derivative we add a Wilson term to re- 
move the doublers) is the interference between helicity flip and nonflip hopping in the case of 
Hamiltonian with forward and backward derivatives. As a consequence of this interference the 
degeneracy structure of the bound state spectrum of the Hamiltonian with forward and backward 
derivatives is different from that with symmetric lattice derivative. 

Since this subject is still under development, there are many unresolved questions one need 
to answer. The transverse dynamics and the structure of the mesons are too constrained by 
the one link approximation and investigations are essential with more than one links. By con- 
struction, quarks in one link approximation are confined in the transverse directions. The true 
nature of confining potential in the transverse directions can only be realized when sufficient 
number of links will be included. These issues can be investigated as an extension of our work. 
A systemic study of the effects of sea quarks also need to be undertaken. In this thesis, all the 
studies are done with light quarks, it will also be highly interesting to address the problem of 
mesons containing one light and one heavy quark in the context of heavy quark effective theory 
on the transverse lattice. 

A major unsettled issue in the transverse lattice formulation is the continuum limit of the 
theory when nonlinear link variables are replaced by linear link variables. The lack of well 
defined transformation rule between linear and nonlinear link variables makes it impossible to 
come back to nonlinear theory once the links are replaced by linear link variables and there 
is no straightforward way to take the continuum {a — > 0) limit. Thus, as an alternative on a 
coarse lattice, one needs to search for a trajectory in the parameter space with minimal viola- 
tion of Lorentz invariance. Another disadvantage of linearization of the links is that one needs 
to include more and more terms in the effective potential when more and more links are in- 
clude and calculations are viable only in the limit of large number of colors which suppress the 
higher order terms in the effective potential. It will be interesting to investigate the light-front 
quantization problem with non-linear constraints. The kinetic energy term of the nonlinear link 
variables is similar to the nonlinear o model. In this respect, the study of nonlinear o model on 
the light-front appears worthwhile. 
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APPENDIX A 



Notations and Conventions 



In this appendix we provide the notations and conventions used in different chapters of this 
thesis. 

For completeness let us start from the definition of light-front coordinates. The light-front 
coordinates in (3+1) dimensions are defined by 



X X I X ^ X 7 } 



(A.l) 



and any four vector in light-front is denoted as 



y^ = (y+,y-,y ± ). 



(A.2) 



The metric tensors are as follows 
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and the scalar product 



1 



1 



x-y=-x r y +-x y^-x -y 



The light- front partial derivatives are 



dx 



(A.3) 



(A.4) 



(A.5) 



(A.6) 



d = 2d+ = 2 



"dx+~ 
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(A.7) 



We define the integral operators 



F /co = \ I dy ~ e{x ~ - r)fiy ~^ (A - 8) 



(A.9) 



d 

f(x-) = l J dy- \x~ -y- | f(y~). 
The 7 matrices are define as 

7 ± = y°±7 3 . (A.10) 
We choose the representation for gamma matrices such that 



r + = 









-2/7 


( 

V 2il 


°) 

o • 


r- S 






d 1 = o 2 and d 2 = — c 1 where (7 r are Pauli matrices and I is a two component identity matrix. 
The projection operators 

A± = I r V = ^yV (A. 12) 



in explicit form are 



So, the fermionic fields 



V / ± =A ± V / (A. 14) 



V + =flO vr=f2 , l (A.15) 



can be written as 



where r\ and % are two component fields. 

In (2+1) dimensions the light-front coordinates are defined as 



x» = ( x ± = x°±x 2 , x l ) (A.16) 
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and accordingly the y matrices 

7 ± = y°±y 2 . (A.17) 
We use the two component representation for y matrices 

/ = - i ),r'= i .3 = (; y.r 1 = ;) <a,s, 



so that 



and 

A±-_ , _ , , _ _ , 

4 ; ' ' v ; ' v 1 y * 

Fermion field operator ty/^ = A ± We have 



A± = W, A + =fi °l A-=f2 ?l- (A.20) 



where £ and 77 are single component fields 



V + =(i), W-=(l) (A.21) 
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APPENDIX B 



Bloch Perturbation Theory for Effective Hamiltonian 

Here, we present the detailed formalism of the Bloch perturbation theory 1 to calculate the 
effective Hamiltonian. 

Consider a Hamiltonian H defined at a cutoff A. Let us try to lower the cutoff to X . In 
general, the cutoff could be in energy and/or particle number. Let us denote by Q the operator 
that projects on to all of the states removed when the cutoff is lowered. Let P = I — Q. We have 

Q? = Q, p2 = p p Q = Q p = (B.1) 

Our purpose is to find an effective Hamiltonian H e /f that produces the same eigenvalues in the 
sub space P as the original Hamiltonian H. 
Introduce an operator R that satisfies 

Q\y)=RP\y) (B.2) 

for all eigenstates of the Hamiltonian that have support in the subspace P. R gives the part of 
| outside the space projected by P in terms of the part of | iff) inside the space. Require that 
R gives zero acting on states outside the subspace. This means R = RP, R = QR, R 2 = 0. From 
R = QR, we have, PR = 0. Note also that i? f ^ R. 

Start from the set of equations (projections 2 of Schrodinger equation by P and Q) 

PHP | y)+PHQ \\jf) = EP | y), (B.3) 

QHP\ ¥ }+QHQ\ ¥ ) = EQ\ ¥ ). (B.4) 

1 C. Bloch, Nucl. Phys. 6, 329 (1958). Here we follow the treatment of R. J. Perry, Ann. Phys. 232 (1994) 116 
[hep-th/9402015] (reader can find here many examples of perturbative calculations); B. D. Jones and R. J. Perry, 
Phys. Rev. D 55, 7715 (1997). 

2 In Bloch-Horowitz formalism [C. Bloch, J. Horowitz, Nucl. Phys. 8, 91 (1958)] also, one has the same projec- 
tion operators P and Q, but not R. Substitution of Q \ y/) from Eq. (B.4) into Eq. (B.3) gives the Bloch-Horowitz 
effective Hamiltonian. 
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From Eq. (B.3), 

RPHP | iff) +RPHQRP | y) = ERP | y). (B.5) 

FromEq. (B.4), 

QHP | i//) + QHQRP | i//) = | iff) . (B.6) 

Subtracting, 

RHpp — HqqR + RHpqR — Hqp = 0. (B.7) 

We have introduced the notations, PHP = H P p and so on. Put H = h + v with [h, Q] = 0. Then 

Rhpp - h QQ R - vqp+Rvpp - vqqR+RvpqR = (B.8) 

which shows that R starts first order in v. 
We start from the eigenvalue equation, 

H(P + Q)\ ¥ )=E(P + Q)\xtr). (B.9) 

i.e., 

H(P+R)P | y) = E(P+R)P | iff). (B.10) 

Multiplying from the left by (P + R^) we have, 

(P + R J( )H(P + R)P | y) =E(P + R t )(P + R)P | i//). (B.ll) 

Using Pi? = 0, J?tp = 0, (P + P + ) (P + R)= P + R^R. Thus we can rewrite the eigenvalue equa- 
tion as 

1 



(P + tf)H(P + R) [y^#x\ '[1 +B*R]iP I V) 



= £[1 +R t R]?P | y/). (B.12) 



i.e. 



tf e// |0)=E|0) (B.13) 
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where 



0) = y) (B.14) 



and 



H, 



eff 



1 



(P + tf)H(P + R) 



1 



\ (B.15) 



Our next task is to generate a perturbative expansion. Denote free eigenstates in P by | a), 
| b), etc. Denote free eigenstates in Q by | /), | j), etc. Then 

&pp \a) =e a \a), 

h QQ \i) =£i\i). (B.16) 

Let us compute R to lowest orders in the perturbation theory. Let us write R = R\ +R2 + . . . 
where the subscript denotes orders in v. A straightforward calculation leads to 

</|*i|«> = (f|v f |fl) , (B.17) 

Vl 1 7 r (e a -£i)(e fc -ei) y (e a - £i )(£ a - £j ) 

Note that the energy denominators in the matrix elements of R involve only the difference 
between free energies of states in P and Q subspaces. This difference may approach to zero and 
give rise the "vanishing energy denominator problem" discussed in chapter 3. 

Our next task is to develop a perturbation theory expansion for the effective Hamiltonian to 
a given order. 

We start from the expression for the effective Hamiltonian (B.15). Remember that Ri ~ 
0(v),R 2 ~0(v 2 ). 

To order v, H e ff = PHP and hence 

(a\H eff \b) = {a\(h + v)\b). (B.19) 

To second order in v, we have 

H ef f = [1-^R][PHP + PHR + R^HP + R^HR][1-^R]. (B.20) 
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From R^HR we get 

In I R^HR I h\ = \ v 

7 (e a -e ; )(e fc -e;) 



From and i^/// 5 terms we get 

£(a | H | /)(/ 1 i?i | b) | /?J | \H\b) (B.22) 

= £ r (a|v|Q(/|v|^) + (a | v | | v | b) (B 23) 

Due to the normalization of the states, the effective Hamiltonian also gets contribution from the 
normalization factors. From the normalization factors we get 



Adding everything, to second order, we have, 



{a | H eff | b) = \YM I v NX* M b) 

1 i 



1 1 

■ + ■ 



(B.25) 



If a = b, this expression reduces to the familiar second order energy shift. 

Why Bloch formalism is preferred over B loch-Horowitz formalism? 

In the former, eigenstates of the effective Hamiltonian are orthonormalized projections of 
the original eigenstates. In the latter, they are not. The Bloch wavefunctions are defined by Eq. 
(B.14) and the Bloch-Horowitz wavefunctions by | 0) = P \ \\f). Consider two ortho normalized 
eigenstates of the original Hamiltonian | and | l/^) with | 1//2) = 0. However, P \ 
and P I i// 2 ) need not be orthogonal, i.e., (y/j \PP\ Wi) = (y\ \ P\ W2) ^ 0. Construct | \jr\) = 
[1 +R^R]?P I y/i), I \j/ 2 ) = [1 +R t R]?P I V2>- Then 

(xffi I Y2) = ( Vi I P I ¥2) + (Vi I PR*RP I W2) = (Yi \ (P+Q)\W2) = (Wi \ Y2}. (B.26) 
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APPENDIX C 



Details of Numerical Procedure to Diagonalize the Effective 

Hamiltonian 

We convert the bound state integral equation to a matrix eigenvalue equation by discretizing 
the integrations using Gauss quadrature points for both Bloch effective and similarity renormal- 
ized Hamiltonians. Here, we elaborate the parametrizations and diagonalization procedure we 
have used. 

Parametrization: The light- front variables are parametrized in the following ways in our 
numerical calculations. The full k- interval is divided into n\ quadrature points, k is defined by 
two different ways. One definition is 



where A is the ultraviolet cutoff and u's are the quadrature points lying between —1 and +1, so 
that k goes from —A to +A. The other definition is 



here k is a parameter that can be tuned to adjust the ultraviolet cutoff. The second definition 
(C.2) of k is very suitable for weak coupling calculations where we need maximum points to be 
concentrated near k = and get better convergence than the first definition (C.l). 

The longitudinal momentum fraction x ranges from to 1 . We divide all x- integrations in 
our calculations into two parts, x ranging from to 0.5 and x ranging from 0.5 to 1 and discretize 
each ^-interval into ri2 quadrature points with the parametrization 



k = 



uAm 



(C.l) 



(l-K 2 )A + m' 




(C.2) 



x — 



_ l + v + 2e(l-v) 

C ~ 4 
3+v-2e(l+v) 

4 ' 



0.5 <x< 1 -e, 



e <x< 0.5, 



(C.3) 



(C4) 
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where v's are the Gauss-quadrature points lying between —1 and +1 and e(— > 0) is introduced 
to handle end-point singularities in x as mentioned in the main text in Chapters 3 and 4. 

To handle the infrared diverging terms we put the cutoff \x — y\ > 8 and at the end we take 
the limit <5 ^ 0. Numerically, it means that the result should converge as one decreases <5 if 
there is no net infrared divergence in the theory. 

Diagonalization: After discretization, solving the integral equation becomes a matrix di- 
agonalization problem. The diagonalization has been performed by using the packed storage 
LAPACK 3 routines DSPEVX for the reduced model (real symmetric matrix) and ZHPEVX for 
the full Hamiltonian (Hermitian matrix). 



3 E. Anderson et ai, LAPACK Users' Guide, third edition (Society for Industrial and Applied Mathematics, 
Philadelphia, 1999). Available on the internet at the URL: http://www.netlib.org/lapack/lug/index.html. 
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APPENDIX D 



Nonrelativistic Bound State Equation 



In order to elucidate the implications of rotational symmetry in the (2+1) dimensional world, 
we review the nonrelativistic bound state equation in this appendix. 

For clarity, in this appendix we restore the superscript to the transverse component, namely, 
k = k l , q = q 1 , etc. We also use the notation k = (k l ,k 2 ), etc. To discuss the nonrelativistic 



limit of the reduced model defined in Sec. 3.6, make the variable change, x = ^1 + 



k 2 

E(k) 



y 



i 



i + 



1 



2 r ^ E(q) 

made. We have 



where E(p) = \J m 2 -\- (p 1 ) 2 + (p 2 ) 2 - So far, no approximations have been 

<k 2 \^ 



1 

x(l—x) 
m(x — y) 

ky — qx 

k(l-y)-q(l-x) 

dy_ 
dq 2 



1- 



E(k) 



2i 



m 



1„2 



2 \E(k) E(q) 
1 (q^ 2 



1 + 

ik 2 q 



m 

2 J2 



(q) E(k) 
\(k l q 2 q l k 2 
E{q)~m 
1 



(D.l) 



2 [E(q)]i 2E(qY 

The ~ equality holds in the nonrelativistic limit | k |, | q |<< ra. Introducing the binding energy 
B by M 2 = Am 2 (I —B), the bound state equation in momentum space in the non-relativistic 
limit is given by 



v(k) = 



An 2 m 



(D.2) 



Fourier transforming to coordinate space, with the momentum in the self energy integral cutoff 

by fermion mass m, one arrives at the coordinate space bound state equation 

f 
2n 
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i d 2 4/ 2 -i r ^ r i 

h -^—^ + 2-C s [Ye +lnmr]jy(r) = E yr(r) 



m 



dr 2 



Amr 2 



(D.3) 



where E = —mB and \j/(r) = r 2 i//(r)e ±!/( ^ and is the Euler constant. We note that rotational 
symmetry implies two-fold degeneracy for / ^ states. 
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APPENDIX E 



Similarity Renormalization Theory for the Effective 

Hamiltonian 



As it was promised in chapter 4, we present the detailed derivation of the effective Hamil- 
tonian using similarity renormalization group approach in this appendix. 

Since the renormalization group transformation based on integrating out the high energy 
states encounters nearly degenerate states, an alternative way of calculating effective Hamilto- 
nian was in demand. The solutions were proposed by Glazek and Wilson and Wegner indepen- 
dently 4 . 

Starting from a cutoff Hamiltonian Hb which includes canonical terms and counterterms we 
wish to arrive at an effective Hamiltonian H a defined at the scale o via a similarity transforma- 
tion 

H a =S a H B si (E.l) 

where S a is chosen to be unitary. 

The boundary condition is Limit H a = Hb- 

Introduce anti-Hermitian generator of infinitesimal changes of scale T a through 

S a = ST e J ° da ' T °' (E.l) 

where & puts operators in order of increasing scale. 

For infinitesimal change (lowering) of scale, S a = 1 — T a do and 5^ = 1 + T a do. Then we 
arrive at the infinitesimal form of the transformation 

^ = [H J \. (E.3) 

4 S. D. Glazek and K. G. Wilson, Phys. Rev. D 48, 5863 (1993); 49, 4214 (1994); F. Wegner, Ann. Phys. 
(Leipzig) 3, 77(1994). 
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This equation which has been called the flow equation of the Hamiltonian is the starting point 
of the investigations. 

The basic goal of the transformation S<y is that H a should be band diagonal relative to the 
scale a. Qualitatively this means that matrix elements of H G involving energy jumps much 
larger than o should be zero. T a still remains arbitrary to a great extent. It is instructive to go 
through the steps of the derivation which leads to the Glazek- Wilson choice. 

We write H B = H B o+Hbi where H B q is the free part and H B i is the interaction part of the bare 
cutoff Hamiltonian. A brute force way of achieving our goal is to define the matrix elements 
Hioij = faijHsiij where we have introduced the function f ai j = f(x C ij) with x a function of o 2 
and AM?-. The function f(x) should be chosen as follows: 

when o 2 » AM?-, f(x) = 1 (near diagonal region); 

when o 2 « AM?-, f(x) = (far off diagonal region); 

in between f(x) drops from 1 to (transition region) . (E.4) 

Here AM?-(= M? — M?) denotes the difference of invariant masses of states i and j. Because of 
the properties of /, Hi a jj is band diagonal. What is wrong with such a choice of inserting form 
factors by hand at the interaction vertices? First of all, we simply discard degrees of freedom 
above o. Secondly, H a will have very strong dependence on o. Thirdly, to ensure that H G has 
no ultraviolet cutoff dependence, H B should contain canonical and counterterms. But, in light 
front Hamiltonian field theory, because of the complexities due to renormalization, a priori we 
do not know the structure of counterterms. 

Note that in the definition of H a given in Eq. (E.3) the form of T a is still unspecified. In 
fact, a wide variety of choices are possible. In the following, we consider the choices made 
by Glazek and Wilson and Wegner. The price we have to pay for the use of flow equations is 
that it will generate complicated interactions even if the starting Hamiltonian has only simple 
interactions. For example, starting with a Hamiltonian which has only 2 particle interaction, the 
transformation will generate 3 particle interactions, 4 particle interactions, etc. 
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E.l Gtazek- Wilson Formalism 



Writing H a = Ho + Hj a , noting that the free Hamiltonian Ho does not depend on o and 
taking matrix elements in free particle states, we have, 



[HcJalij = {Pf -P;)Toij+[HloJo\ij 
where Hq \ i) = \ i), etc. . i.e., 

Jaij "<J Jaij Jaij 

Since we want Hiaij to be band diagonal, it is advantageous to trade j—. dH d '°' J for ^ 
which on integration has the chance to ensure that H I(y ij is band diagonal, we use 



(E.5) 



d 

do 



.ft 



-H 



aij 



'Iaij 



1 dfcjj 1 dH/gjj 

flij d ° ^ ~ foii do 



(E.6) 



fk Hloij . 



(E.7) 



and arrive at 



d 

do 



-Hioij 



aij 



faij 



(pr-p;)Tdj 



1 df ai 



Hioij 



(E.8) 



+ -F—[ H la,T a ]ij--M 

Jaij Jaij ~~ 

Still T ij is not defined. We next convert this equation into two equations, one defining the flow 
of Hiaij and other defining T Gi j. Recalling the starting equation Eq. (E.3) we add and subtract 
[Hi a , T a ]jj to the r.h.s. and arrive at 



d 

do 



.foij 



Hiaij 



= [H Ia ,T a ] lJ + ^-(Pr-P7)T (7 
Joij 



'.! 



+ J~M ~ frfi) ^ /(T ' T<y ^j ~ J2~ do i Hlaij ' 



(E.9) 



Glazek and Wilson choose T a to be 



Taij — ~ 



I 



laij pr-pr 
j i 



( 1 - faij) [H Ia , Ta]ij - ^ (In faij)H Iai j 



(E.10) 



Then from Eq. (E.9), we have, 



d_ 
do 



faij 



Hiaij 



— [Hia,T a ]ij. 
137 



(E.ll) 



Integrating Eq. (E. 1 1) from o to oo, we arrive at, 



Hloij — faij 



POO 

HiBij- / do'[H Ia i,T a i]ij 
J a 



(E.12) 



Note that Hj a ij is zero in the far off-diagonal region. This is clear from the solution given 
in Eq. (E.12) since f(x) vanishes when x > 2/3. 

T G ij vanishes in the near diagonal region. When i is close to j, fdj = 1 and both (1 — f G ij) 
and j^faij vanishes. It follows, then, from Eq. (E.10) that T C ij vanishes in the near-diagonal 
region. This guarantees that a perturbative solution to H/ ai j in terms of H B nj will never involve 
vanishing energy denominators. 

The effective Hamiltonian can be calculated up to any order of perturbation theory from 
Eq. (E.12) by iterative method. Here, we derive the effective Hamiltonian to second order in 
perturbation theory. Using 

H ioik- foikHeiik (E.13) 



and 



?*kj * {"^ (ln f<*j)f°*J H Bikj} (E.14) 



in Eq. (E.12), a straightforward calculation leads to 



H Iaij = ~Y< HBI ik H BIkj 



Soijk 8 a jik 



p k -PJ l\ I'i 



(E.15) 



where 



f°° d 
Saijk = faij J do' f a i ik ^—jfa'jki 

r i d 

8a jik =faij do fo'jkj^jfa'ik- (E.16) 

We find that the effective Hamiltonian in similarity perturbation theory is a modification of the 
effective Hamiltonian in Bloch perturbation theory 5 . 

5 Detail discussion of Bloch effective perturbation theory is given in Appendix B. 
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E.2 Wegner Formalism 

In the Wegner formalism 6 , the flow equation is given by 

dH{l) 



dl 

Wegner chooses 



[t(/),H(/)]. (E.17) 



x(l) = [H d ,H] = [H d ,H r ] (E.18) 

where H d is the diagonal part of the Hamiltonian and H r is the rest, i.e., H = H d + H r . Here the 
word diagonal is used in the particle number conserving sense. It is important to note that H d is 
not the free part of the Hamiltonian and both H d and H r depend on the length scale /. 

The light front Hamiltonian has dimension of {mass) 2 and hence x has the dimension of 
(mass) 4 , 1 has dimension of -, — ^r. 

v ' (massy 

Expanding in powers of the coupling constant, 

H = H^+H^ + H d 2) +h!: 2) + ... (E.19) 
where the superscript denotes the order in the coupling constant, 

x(l) = [Hf\H^\ + [Hf\H? ] ] + .... (E.20) 

Then, to second order, 

f = i[H?> ,H<Vf ] + Pf A (1 WV[[»f ,«?>],<] + ... ■ (E.2D 
Introduce the eigenstates of Hf\ 

Hf ) \i)=PT\i). (E.22) 

Then, to second order, 

^ = -(/T-/7) 2 < + ^ 1) ^ 1) ]^-(^-^) 2 < + --- • (E.23) 

6 For applications of Wegner formalism in condensed matter physics and quantum field theory see F. J. Wegner, 
Nucl. Phys. Proc. Suppl. 90, 141 (2000); E. L. Gubankova and F. Wegner, hep-th/9708054; E. L. Gubankova and 
F. Wegner, Phys. Rev. D 58, 025012 (1998) [hep-th/9710233]. 
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To first order in the coupling, 

d ^ = -(Pr-P J fH^ (E.24) 



(Pi-pj) 2 



which on integration yields 

H%](o) = e-^-H} l j{A) (E.25) 

where we have introduced the energy scale o via / = and used the fact that / = corresponds 
to the original bare cutoff. We notice the emergence of the similarity factor f a ij = e 



If we are interested only in particle number conserving (diagonal) part of the effective inter- 
action, to second order we have, 

d -^ = l4 l) M%j (E.26) 

Using 

$ = (Pi-PnH%, (E.27) 
the effective interaction generated to second order in the diagonal sector is 



h y h b h b vr-v+vr-v 



l _-{(Pi-Pk) 2 HPj-P^) 2 }l^ 



(E.28) 



r * KJ (pr-p k ) 2 +(p;-p k ) : 

Even though the second order formula is very similar to the one in Glazek- Wilson formalism 
when an exponential form is chosen for the similarity factor (see Sec. IV), we note a slight 
difference. In the Glazek- Wilson formalism, since the purpose is to bring the Hamiltonian 
into a band diagonal form, even in the particle number conserving sectors the large jumps in 
energies do not appear by construction. In the version of the Wegner formalism presented here 
the purpose is to bring the Hamiltonian in the block diagonal form in particle number sector so 
that large jumps in energies are allowed by the effective Hamiltonian. Note that small energy 
denominators do not appear in both formalisms. 
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APPENDIX F 



Violations of Hypercubic Symmetry on Transverse Lattice 



The canonical helicity non-flip interactions given in Eq. (6.23) for r ^ s break the hypercu- 
bic symmetry on the transverse lattice. For interacting theory this is also true for the Hamilto- 
nian with symmetric derivative. In the free field limit they do not survive for Hamiltonian with 
symmetric derivative but for forward-backward derivative they survive. In that case, in the free 
field limit they reduce to 

\ I dx ~mi n\* + at)d r d s ^-T}{n) 
+T7 t (x)d r d,-^T](x + as) 
-77 f (x + at) a r a s r\ (x + as) 
Going to the transverse momentum space via 

7](x-,x ± ) = j d 2 k ± e ikL - xL M-O 

we get 

J dx~ J d 2 k ± (j)l ± (x-)a3j^ T (j) k ±{x~) 
sin(^ y a) — sin(k x a) + sin(k x a — k y a) 



(F.l) 



(F.2) 



(F.3) 



Thus the violations of hypercubic symmetry are of the order of the lattice spacing a. Sign in 
front of this term changes if we switch forward and backward derivatives. 

In case of interacting theory with symmetric lattice derivative, these terms do not come in 
one link approximation but will come in if one considers more than one link. But in case of 
forward and backward derivative, they also appear in one link approximation. In our numerical 
studies presented in this work, we have set the coefficients of hypercubic symmetry violating 
terms to zero. 
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APPENDIX G 



Fermions With Forward-Backward Derivatives in 
Conventional Lattice Theory 



In chapter 5 we have discussed the light-front transverse lattice formulation of fermions 
with forward and backward lattice derivatives. Let us discuss the situation in conventional 
lattice gauge theory in this appendix. 

In discretizing the Dirac action in conventional lattice theory the use of forward or back- 
ward derivative for leads to non-hermitian action. The hermiticity can be preserved in the 
following way 1 . 

In the chiral representation 



o -/ 
-/ o 



o 



-a 1 



I 

-/ 



(G.l) 



The Dirac operator in Minkowski space 



<7% 







-io^d^ 



(G.2) 



where, = (7, a), = (/, —o). For massive Dirac fermions, this leads to the structure 

-ioPdnYR-myfL (G.3) 

-ia^dnYL-myfR. (G.4) 

For discretization we replace in Eq. (G.3) by forward derivative 

4 (8 y . x . fl 8 y . x )/a (G.5) 

and in Eq. (G.4) by backward derivative 

*i = (Sy*-Syw)/a. (G-6) 

7 In this appendix we follow the treatment of H. Banerjee and Asit K. De, Nucl. Phys. B (Proc. Suppl.) 53, 641 
(1997). 
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This leads to the structure 



z'y% - m = iy M A* - ifyftAj - 



(G.7) 



which results in hermitian action. Here, 




(G.8) 



Note that irrelevant helicity nonflip second order derivative term is produced in this method of 
discretization. In contrast, the corresponding term in the transverse lattice depends linearly on 
m and flips helicity. One can trace this difference to the presence of the constraint equation in 
the light front theory. 

Writing lattice derivatives in this fashion in the conventional lattice theory eliminates the 
doublers from the edges of the Brillouin zone but some non-covariant doublers prop up from 
other places. 
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APPENDIX H 



Mixed Covariant Derivative 



Here we derive the form of mixed covariant derivative used in the transverse lattice formu- 
lation. 

In a local field theory, we cannot simply compare objects at distances. As Feynman reminds 
us, "We must take into account of the rotation of the frame by transporting U (x! 1 +Ax^) back 
to before making a comparison'^ . Hence the total change is given by 



8 R. P. Feynman, in Weak and Electromagnetic Interactions at High Energies, (eds.) Roger Balian, Christopher 
H. Llewellyn Smith, (Elsevier North-Holland Pub. Co., Amsterdam 1977). 



D^U r {x^ ,x)Ajc j 



R{xy l U r (x^ +Ax^,x)R(x + ar) -U r (x^,x) 

[l+AP (x) Ajc m ] U r (x^ + AxV , x) [1 - A» (x + ar) A* M ] - U r (xP , x) 

C/ r (x M ,x) + ^C/ r (x^,x)Ajc M -C/ r (x M ,x) 

(x) U r (x^, x) Ajc m - U r (x? , x)A^ (x + ar)Ax^ 




(H.l) 
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APPENDIX I 



Transverse Gauge Invariance 

Let us explicitly verify the transverse gauge invariance of the lattice theory. The theory is 
invariant under the gauge transformations 

r 7 (x)^r 7 / (x)=G t (x)r 7 (x) (1.1) 

and 

Af r (x) -> A#(x) = G t (x)M r (x)G(x + ar) (1.2) 

where 

G(x)=e- iTa0a M. (1.3) 

For infinitesimal transformation, 

G(x) « l-iT fl a (x) (1.4) 

and 

77 (x) -> T7 / (x) = r] (x) + iT a e a (x)ri (x) (1.5) 

and 

M r (x) M ^M;(x) M =M r (x) M + /r; / M r (x) / ,0 fl (x) -/M r (x) p/ 7^e a (x + ar). (1.6) 

In quantum theory the gauge transformations are generated by the operator 

^ = e iLyQ a (y)0 a (y) (L7) 
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with 



Q a (y) = jdy- rr{r fl £(M^y)/r+Mj,(y)+M;(y-arO/r+M ; .,(y-^ / ))} 



-2n\y)T a 7]{y) 



so that 



and 



tj(x) — > 7]'(x) = Sfr](x 



M r (x) ->M£(x) = ^M r (x)Sf. 



(1.8) 



(1.9) 



(1. 10) 



For infinitesimal G a , using [A,5C] = A{5,C} — {A,C}5 and the canonical commutation 
relations for the fermion field operator, we readily verify Eq. (1.5). Using [A,5C] =A[B,C] + 
[A,C]fi and canonical commutation relation for the link field, we also verify Eq. (1.6). 

Next we look at the behavior of fermion and link creation and annihilation operators under 
transverse gauge transformations in order to construct gauge invariant multiparticle states. 

The gauge transformation on the link variable is 



M r (x) -^M' r (x) = G + (x)M r (x)G(x + af) 
where G belongs to SU(N). Thus 

A r (x) -^A' r (x) = G + (x)A r (x)G(x + af), fij(x) -> B^' r (x) = G t (x)5j(x)G(x + af). 



(1.11) 



(1.12) 



Thus Tr(A^B^), Tr(AB), 7V(AtA), Tr(B^B) are locally gauge invariant operators. A locally 
gauge invariant two link state is Tr(A^B^) | 0). 
Recall that the current has the structure 



j+ a ~ Tr(T a MM^) with M ~ A + 5 t . 
The interaction term J +a (-j T ) 2 J +a has many terms. Consider one term 



(1.13) 



Tr(T a AA^) Tr{T a B*B) = Tr({A^T a A) Tr{BT a B^) = T^ p T s a t Al m A pn B rs Bt 

T« p T? t 8 nr 8 mt 8 ps 8 nr =(T a T a ) mm . (1.14) 
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On the other hand, consider the "pair creation term" 

Tr(T a AB)Tr(T a B^A f ) = Tr(BT a A)Tr(A f T a B^) 

— 1 mp 1 st D nm^pn^rs D tr 

T« p T*8 mp 8 ts = 0. (1.15) 

Thus pair creation or pair destruction terms (even if they conserve particle number) do not 
contribute if we restrict ourselves to a two link gauge invariant sector. 

Next let us look at the two component fermion field r\ (x) . The transformation of r\ (x) is 

r 7 (x)^r 7 (x) / = G t (x)r 7 (x). (1.16) 

Since rj(x) m b(x) + d^(x), we have 

b(x) ^b(xY = G f (x)b(x) and d^(x) -> d*(x)' = G t (x)d' t (x). (1.17) 

Next consider how to form gauge invariant two-particle (qq) and three-particle (qq link) states. 
A gauge invariant qq state is b^(x)d^(x) | 0). Gauge invariant three particle states are 
& t (x)fi+(x)J t (x + of) | 0) andM(x + ar)Aj(x)J t (x) | 0) (= b^(x + ar)Bl r (x + ar)d^(x) | 0)). 
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APPENDIX J 



Hamiltonian Matrix Elements in One Link Approximation 



J.l Structure of terms in DLCQ 

We use DLCQ for the longitudinal dimension (— L <x < +L) and implement anti periodic 
boundary condition for the two component fermion field, 

7] c (*-,x) = 4fEj& £ [b c (m,xA)e- i ^^+4(m,x,-X)e 1 ^/^] (J.l) 

V lL X m=l,3,5,... 

with 

{b c (m,x,X),bl(m',x',X')} = {J c (m,x, A),^(m',x', A')} = 5 mm /5 x y 5 C;C /5 A;A /. (J.2) 
The link field has periodic boundary condition (with the omission of the zero momentum mode), 
M rpq (x-,x) = ^= £ -^[5_ rp9 (m,x + af) e -^-/ L + Bt ( m5X ) e ^-A)] (j. 3 ) 

V™ m= L2,3,...V m 

with 

M (m, x) , B l fs (m', x') ] = 8 mn j 8 X ^ 8 r y 8 ps 8 qt . (J.4) 
The Hamiltonian P =±H. 

In the following subsection we give the explicit structure of terms in the Hamiltonian in the 
forward-backward case in DLCQ restricting to those relevant for the one link approximation. 

Mass terms 

Mass terms: 

H f free = "^Y.Hllll' b c ("> x > *)M W > x, X) + 4(n,x, X) d c (n,x, A)] . (J.5) 



x c X n 
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HUNK free 



J" 



,EEEj 



5 r (m, x)B r (m, x) + B_ r (m, x + ar)B- r (m, x + ar) 



(J.6) 



Four fermion instantaneous term 



The four fermion instantaneous term which gives rise to a linear potential in the color singlet 
state 

g 2 

2 ^2 E E E 5 AA'$i"A'" E 

x Z?J(mi,x, A)<iJ,„(m4,x, —X"')b c i (rti2,x, X')d c ii (m?, ,x, X'") 

1 



X 



(m3 — m4) 2 



m 2 +m 3 



(J.7) 



Helicity flip terms 



Particle number conserving terms: 



mg 

(2 



EE E 4 £ 



t/ 



mi 



mi 



^ (mi , x, X\)b c (m\ , x, X 2 ) + dj. (mi , x, -A 2 )d c (mi , x, -A 2 ) 



Particle number non conserving terms: a typical term is 



mg 1 



EE E A 6rX ^ E 



i 



i 



a r x 1 1; A 2 



m\m2mj, 



m.2, 2m3 + m 2 



3m!— ni2,2m 3 



b\ (mi , x, Ai )B_ rcc i (m3 , x + ar)b c i (m 2 , x + ot, A 2 ) . 



(J.8) 



(J.9) 



Helicity non flip terms 



Two operators: 



2 1 r 

~EEE~ ^ t («,x,A)^ c .(n,x,A)+4(«,x,A)J c (n,x,A 



a x X n H 



(J. 10) 



Three operators: 
A typical term is 



1 1 



^^/z^EEE E 

" Vtl r x Jl minion 



1 1 



l_i l_i l_i z_ i '^m l -m 2 ,2m 3 

r x A m»3 V^3 2m 3+ m 2 

b\{m\ 1 x, X)B_ rc( j (m-i , x + af )Zv (m 2 , x + ar, A ) . 



(J. 11) 
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Fermion - link instantaneous term 

A typical term is 



2 L fl 2EE E Yjcd T dd< E r— r— 

* !la x rcdd'ddl m l m 2 m 3 m 4 V m 3 V m 4 

b\{m\ , x, Ai )b d > (m 2 , x, A 2 )B_ rc > c « (m 3 , x + ar)Bt rc „ c (ra 4 ,x + of) 

(-)(m 3 + m 4 )/(mi -m 2 ) 2 <5 mi _ m2i 2m 3 -2m 4 • (J-12) 



J.2 States in DLCQ 



We will consider states of zero transverse momentum. In the one - link approximation, the 
gauge invariant states are qq state 

1 1 

fcj(ni,y(«),ai)^(n 2 ,y(9),^) | 0) (J. 13) 

and the qq link states 

111 

' 3 ^ = wlr^III I I 5 y(0,y(?) fy?),y(«)-flS 
iV V2 ^ , y((?) y(§) y(z) 



Wi>y(«)>tfi) fi L d '("3,y(0) dl(n 2 ,y(q),oz) 1 0) 



and 



' 3 ^ = i^^III I I 5 y(/),y(?) 5 y(?),y(?)+«s 

^ , y (9 ) y (§) y (/) 

^(ni,y(9),CTi) <„ w ,(n 3 ,y(/)) 4(n 2 ,y(?),(T 2 ) | 0). (J.14) 
We shall consider transition from these initial states to the following final states: The qq state 

(0 I c/ e (« 2 ,z(<7),(7 2 ) fc e (ni,z(^),oi) (J.15) 

and the <7^ link states 

^ V V ee > t z ( 9 ) z (^) z (/) 

(0 I d e {n 2 ,z(q),o' 2 ) B tee/ (n 3 ,z(l)) b e/ (n\,z(q),o[) (J.16) 
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and 

111 

yv yi ee , t z ^ z ^ z ^ 

(0 | d e (n 2 ,z(q),<T2) B_ tee ,{n'^x{l)) b e >(n[,z(q),a[) 

(J. 17) 

J.3 Forward-backward derivatives: Matrix Elements in DLCQ 

J.3.1 Transitions from two particle state 
To two particle state 

Let us consider transitions to the two particle state: We have, from the free particle term, 

(2 f | Hf free \2)=m 2 (— + —) (J. 18) 



where 



^1 = $n u n' Vtf 8 n 2 ,n< 8 a 2 ,a' • ( J - 19 ) 



From the four fermion instantaneous term we get 



g 2 



(2' | Hqq C | 2) - 2^ 2 C f 5 m+niA+n , 2 ^ _ w , )2 

5^, a; 5 CT2; ^ (J.20) 



where C/ = 4\T ■ 

To implement the regulator prescription for we add the counterterm matrix elements 

g 2 £ 1 

(2' I H CT I 2) = 2^ C/ S MlWi+n , {ni - nioop) 2 8 °v4 5 ^ 2 - (J-21) 

Here the term n; oop = «i is dropped from the sum. 
From the helicity flip term we get 

(2' \H hfl I 2) = -2^£ ^xl^sXa^^ + ^xl^sX-d,^ ^hf 

(J.22) 

with 

«^/ = ^m,^ 5 n 2 ,4 • (J- 23 ) 
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From the helicity non-flip term we get 



(2' | H hnf {\) | 2) = 2\(- + -) ,y^ 2 ■ (J-24) 

a 1 \n\ n 2 J 



To three particle state 
To the state | 3a) 

From the helicity flip term we get 



d «2,«2 n' 7^ L L d z(q),y(q)-ai 

H l V H 3 z(q) y(q) 



mg 
a 



8n' 2 +2n' v n 2 1 V V ^ / T TO 
ni ' n 'i ^ 7^ ^- ^- d z®,y(?)+«t- ( J -^ 

" 2 V n 3 z(§) y(?) 



From the helicity non-flip term we get 

(3a' | H hnf (2) | 2) = -^^^-^^.j^ 

„ ^n' 1 +2n' 3 ,«i 1 V V V s; 
"l V"3 / z( 9 ) y( 9 ) 
1 /i? 1 1 1 X X 



" 2 V"3 * z(«) y(§) 



To the state | 3b) 



From the helicity flip term we get 

mg /- 1 1 I 



(3^|^ /2 |2) = 

" 2 '"2 m 7^ 2- L °z(q),y(q)+ai 

V"3 z( 9 ) y( 9 ) 

fflf / — 11 1 + ^ 



Sn' 2 +2n' v n 2 1 V V £ rT ^ 

^ 7=f 2- L °z(q),y(q)-ai ■ 

H 2 V H 3 z(q) y(q) 

152 



From helicity non-flip term we get 

(3b' | H hnf (3) | 2) = -g^ 2 VN l -^^=8 aua[ 8 a2 ^ 



^+2n^ni 1 V V V £ 

"2,4 „, /7L L L °z( 9 ),y( 9 )+at 

S ni ,n\ " 2 ^" 3 '" 2 -Ff £ £ £ 5 z( ? -),y(g)- fl t • (J- 28 ) 

"2 V"3 * z(q) y(q) 



J.3.2 Transitions from three particle (q q link) state | 3a) 
To three particle state 

From the free particle term, we get 



{3a' | H free \ 3d) = I m 2 ( — + —) + l -n 2 - ] ^ (J-29) 



with 



^ = 5n u n\ S n 2 ,n> 2 ^ 3 X 3 5 <7l,cr{ ^cr 2 ,^ • (J-30) 

Diagonal contribution from the four fermion instantaneous term to the three particle state 
vanishes due to the vanishing trace of the generators of SU (N). 
Contribution from the fermion - link instantaneous term 

g 2 1 

(3d | H qgc (\) | 3a) = - — -^Cf <5„ 1+2 „ 3 / i+2 „' 3 8 n2 ,n' 2 

1 (ni-n[+4n 3 ) J_ 

„2 l 



K a 2 C f 8n 2 +2n 3 ,ri 2 +2ri 3 ^n u n[ 



1 (n 2 -n' 2 +4n 3 ) 1 



(J.31) 

Counterterm matrix elements in DLCQ to implement the regulated prescription for jj^y 

g 2 1 

(3a | #cr(l) I 3a ) = — ^2 C / ^ 1 +2« 3 ,«' 1 +2«' 3 <5„ 2 / 2 5 CTlj0 j 5^ 
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E 

wp — 

+ E 

n loop~ 



1 



(ni -ni oop + 4n 3 ) 1 



,,,.,„, -i V^3\/ n i -nioop + ln?, (ni-ni oop ) 2 y/2 

nimax 1 {ni-n loop + Anj) 1 

l V"3\/ n 2 ~ nioop + 2n 3 (n 2 -ni oop ) 2 y/2 



(J.32) 



where n lmax < n\ + 2n 3 and n 2max <n 2 + 2n 3 . 

The contribution from the helicity flip term that conserves particle number is 



{3a' | H hf {\) | 3a) = -2 — 8^ 8 ni ^ 8 n ^ 

- E xl Zr Xa x 8o 2 ,oi + E X-cM-^ 

"■If 1 ,l 2 y 

The contribution from the helicity non-flip term that conserves particle number is 



(3a'\H hnf (\) \3a) = + 

a 1 \n\ n 2 J 



To two particle state 



From the helicity flip term we get 

mg 



1 1 1 



^n\,ni+2n 3 1 



J n 2 ,n' 2 



E E 5 <q), 



y(?)+as 



1 1 1 



& 2 
U «l,n'l 



E E 5 S 

From the helicity non-flip term we get 

(2' | H hnf | 3a) = -gl VN 1 1 1 



z(?),y(?)-«s 



V V2 ^ ^ 
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y(q)+as 



(J.33) 



(J.34) 



(J.35) 



o ^n' 2 ,n 2 +2n } 1 Y" 1 V £ 
0„, „' ; — 2^ L d z( ? -),y(g)-as 



" 3 y(?) 



J.3.3 Transitions from three particle (q q link) state | 3b) 
To three particle state 

From the free particle term, we get 



(3b' | H free | 3b) = (m 2 (- + -)+ U 2 -) 
\ \ni ri2/ 2 n 3 ] 



with 



(J.36) 



(J.37) 



(J.38) 



The diagonal contribution from the four fermion instantaneous term to the three particle 
state vanishes due to the vanishing trace of the generators of SU (N). 
The contribution from the fermion - link instantaneous term is 



(3b'\H qgc {\)\3b) 



n 



2 1 

^2 C f ^ni+2n 3 ,n[+2n' 3 $n 2 ,n 2 

1 (ni-n[+4ri3) 1 



S 2 1 



7T a 



2^f 8?i 2 +2n 3 ,n 2 +2n' 3 ^n\,n\ 



1 



^55- ^ V« ^ (J.39) 



^/n 3 ~^n 2 -n' 2 + 2n 3 (n 2 

Here also we have the counterterm matrix elements given in Eq. (J. 32). 

The contribution from the helicity flip term that conserves particle number is 



(3b'\H hf {\)\3b) 



a U "'l w 2,«2 »3,« 3 



Hi r 

The contribution from the helicity non-flip term that conserves particle number is 

2/1 1 \ 



.(J.40) 



(3b'\H hnf (l) I 3b) 



(J.41) 
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To the two particle state 



From the helicity flip term we get 



« ^n' 1 ,ni+2n3 1 V V ^ 
"1 V"3 z((?) y(?) 
tn 1 1 1 



V y/2 y^4n 

^n' 2 ,2n 3 +n 2 1 

From the helicity non-flip term we get 



^' 1 ,wi+2ra 3 1 V V V ^> 
"1 V"3 s z ( a ) v( a ) 



^n' v n\+2riT, 1 

<q) y(q) 

Sn u n[ " 2 '" 2+2 " 3 — £ £ £ 5 z(9),y( 9 )+as • (J-43) 

J.4 Symmetric derivatives and Wilson term: Matrix elements in DLCQ 

In this section, we list only those matrix elements that differ from the forward-backward 
case. 

J.4.1 Transitions from the two particle state 
To the state | 3a) 

Helicity flip: 

< 3a < i P- bf 1 2) = („ + £ ) 1 Vn I -L J= £ d< ^ s<m5 

£ £ 5 z( ? ),y(g)-«t 

y(«) 

'I n-, n . n 



H2,n 2 n l +2n 3 ,n 1 
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Helicity non-flip: 



E E S z(q),y(q)+d 

y(<?) Aq) 

1 a44) 



K\ K 1 1 1 

AT = —=o i 5 



E E H 5 z( q ),y( q )-ai 

1 y{q) z («) 



1 



«1 «j / " 2 '"2 «i+2n 3 ,«l 



E E E 5 z®,y(^)+at 

f y(<?) 

1 ( 1 l \ c 

— r H o / o i . - / 

n' n 9 / "i- n i "2+ 2 "3-"2 



"3 v * 2 



(J.45) 



To the state | 3b) 



Helicity flip: 



E E * 



y{q) <q) 
1 



z(q),y(q)+at 



«j n 'j / »2,n 2 n 1 +2n 3 ,n, 



a/ 2a VV2V5t 

E E 5 z( ? -),y( ? -)-at 



1 



n' n 2 / Wl ' n i "2+ 2 «3'«2 



(J.46) 
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Helicity non-flip: 

(3a'\p- n \2) = -(m + 4-)-VN-^=-^8 ,8 , 

E E E 5 z(q),y(q)+ai 



' I "1 «' / " 2 '" 2 n l+ 2w 3' n l 

n 3 \ 1/ 

- + - v 7 ^— -4= — ; =<5 ' <5 / 
V a) a V y/2 a >' a i 

E E E 5 z(q),y(q)-ai 

1 ^ 1 O s 

-r H o i o i.r, i 

' ni J n u n x n 2 +2n 3 ,n 2 



n' \' l 2 



J.4.2 Transitions from three particle state | 3a) to two particle state 

Helicity flip: 

(2'|P; t/ |3«> = K^^^I^te,^; 

EE 

z(g) z(g) 

4-- I 5 ' 5 



>/2 V4tt' 

E E <V?),y (<?)-«§) 
z(?)y(?) 

'I 17 . m 17- - 



\7«3 l«2 Wo/ ni ' Bl " 2+2 " 3 '" 2 ' 



Helicity non-flip: 



(2'|p- f1 |3a) = -fm + 4- > )-v / iV-4= ' 5 ' 

E EE 

- s ' y(?) 

1 / 1 1 

+ - 5 , S 



V«3 I ni n'J n2 '"2 »i+2»3.»i 
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(m + 4-) -VN^--^= —^5^ td n > 
\ a) a V y/2 \[4% ° l > a i 

E E Sz(q),y(q)+as) 

' - I Hi n 



yn' 2 ^2 / n l' n l n 2 +2n 3 ,n 2 

J.4.3 Transitions from three particle state | 3b) to two particle state 

Helicity flip: 

E E 5 z( 9 ),y( 9 )-«s) 

y(?) 



EE 

z(?)y(?) 

1 / 1 1 



5 / 5 



Helicity non-flip: 



(2'|p- fc -|3fc) = -(m + 4-)-VN--^=-^S ,8 , 

EEL ^y^-as) 

- s ' z(<?) y(<?) 

1 + J r) 5 „: 5, 



v /njln 1 w ' / "2,n 2 ni-t-2n 3 ,n, 
/ A K\ K i— 1 1 1 

V a/ a ^ ai,<Tl 

E E E S z(q),y(q)+as) 

s y(?) 



1/1 1 

- + - 5 , 5 



V^3 \ M 2 ^ 2 / Wl ' ni "2+ 2n 3' W 2 
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J.5 Self energy counterterms 



In this section we list the self energy counterterms. 
Symmetric derivatives case 

The counterterm for self energy for a quark or an antiquark with longitudinal momentum n\ 
due to double helicity flip hops 

CT\ = — V — — = — (J.52) 

n\ n\ n z n\n\ +m z {ni — n\) z 

The counterterm for self energy for a quark or an antiquark with longitudinal momentum n\ due 
to double helicity non-flip hops 

CTl = Y" — ~ ^— 7T7- T\~l ' (J. 53) 

n\ f Jr± 1 n\ \l L n\n x + m z (ni —n\) A 

Forward and backward derivative case 

In this case we have three types of contributions: (1) helicity flip acting twice, (2) helicity 
non-flip acting twice and (3) interference of helicity flip and helicity non-flip hops. The first 
two are diagonal in helicity space but the last one is off-diagonal in helicity space. 

The transition from state | 2) to state | 3d) and back due to a quark hop gives rise to longitu- 
dinal infrared divergence. In this case the counterterm due to double helicity flip is 

CT 3 = lY -r—. -. \r-. j-,. (J.54) 

n \ li n\n\ +m z (n\ — n\) A 

The counterterm due to double helicity non-flip is the same without the factor of 2. The transi- 
tion from state | 2) to state | 3b) and back due to a quark hop does not give rise to longitudinal 
infrared divergence. Similarly the transition from state | 2) to state | 3d) and back due to an 
antiquark hop does not give rise to longitudinal infrared divergence. The transition from state 
| 2) to state | 3b) and back due to an antiquark hop gives rise to longitudinal infrared divergence 
which requires counterterms the explicit forms of which are the same as in the quark case for the 
transition from | 2) to state | 3d). Lastly we consider counterterms for self energy contributions 
arising from the interference of helicity flip and helicity non-flip hopping. The counterterms 
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have the same structure as in the case of helicity non-flip transitions accompanied by the fol- 
lowing extra factors. Since we have two possibilities namely helicity flip followed by helicity 
non-flip and vice versa and these two contributions are the same, we get a factor of two. We 
also get a factor Xj&^Xs where s(s') is the initial (final) helicity and d 1 = o 2 , 6 2 = —a 1 . 
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